I was able to compile the SPheno output for the B-L-SM Model of Sarah in both Mathematica 11 and 13.
However, when I change the B-L charge assignment with a new anomaly free set which involves a non-zero B-L charge for the SM Higgs (attached), I am still able to compile the SPheno output in Mathematica 11, but not for Mathematica 13.2. In fact:
1) I get the following error during the MakeSPheno[] part:
Code: Select all
Transpose::nmtx:
The first two levels of {{{Ah[{gt1}], Ah[{gt2}], VP[{lt3}], VP[{lt4}]},
I 2
{-- (((4 g1 + g1p1) Cos[ThetaW] - 4 g2 Sin[ThetaW]) ZA[gt1, 1]
32
2 2
ZA[gt2, 1] + 256 g1p1 Cos[ThetaW] ZA[<<2>>] ZA[gt2, 2]), <<1>>}},
<<32>>} cannot be transposed.
Transpose::nmtx:
The first two levels of {{{Ah[{gt1}], Ah[{gt2}], VP[{lt3}], VP[{lt4}]},
I 2
{-- (((4 g1 + g1p1) Cos[ThetaW] - 4 g2 Sin[ThetaW]) ZA[gt1, 1]
32
2 2
ZA[gt2, 1] + 256 g1p1 Cos[ThetaW] ZA[<<2>>] ZA[gt2, 2]), <<1>>}},
<<32>>} cannot be transposed.
Transpose::nmtx:
The first two levels of {{{Ah[{gt1}], Ah[{gt2}], VP[{lt3}], VP[{lt4}]},
I 2
{-- (((4 g1 + g1p1) Cos[ThetaW] - 4 g2 Sin[ThetaW]) ZA[gt1, 1]
32
2 2
ZA[gt2, 1] + 256 g1p1 Cos[ThetaW] ZA[<<2>>] ZA[gt2, 2]), <<1>>}},
<<32>>} cannot be transposed.
General::stop: Further output of Transpose::nmtx
will be suppressed during this calculation.
Part::partw: Part 2 of Transpose[{{{Ah[{gt1}], Ah[{gt2}], VP[{lt3}],
I 2
VP[{lt4}]}, {-- (((4 g1 + g1p1) Cos[ThetaW] - 4 g2 Sin[ThetaW])
32
2
ZA[gt1, 1] ZA[gt2, 1] + 256 g1p1 <<2>> ZA[gt2, 2]), g[<<2>>]}},
<<32>>}] does not exist.
Part::partw: Part 2 of Transpose[{{{Ah[{gt1}], Ah[{gt2}], VP[{lt3}],
I 2
VP[{lt4}]}, {-- (((4 g1 + g1p1) Cos[ThetaW] - 4 g2 Sin[ThetaW])
32
2
ZA[gt1, 1] ZA[gt2, 1] + 256 g1p1 <<2>> ZA[gt2, 2]), g[<<2>>]}},
<<32>>}] does not exist.
Part::partw: Part 2 of Transpose[{{{hh[{gt1}], VP[{lt2}], VZ[{lt3}]},
{<<2>>}, {(-8 I) g1p1 vX Cos[ThetaW]
(g1p1 Cos[ThetaWp] Sin[ThetaW] - g1p Sin[ThetaWp]) ZH[gt1, 2],
g[lt2, lt3]}}, {<<2>>}, <<9>>, {<<2>>}}] does not exist.
General::stop: Further output of Part::partw
will be suppressed during this calculation.Code: Select all
CouplingsForDecays_BLSM.f90:818:37:
818 | rHB_S_VZ(i1) = Abs(-0.5_dp*cplhhVZVZ(i1)*vev/MVZ2)
| 1
Error: Syntax error in argument list at (1)
CouplingsForDecays_BLSM.f90:819:41:
819 | rHB_S_VWm(i1) = Abs(-0.5_dp*cplhhcVWmVWm(i1)*vev/MVWm2)
| 1
Error: Syntax error in argument list at (1)
CouplingsForDecays_BLSM.f90:810:29:
810 | cplHiggsWWvirt = cplhhcVWmVWm/vev
| 1
Error: Symbol ‘cplhhcvwmvwm’ at (1) has no IMPLICIT type
CouplingsForDecays_BLSM.f90:811:26:
811 | cplHiggsZZvirt = cplhhVZVZ*Sqrt(7._dp/12._dp-10._dp/9._dp*Sin(TW)**2+40._dp/27._dp*Sin(TW)**4)/vev
| 1
Error: Symbol ‘cplhhvzvz’ at (1) has no IMPLICIT type; did you mean ‘cplhhhhhh’?