[5428233] | 1 | # -*- coding: utf-8 -*- |
---|
| 2 | """ |
---|
| 3 | Created on Tue Feb 17 11:43:56 2015 |
---|
| 4 | |
---|
| 5 | @author: David |
---|
| 6 | """ |
---|
| 7 | |
---|
| 8 | import unittest |
---|
| 9 | import warnings |
---|
| 10 | import numpy as np |
---|
| 11 | |
---|
| 12 | from os.path import basename, dirname, join as joinpath |
---|
| 13 | from glob import glob |
---|
| 14 | |
---|
| 15 | try: |
---|
| 16 | from .kernelcl import load_model |
---|
| 17 | except ImportError,exc: |
---|
| 18 | warnings.warn(str(exc)) |
---|
| 19 | warnings.warn("using ctypes instead") |
---|
| 20 | from .kerneldll import load_model |
---|
| 21 | |
---|
| 22 | def load_kernel(model, dtype='single'): |
---|
| 23 | kernel = load_model(model, dtype=dtype) |
---|
| 24 | kernel.info['defaults'] = dict((p[0],p[2]) for p in kernel.info['parameters']) |
---|
| 25 | return kernel |
---|
| 26 | |
---|
| 27 | def get_weights(model, pars, name): |
---|
| 28 | from . import weights |
---|
| 29 | |
---|
| 30 | relative = name in model.info['partype']['pd-rel'] |
---|
| 31 | disperser = pars.get(name+"_pd_type", "gaussian") |
---|
| 32 | value = pars.get(name, model.info['defaults'][name]) |
---|
| 33 | width = pars.get(name+"_pd", 0.0) |
---|
| 34 | npts = pars.get(name+"_pd_n", 30) |
---|
| 35 | nsigma = pars.get(name+"_pd_nsigma", 3.0) |
---|
| 36 | v,w = weights.get_weights( |
---|
| 37 | disperser, npts, width, nsigma, |
---|
| 38 | value, model.info['limits'][name], relative) |
---|
| 39 | return v,w/np.sum(w) |
---|
| 40 | |
---|
| 41 | def eval_kernel(kernel, q, pars, cutoff=1e-5): |
---|
| 42 | input = kernel.make_input(q) |
---|
| 43 | finput = kernel(input) |
---|
| 44 | |
---|
| 45 | fixed_pars = [pars.get(name, finput.info['defaults'][name]) |
---|
| 46 | for name in finput.fixed_pars] |
---|
| 47 | pd_pars = [get_weights(finput, pars, p) for p in finput.pd_pars] |
---|
| 48 | return finput(fixed_pars, pd_pars, cutoff) |
---|
| 49 | |
---|
| 50 | def annotate_exception(exc, msg): |
---|
| 51 | """ |
---|
| 52 | Add an annotation to the current exception, which can then be forwarded |
---|
| 53 | to the caller using a bare "raise" statement to reraise the annotated |
---|
| 54 | exception. |
---|
| 55 | Example:: |
---|
| 56 | >>> D = {} |
---|
| 57 | >>> try: |
---|
| 58 | ... print D['hello'] |
---|
| 59 | ... except Exception,exc: |
---|
| 60 | ... annotate_exception(exc, "while accessing 'D'") |
---|
| 61 | ... raise |
---|
| 62 | Traceback (most recent call last): |
---|
| 63 | ... |
---|
| 64 | KeyError: "hello while accessing 'D'" |
---|
| 65 | """ |
---|
| 66 | args = exc.args |
---|
| 67 | if not args: |
---|
[5ca9762] | 68 | exc.args = (msg,) |
---|
[5428233] | 69 | else: |
---|
[5ca9762] | 70 | try: |
---|
| 71 | arg0 = " ".join((args[0],msg)) |
---|
| 72 | exc.args = tuple([arg0] + list(args[1:])) |
---|
| 73 | except: |
---|
| 74 | exc.args = (" ".join((str(exc),msg)),) |
---|
[5428233] | 75 | |
---|
| 76 | def suite(): |
---|
| 77 | root = dirname(__file__) |
---|
| 78 | files = sorted(glob(joinpath(root, 'models', "[a-zA-Z]*.py"))) |
---|
| 79 | models_names = [basename(f)[:-3] for f in files] |
---|
| 80 | |
---|
| 81 | suite = unittest.TestSuite() |
---|
| 82 | |
---|
| 83 | for model_name in models_names: |
---|
| 84 | module = __import__('sasmodels.models.' + model_name) |
---|
| 85 | module = getattr(module, 'models', None) |
---|
| 86 | |
---|
| 87 | model = getattr(module, model_name, None) |
---|
[5ca9762] | 88 | smoke_tests = [[{},0.1,None],[{},(0.1,0.1),None]] |
---|
| 89 | tests = smoke_tests + getattr(model, 'tests', []) |
---|
[5428233] | 90 | |
---|
| 91 | if tests: |
---|
| 92 | #print '------' |
---|
| 93 | #print 'found tests in', model_name |
---|
| 94 | #print '------' |
---|
| 95 | |
---|
| 96 | kernel = load_kernel(model) |
---|
| 97 | suite.addTest(ModelTestCase(model_name, kernel, tests)) |
---|
| 98 | |
---|
| 99 | return suite |
---|
| 100 | |
---|
| 101 | class ModelTestCase(unittest.TestCase): |
---|
| 102 | |
---|
| 103 | def __init__(self, model_name, kernel, tests): |
---|
| 104 | unittest.TestCase.__init__(self) |
---|
| 105 | |
---|
| 106 | self.model_name = model_name |
---|
| 107 | self.kernel = kernel |
---|
| 108 | self.tests = tests |
---|
| 109 | |
---|
| 110 | def runTest(self): |
---|
| 111 | #print '------' |
---|
| 112 | #print self.model_name |
---|
| 113 | #print '------' |
---|
| 114 | try: |
---|
| 115 | for test in self.tests: |
---|
| 116 | params = test[0] |
---|
| 117 | Q = test[1] |
---|
| 118 | I = test[2] |
---|
| 119 | |
---|
| 120 | if not isinstance(Q, list): |
---|
| 121 | Q = [Q] |
---|
| 122 | if not isinstance(I, list): |
---|
| 123 | I = [I] |
---|
| 124 | |
---|
| 125 | if isinstance(Q[0], tuple): |
---|
[5ca9762] | 126 | Qx,Qy = zip(*Q) |
---|
| 127 | Q_vectors = [np.array(Qx), np.array(Qy)] |
---|
[5428233] | 128 | else: |
---|
[5ca9762] | 129 | Q_vectors = [np.array(Q)] |
---|
[5428233] | 130 | |
---|
[5ca9762] | 131 | self.assertEqual(len(I), len(Q)) |
---|
[5428233] | 132 | |
---|
[5ca9762] | 133 | Iq = eval_kernel(self.kernel, Q_vectors, params) |
---|
[5428233] | 134 | |
---|
| 135 | self.assertGreater(len(Iq), 0) |
---|
| 136 | self.assertEqual(len(I), len(Iq)) |
---|
| 137 | |
---|
| 138 | for q, i, iq in zip(Q, I, Iq): |
---|
[5ca9762] | 139 | if i is None: continue # smoke test --- make sure it runs |
---|
[5428233] | 140 | err = np.abs(i - iq) |
---|
| 141 | nrm = np.abs(i) |
---|
| 142 | |
---|
| 143 | self.assertLess(err * 10**5, nrm, 'q:%s; expected:%s; actual:%s' % (q, i, iq)) |
---|
| 144 | |
---|
| 145 | except Exception,exc: |
---|
| 146 | annotate_exception(exc, '\r\nModel: %s' % self.model_name) |
---|
| 147 | raise |
---|
| 148 | |
---|
| 149 | def main(): |
---|
| 150 | #unittest.main() |
---|
| 151 | runner = unittest.TextTestRunner() |
---|
| 152 | runner.run(suite()) |
---|
| 153 | |
---|
| 154 | if __name__ == "__main__": |
---|
| 155 | main() |
---|