Changeset a7f909a in sasmodels

Timestamp:

Jan 7, 2016 11:47:25 AM (8 years ago)

Author:

Paul Kienzle <pkienzle@…>

Branches:

master, core_shell_microgels, costrafo411, magnetic_model, release_v0.94, release_v0.95, ticket-1257-vesicle-product, ticket_1156, ticket_1265_superball, ticket_822_more_unit_tests

Children:

82c299f, 55a8316

Parents:

3f832f9

Message:

fix 2d multi-compare

File:

• sasmodels/compare_many.py (modified) (4 diffs)

sasmodels/compare_many.py

@@ -1 +1 @@
 #!/usr/bin/env python
+from __future__ import print_function
+
 import sys
 import traceback
@@ -6 +8 @@
 
 from . import core
-from .kernelcl import environment
-from .compare import (MODELS, randomize_pars, suppress_pd, eval_sasview,
-                      eval_opencl, eval_ctypes, make_data, get_demo_pars,
-                      columnize, constrain_pars, constrain_new_to_old,
-                      make_engine)
+from . import generate
+from .compare import (MODELS, randomize_pars, suppress_pd, make_data,
+                      make_engine, get_demo_pars, columnize,
+                      constrain_pars, constrain_new_to_old)
 
 def calc_stats(target, value, index):
@@ -47 +48 @@
 def compare_instance(name, data, index, N=1, mono=True, cutoff=1e-5,
                      base='sasview', comp='double'):
+    is2D = hasattr(data, 'qx_data')
     model_definition = core.load_model_definition(name)
     pars = get_demo_pars(model_definition)
-    header = '\n"Model","%s","Count","%d"'%(name, N)
+    header = ('\n"Model","%s","Count","%d","Dimension","%s"'
+              % (name, N, "2D" if is2D else "1D"))
     if not mono: header += ',"Cutoff",%g'%(cutoff,)
     print(header)
+
+    if is2D:
+        info = generate.make_info(model_definition)
+        partype = info['partype']
+        if not partype['orientation'] and not partype['magnetic']:
+            print(',"1-D only"')
+            return
 
     # Some not very clean macros for evaluating the models and checking the
@@ -88 +98 @@
     max_diff = [0]
     for k in range(N):
-        print("%s %d"%(name, k))
+        print("%s %d"%(name, k), file=sys.stderr)
         seed = np.random.randint(1e6)
         pars_i = randomize_pars(pars, seed)