Add a model here when you are working on it. Mark it as converted when it is tested and pushed.

The err columns are the result of running multi_compare.sh, which generates random parameters for comparison against sasview. double should be at most 5e-14 relative error, single should be 5e-5.

./multi_compare.sh MODEL 200 1d200 mono double
./multi_compare.sh MODEL 200 1d200 mono single

Some models still have issues when run using single precision (lamellarPC for example).

We are not dealing with magnetism at this point as it requires some work on the wrapper functionality. Thus the *Magnetic* column below is for reference as to which models need attention. The *SESANS* column should be used to indicate if the model is a priority (1 highest, 4 lowest).

SasView name	sasmodels name	P(Q)	S(Q)	Multiply	Multiplicity	Magnetic	SESANS	Assigned	Converted	double err	single err
HayterMSAStructure	HayterMSAsq	N	Y	N	N	N	1	RKH	Y	0, 3/200 NaN	FAIL
BarBellModel	barbell	Y	N	N	N	N	3	AJJ	Y	good	good
BCCrystalModel	bcc_paracrystal	Y	N	N	N	N	3	AJJ	Y	good	FAIL
BinaryHSModel	binary_hard_sphere	Y	Y	N	N	N	1	PDB	N
BroadPeakModel	broad_peak	Y	N	N	N	N	4	PDB	Y	good	good
CappedCylinderModel	capped_cylinder	Y	N	N	N	N	3	HP	Y	good	good
CoreShellCylinderModel	core_shell_cylinder	Y	N	Y	N	N	3	HP	Y	good	good
CylinderModel	cylinder	Y	N	Y	N	N - but doc says so	1	HP	Y	good	good
!DABModel	dab	Y	N	N	N	N	1	DM	Y	good	good
EllipsoidModel	ellipsoid	Y	N	Y	N	N	1	HP	Y	good*	good
FCCrystalModel	fcc_paracrystal	Y	N	N	N	N	3	AJJ	Y	good	FAIL
PeakGaussModel	gaussian_peak						4	AJJ	N	good	1e-3
!Guinier	guinier	Y	N	N	N	N	1	JRK	Y	good	good
HardsphereStructure	hardsphere	N	Y	N	N	N	1	RKH	Y needs utest	good	FAIL
HollowCylinderModel	hollow_cylinder	Y	N	N	N	N	2	JRK	Y	good	2e-3
LamellarModel	lamellar	Y	N	N	N	N	2	PAK	Y	good	good
LamellarPSModel	lamellarCaille	Y	N	N	N	N	3	RKH	Y needs utest	good	5e-2
LamellarPSHGModel	lamellarCailleHG	Y	N	N	N	N	3	RKH	Y needs utest	3e-12	5e-3
LamellarFFHGModel	lamellarFFHG	Y	N	N	N	N	3	RKH	Y needs utest	good	good
LamellarPCrystallModel	lamellarPC	Y	N	N	N	N	3	RKH	Y needs utest	good	FAIL
!Lorentz	lorentz	Y	N	N	N	N	1	JRK	Y	good	good
ParallelepipedModel	parallelepiped	Y	N	Y	N	N - but doc says so	2	MAG	Y	good	good
!AbsolutePower_Law	power_law						3	SMK	N	good	good
SphereModel	sphere	Y	N	Y	N	Y	1	HP	Y	good*	good
StickyHSStructure	stickyhardsphere	N	Y	N	N	N	1	RKH	Y needs utest	good	FAIL
TeubnerStrey	teubner_strey						1	PAK	Y	good*	good
TriaxialEllipsoidModel	triaxial_ellipsoid	Y	N	Y	N	Y	3	HP	Y	good*	good
FuzzySphereModel							2	TRN	N
RaspBerryModel							2	AJJ
CoreShellModel						doc says so	2	MD	N
CoreMultiShellModel						doc says so	3
!Core2ndMomentModel							3	SMK	N
MultiShellModel							3	PAK
OnionExpShellModel							3
VesicleModel							2	PDB	N
!SphericalSLDModel							3
LinearPearlsModel							3
PearlNecklaceModel							3	JRK	N
EllipticalCylinderModel							3	MD	N
FlexibleCylinderModel							2		N
!FlexCylEllipXModel							3
CoreShellBicelleModel							2
StackedDisksModel							2	TRN	N
PringleModel							3	AJJ
CoreShellEllipsoidModel							3
!CoreShellEllipsoidXTModel							3
!SCCrystalModel							3
!CSParallelepipedModel							3	MAG	N
RectangularPrismModel							2	MAG	N
RectangularHollowPrismModel							3	MAG	N
RectangularHollowPrismInfThinWallsModel							3	MAG	N
!Debye							1	SMK	N
CorrLength							3
FractalModel							1	PDB
MassFractalModel							1
SurfaceFractalModel							1
MassSurfaceFractal							1
FractalCoreShell							3
GaussLorentzGel							2
!BEPolyelectrolyte							3
GuinierPorod							2	MD	N
PorodModel							3	MD	N
PeakLorentzModel							4	AJJ	N
!Poly_GaussCoil							2	SMK	N
PolyExclVolume							3
!RPA10Model							3	PAK
TwoLorentzian							4
TwoPowerLaw							4
UnifiedPowerRg							4
LineModel							4
GelFitModel							1
StarPolymer							2
ReflectivityModel							4
!ReflectivityIIModel							4
SquareWellStructure							1	PDB
!testmodel							3
!testmodel_2							3
!sum_p1_p2							3
!sum_Ap1_1_Ap2							3
!polynomial5							3
!sph_bessel_jn							3

^* using taylor expansion near 0 for fq = 3 (sin(qr)/qr - cos(qr))/qr²; using quad precision as a cross check, this is more accurate than computing the direct expression for qr < 0.1 in double precision, and so the values will differ from those in sasview. It is much more accurate in single precision. This affects many models, including sphere and ellipsoid models. The cylinder models are using J1(qr)/qr, which has its own numerical expansion and does not suffer from this problem. Should probably settle on one or the other and move it into the library. So even though the double precision values for sphere, ellipsoid and triaxial ellipsoid differ by more than 5e-14, they are still marked good. Similarly, the Teubner-Strey model checks out against the quad precision version and so is marked good.

1GList of standard parameters (name = meaning):

• radius = radius
• sld = scattering length density of particle
• solvent_sld = scattering length density of matrix
• cor_length = correlation length
• exp = exponent (example: porod_exp)
• peak_pos = q_peak or q0 etc
• theta, phi, psi