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
./multi_compare.sh MODEL 200 2d32 mono double
./multi_compare.sh MODEL 200 2d32 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). The Multiply column indicates whether the model can be multiplied by a structure factor. This is currently hardcoded into src/sas/perspectives/fitting/models.py

SasView name	sasmodels name	P(Q)	S(Q)	Multiply	Multiplicity	Magnetic	SESANS	Assigned	double err	single err	double 2D	single 2D
HayterMSAStructure	HayterMSAsq	.	Y	.	.	.	1	RKH	NaNs	FAIL	-	-
BarBellModel	barbell	Y	.	.	.	.	3	AJJ	good	good	FAIL	FAIL
BCCrystalModel	bcc_paracrystal	Y	.	.	.	.	3	AJJ	good	FAIL	9e-14	2e-2
BroadPeakModel	broad_peak	Y	.	.	.	.	4	PDB	good	good	-	-
CappedCylinderModel	capped_cylinder	Y	.	.	.	.	3	HP	good	good	2e-11	3e-3
CoreShellCylinderModel	core_shell_cylinder	Y	.	Y	.	.	3	HP	good	good	9e-14	4e-4
CylinderModel	cylinder	Y	.	Y	.	Y	1	HP	good	good	2e-13	1e-4
DABModel	dab	Y	.	.	.	.	1	DM	good	good	-	-
EllipsoidModel	ellipsoid	Y	.	Y	.	.	1	HP	good*	good	3e-14	1e-4
FCCrystalModel	fcc_paracrystal	Y	.	.	.	.	3	AJJ	good	FAIL	FAIL	FAIL
PeakGaussModel	gaussian_peak	Y	.	.	.	.	4	AJJ	good	1e-3	-	-
Guinier	guinier	Y	.	.	.	.	1	JRK	good	good	-	-
HardsphereStructure	hardsphere	.	Y	.	.	.	1	RKH	good	FAIL	-	-
HollowCylinderModel	hollow_cylinder	Y	.	Y	.	.	2	JRK	good	2e-3	4e-13	2e-3
LamellarModel	lamellar	Y	.	.	.	.	2	PAK	good	good	-	-
LamellarPSModel	lamellarCaille	Y	.	.	.	.	3	RKH	good	5e-2	-	-
LamellarPSHGModel	lamellarCailleHG	Y	.	.	.	.	3	RKH	3e-12	5e-3	-	-
LamellarFFHGModel	lamellarFFHG	Y	.	.	.	.	3	RKH	good	good	-	-
LamellarPCrystallModel	lamellarPC	Y	.	.	.	.	3	RKH	good	FAIL	-	-
Lorentz	lorentz	Y	.	.	.	.	1	JRK	good	good	-	-
ParallelepipedModel	parallelepiped	Y	.	Y	.	Y	2	MAG	good	good	2e-13	8e-5
PearlNecklaceModel	pearl_necklace	Y	.	.	.	.	3	JRK	fail	fail	-	-
!AbsolutePower_Law	power_law	Y	.	.	.	.	3	SMK	good	good	-	-
SphereModel	sphere	Y	.	Y	.	Y	1	HP	good*	good	-	-
StickyHSStructure	stickyhardsphere	.	Y	.	.	.	1	RKH	good	FAIL	-	-
TeubnerStrey	teubner_strey	Y	.	.	.	.	1	PAK	good*	good	-	-
TriaxialEllipsoidModel	triaxial_ellipsoid	Y	.	Y	.	.	3	HP	good*	good	good	good
TwoLorentzian	two_lorentzian	Y	.	.	.	.	4	PR	good	good	-	-
BEPolyelectrolyte		Y	.	.	.	.	3
BinaryHSModel	binary_hard_sphere	Y	.	.	.	.	1	PDB
Core2ndMomentModel		Y	.	.	.	.	3	SMK
CoreMultiShellModel		Y	.	Y	Y	Y	3
CoreShellModel		Y	.	Y	.	Y	2	MD
CoreShellEllipsoidModel		Y	.	Y	.	.	3
CoreShellEllipsoidXTModel		Y	.	Y	.	.	3
CoreShellBicelleModel		Y	.	Y	.	.	2
CorrLength		Y	.	.	.	.	3
CSParallelepipedModel		Y	.	Y	.	.	3	MAG
Debye		Y	.	.	.	.	1	SMK
EllipticalCylinderModel		Y	.	Y	.	.	3	MD
FlexibleCylinderModel		Y	.	.	.	.	2	PR
FlexCylEllipXModel		Y	.	.	.	.	3
FractalModel		Y	.	.	.	.	1	PDB
FractalCoreShell		Y	.	.	.	.	3
FuzzySphereModel		Y	.	.	.	.	2	TRN
GaussLorentzGel		Y	.	.	.	.	2
GelFitModel		Y	.	.	.	.	1
GuinierPorod		Y	.	.	.	.	2	MD
LineModel		Y	.	.	.	.	4
LinearPearlsModel		Y	.	.	.	.	3
MassFractalModel		Y	.	.	.	.	1
MassSurfaceFractal		Y	.	.	.	.	1
MicelleSphCore		Y	.	Y	.	.
MultiShellModel		Y	.	Y	Y	.	3	PAK
MultiplicationModel		.	.	.	.	.		PAK
OnionExpShellModel		Y	.	Y	Y	.	3
PeakLorentzModel		Y	.	.	.	.	4	AJJ
Poly_GaussCoil		Y	.	.	.	.	2	SMK
PolyExclVolume		Y	.	.	.	.	3	PR
PorodModel		Y	.	.	.	.	3	MD
PringleModel		Y	.	.	.	.	3	AJJ
RaspBerryModel		Y	.	.	.	.	2	AJJ
RectangularPrismModel		Y	.	Y	.	.	2	MAG
RectangularHollowPrismModel		Y	.	Y	.	.	3	MAG
RectangularHollowPrismInfThinWallsModel		Y	.	.	.	.	3	MAG
ReflectivityModel		Y	.	.	Y	.	4
ReflectivityIIModel		Y	.	.	Y	.	4
RPA10Model		Y	.	.	Y	.	3	PAK
ReflectivityModel		Y	.	.	.	.	4
SCCrystalModel		Y	.	.	.	.	3
SphericalSLDModel		Y	.	Y	Y	.	3
SquareWellStructure		.	Y	.	.	.	1	PDB
StackedDisksModel		Y	.	Y	.	.	2	TRN
StarPolymer		Y	.	.	.	.	2
SurfaceFractalModel		Y	.	.	.	.	1
TwoPowerLaw		Y	.	.	.	.	4	WP
UnifiedPowerRg		Y	.	.	Y	.	4
VesicleModel		Y	.	Y	.	.	2	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