Page 1 of 1

SPheno compilation errors with unmixed fermions

Posted: 19. Apr 2016, 10:49
by ForumImport
Thread imported from first SARAH forum


Hi Florian,

I have encountered a problem in generating the SPheno output for models with multiple unmixed fermions in SARAH-4.8.2, which seems to be a bug in the most recent version? For example, the generated SPheno code for the NMSSM+VL/10plets model fails to compile, with the errors:

Code: Select all

SusyMasses_NMSSM10.f90:56.16:

Complex(dp) :: 0.C 
               1
Error: Invalid character in name at (1)
SusyMasses_NMSSM10.f90:229.16:

Complex(dp) :: 0.C 
               1
Error: Invalid character in name at (1)
SusyMasses_NMSSM10.f90:268.5:

MFepC = (sqrt(2._dp)*vS*Conjg(le10) + 2*Conjg(mve))/2._dp 
    1
Error: Symbol 'mfepc' at (1) has no IMPLICIT type
SusyMasses_NMSSM10.f90:95.5:

MFepC = (sqrt(2._dp)*vS*Conjg(le10) + 2*Conjg(mve))/2._dp 
    1
Error: Symbol 'mfepc' at (1) has no IMPLICIT type
make[2]: *** [../lib/libSPhenoNMSSM10.a(SusyMasses_NMSSM10.o)] Error 1
This seems to be because the left-handed neutrinos are being included in the the list of unmixed fermions that are used for generating the necessary type declarations. That is, in SPhenoTreeMasses.m, should the lines

Code: Select all

nmixF=Select[listNotMixedMasses,getType[#[[1]]]===F&&#[[2]]=!=0&];
For[i=1,i<=Length[nmixF],
WriteString[sphenoTree,"Complex(dp) :: "<>ToString[SPhenoMass[listNotMixedMasses[[i,1]]]]<>"C \n"];
i++];
in WriteCalcAllTreeMasses instead be

Code: Select all

nmixF=Select[listNotMixedMasses,getType[#[[1]]]===F&&#[[2]]=!=0&];
For[i=1,i<=Length[nmixF],
WriteString[sphenoTree,"Complex(dp) :: "<>ToString[SPhenoMass[nmixF[[i,1]]]]<>"C \n"];
i++];
? Making this change does then mean that the generated SPheno code successfully compiles, at least.

Thanks very much for your help!

Cheers,
Dylan

Re: SPheno compilation errors with unmixed fermions

Posted: 19. Apr 2016, 10:59
by FStaub
Hi Dylan,

yes, you are correct. That's a bug which I entered in the last update. It will be solved (exactly as you said) in the next update (which is hopefully online today or tomorrow).

Best,
Florian