Testing `galactic_expanding_shell` branch

[sona-patel]

Hi @andreatramacere,

Very glad to see the JetSet evolving so fast and now also having preliminary implementation for Galactic objects! 👏🏻

Isn't that the very first line on the documentation page should read as 1) GalacticBeamed to hadle galatcit objects that are not beamed 2) GalacticUnbeamed to hadle galatcit objects that are not beamed ?

Thanks!

[andreatramacere]

Hi @sona-patel, yes you are right, in any case, I did some more improvements not yet released, such as adding an expanding shell geometry. If you are interested you can try the branch 'galactic_expanding_shell', but it is still experimental, even though a postdoc in Japan used it to reproduce the RS Ophiuchi results published in the MAGIC paper. I guess I will have something ready to experiment early next year, but if you are interested I can provide you with some notebooks after the Xmas break.

[sona-patel]

That's nice! I will also give more serious try early next year! - Thanks!

[sona-patel]

Hi @andreatramacere , I tested the 'galactic_expanding_shell' branch of Jetset! I can run smoothly galactic object models. However, I get the following errors while running the test:

Failures in '$ pytest --pyargs -vvv jetset.tests.test_users::TestUser':

==================================== test session starts ==================================== platform darwin -- Python 3.9.16, pytest-7.2.1, pluggy-1.0.0 -- /Users/sonal/anaconda3/envs/jetset_expanding_shell/bin/python3.9 cachedir: .pytest_cache rootdir: /Users/sonal/test_jetset plugins: anyio-3.5.0 collected 5 items ::TestUser::test_jet <- ../anaconda3/envs/jetset_expanding_shell/lib/python3.9/site-packages/jetset/tests/test_users.py -------------------------------------------------------------------------------- model description: -------------------------------------------------------------------------------- type: Jet FAILED [ 20%] ::TestUser::test_emitters <- ../anaconda3/envs/jetset_expanding_shell/lib/python3.9/site-packages/jetset/tests/test_users.py PASSED [ 40%] ::TestUser::test_dep_pars <- ../anaconda3/envs/jetset_expanding_shell/lib/python3.9/site-packages/jetset/tests/test_users.py FAILED [ 60%] ::TestUser::test_composit_model <- ../anaconda3/envs/jetset_expanding_shell/lib/python3.9/site-packages/jetset/tests/test_users.py PASSED [ 80%] ::TestUser::test_composit_ebl <- ../anaconda3/envs/jetset_expanding_shell/lib/python3.9/site-packages/jetset/tests/test_users.py PASSED [100%] ========================================= FAILURES ========================================== _____________________________________ TestUser.test_jet _____________________________________ self = plot = def test_jet(self,plot=plot): t=TestJets() t.test_jet(plot=plot) t.test_set_N_from_nuFnu(plot=plot) t.test_EC(plot=plot) > t.test_hadronic_jet(plot=plot) ../anaconda3/envs/jetset_expanding_shell/lib/python3.9/site-packages/jetset/tests/test_users.py:29: _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ ../anaconda3/envs/jetset_expanding_shell/lib/python3.9/site-packages/jetset/tests/test_jet_model.py:18: in test_hadronic_jet j = Jet(proton_distribution='plc') _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ self = <[AttributeError("'Jet' object has no attribute 'name'") raised in repr()] Jet object at 0x7fda136abaf0> cosmo = None, name = None, emitters_type = 'electrons', emitters_distribution = 'plc' emitters_distribution_log_values = False, beaming_expr = 'delta', T_esc_e_second = None jet_workplace = None, verbose = None, clean_work_dir = True, electron_distribution = None proton_distribution = 'plc', electron_distribution_log_values = None proton_distribution_log_values = None, geometry = 'spherical' def __init__(self, cosmo=None, name=None, emitters_type='electrons', emitters_distribution='plc', emitters_distribution_log_values=False, beaming_expr='delta', T_esc_e_second=None, jet_workplace=None, verbose=None, clean_work_dir=True, electron_distribution=None, proton_distribution=None, electron_distribution_log_values=None, proton_distribution_log_values=None, geometry='spherical'): """ Parameters ---------- cosmo name emitters_type emitters_distribution emitters_distribution_log_values beaming_expr T_esc_e_second jet_workplace verbose clean_work_dir electron_distribution proton_distribution electron_distribution_log_values proton_distribution_log_values """ if electron_distribution is not None and emitters_type == 'protons': raise RuntimeError("if you provide an electron_distribution, you can't set emitters_type=",emitters_type, 'use emitters_distribution instead') if proton_distribution is not None and emitters_type == 'electrons': > raise RuntimeError("if you provide a proton_distribution, you can't set emitters_type=",emitters_type, 'use emitters_distribution instead') E RuntimeError: ("if you provide a proton_distribution, you can't set emitters_type=", 'electrons', 'use emitters_distribution instead') ../anaconda3/envs/jetset_expanding_shell/lib/python3.9/site-packages/jetset/jet_model.py:1625: RuntimeError ----------------------------------- Captured stdout call ------------------------------------ --------> test_jet name type units val ----------------- ----------------------------- --------- ------------- BulkLorentzFactor jet-bulk-factor 1.000000e+01 U_e Energy dens. blob rest. frame erg / cm3 1.299489e-03 U_p_cold Energy dens. blob rest. frame erg / cm3 1.503276e-02 U_B Energy dens. blob rest. frame erg / cm3 3.978874e-04 U_Synch Energy dens. blob rest. frame erg / cm3 1.187517e-06 U_Synch_DRF Energy dens. disk rest. frame erg / cm3 1.187517e-02 U_Disk Energy dens. blob rest. frame erg / cm3 0.000000e+00 U_BLR Energy dens. blob rest. frame erg / cm3 0.000000e+00 U_DT Energy dens. blob rest. frame erg / cm3 0.000000e+00 U_CMB Energy dens. blob rest. frame erg / cm3 0.000000e+00 U_Star Energy dens. blob rest. frame erg / cm3 0.000000e+00 U_Disk_DRF Energy dens. disk rest. frame erg / cm3 0.000000e+00 U_BLR_DRF Energy dens. disk rest. frame erg / cm3 0.000000e+00 U_DT_DRF Energy dens. disk rest. frame erg / cm3 0.000000e+00 U_CMB_DRF Energy dens. disk rest. frame erg / cm3 0.000000e+00 U_Star_DRF Energy dens. disk rest. frame erg / cm3 0.000000e+00 U_seed_tot Energy dens. blob rest. frame erg / cm3 1.187517e-06 L_Sync_rf Lum. blob rest. frame. erg / s 1.118434e+37 L_SSC_rf Lum. blob rest. frame. erg / s 7.050994e+34 L_EC_Disk_rf Lum. blob rest. frame. erg / s 0.000000e+00 L_EC_BLR_rf Lum. blob rest. frame. erg / s 0.000000e+00 L_EC_DT_rf Lum. blob rest. frame. erg / s 0.000000e+00 L_EC_CMB_rf Lum. blob rest. frame. erg / s 0.000000e+00 L_EC_Star_rf Lum. blob rest. frame. erg / s 6.953119e-310 jet_L_Sync jet Lum. erg / s 2.782070e+38 jet_L_SSC jet Lum. erg / s 1.753913e+36 jet_L_EC_Disk jet Lum. erg / s 0.000000e+00 jet_L_EC_BLR jet Lum. erg / s 0.000000e+00 jet_L_EC_Star jet Lum. erg / s 6.945305e-310 jet_L_EC_DT jet Lum. erg / s 0.000000e+00 jet_L_EC_CMB jet Lum. erg / s 0.000000e+00 jet_L_pp_gamma jet Lum. erg / s 0.000000e+00 jet_L_rad jet Lum. erg / s 2.799609e+38 jet_L_kin jet Lum. erg / s 3.826258e+42 jet_L_tot jet Lum. erg / s 3.919754e+42 jet_L_e jet Lum. erg / s 3.044393e+41 jet_L_B jet Lum. erg / s 9.321554e+40 jet_L_p_cold jet Lum. erg / s 3.521819e+42 NH_cold_to_rel_e cold_p_to_rel_e_ratio 1.000000e-01 model name name par type units val phys. bound. min phys. bound. max log frozen ------------ ---------------- --------------------- --------------- ------------ ---------------- ---------------- ----- ------ jet_leptonic gmin low-energy-cut-off lorentz-factor* 2.000000e+00 1.000000e+00 1.000000e+09 False False jet_leptonic gmax high-energy-cut-off lorentz-factor* 1.000000e+06 1.000000e+00 1.000000e+15 False False jet_leptonic N emitters_density 1 / cm3 1.000000e+02 0.000000e+00 -- False False jet_leptonic gamma_cut turn-over-energy lorentz-factor* 1.000000e+04 1.000000e+00 1.000000e+09 False False jet_leptonic p LE_spectral_slope 2.000000e+00 -1.000000e+01 1.000000e+01 False False jet_leptonic R region_size cm 5.000000e+15 1.000000e+03 1.000000e+30 False False jet_leptonic R_H region_position cm 1.000000e+17 0.000000e+00 -- False True jet_leptonic B magnetic_field gauss 1.000000e-01 0.000000e+00 -- False False jet_leptonic NH_cold_to_rel_e cold_p_to_rel_e_ratio 1.000000e-01 0.000000e+00 -- False True jet_leptonic beam_obj beaming lorentz-factor* 1.000000e+01 1.000000e-04 -- False False jet_leptonic z_cosm redshift 1.000000e-01 0.000000e+00 -- False False name type units val ----------------- ----------------------------- --------- ------------- BulkLorentzFactor jet-bulk-factor 1.000000e+01 U_e Energy dens. blob rest. frame erg / cm3 1.299489e-03 U_p_cold Energy dens. blob rest. frame erg / cm3 1.503276e-02 U_B Energy dens. blob rest. frame erg / cm3 3.978874e-04 U_Synch Energy dens. blob rest. frame erg / cm3 1.187517e-06 U_Synch_DRF Energy dens. disk rest. frame erg / cm3 1.187517e-02 U_Disk Energy dens. blob rest. frame erg / cm3 0.000000e+00 U_BLR Energy dens. blob rest. frame erg / cm3 0.000000e+00 U_DT Energy dens. blob rest. frame erg / cm3 0.000000e+00 U_CMB Energy dens. blob rest. frame erg / cm3 0.000000e+00 U_Star Energy dens. blob rest. frame erg / cm3 0.000000e+00 U_Disk_DRF Energy dens. disk rest. frame erg / cm3 0.000000e+00 U_BLR_DRF Energy dens. disk rest. frame erg / cm3 0.000000e+00 U_DT_DRF Energy dens. disk rest. frame erg / cm3 0.000000e+00 U_CMB_DRF Energy dens. disk rest. frame erg / cm3 0.000000e+00 U_Star_DRF Energy dens. disk rest. frame erg / cm3 0.000000e+00 U_seed_tot Energy dens. blob rest. frame erg / cm3 1.187517e-06 L_Sync_rf Lum. blob rest. frame. erg / s 1.118434e+37 L_SSC_rf Lum. blob rest. frame. erg / s 7.050994e+34 L_EC_Disk_rf Lum. blob rest. frame. erg / s 0.000000e+00 L_EC_BLR_rf Lum. blob rest. frame. erg / s 0.000000e+00 L_EC_DT_rf Lum. blob rest. frame. erg / s 0.000000e+00 L_EC_CMB_rf Lum. blob rest. frame. erg / s 0.000000e+00 L_EC_Star_rf Lum. blob rest. frame. erg / s 6.953119e-310 jet_L_Sync jet Lum. erg / s 2.782070e+38 jet_L_SSC jet Lum. erg / s 1.753913e+36 jet_L_EC_Disk jet Lum. erg / s 0.000000e+00 jet_L_EC_BLR jet Lum. erg / s 0.000000e+00 jet_L_EC_Star jet Lum. erg / s 6.945305e-310 jet_L_EC_DT jet Lum. erg / s 0.000000e+00 jet_L_EC_CMB jet Lum. erg / s 0.000000e+00 jet_L_pp_gamma jet Lum. erg / s 0.000000e+00 jet_L_rad jet Lum. erg / s 2.799609e+38 jet_L_kin jet Lum. erg / s 3.826258e+42 jet_L_tot jet Lum. erg / s 3.919754e+42 jet_L_e jet Lum. erg / s 3.044393e+41 jet_L_B jet Lum. erg / s 9.321554e+40 jet_L_p_cold jet Lum. erg / s 3.521819e+42 NH_cold_to_rel_e cold_p_to_rel_e_ratio 1.000000e-01 --------> test_set_N_from_nuFnu --------> test_EC --------> test_hadronic_jet __________________________________ TestUser.test_dep_pars ___________________________________ args = (,), kwargs = {'val': 2.0} message = 'the jetkernel failed\n\n exception message: \ns1+1 argument after ** must be a mapping, not list' def func_wrapper(*args, **kwargs): try: > return func(*args, **kwargs) ../anaconda3/envs/jetset_expanding_shell/lib/python3.9/site-packages/jetset/utils.py:93: _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ self = , keywords = {'val': 2.0} @safe_run def set(self,**keywords): """ overrides the :meth:`.ModelParameter.set` method in order to propagate the parameter value to the BlazarSED object instance """ > super(JetParameter,self).set(**keywords ) ../anaconda3/envs/jetset_expanding_shell/lib/python3.9/site-packages/jetset/jet_paramters.py:110: _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ self = skip_dep_par_warning = False, args = (), keywords = {'val': 2.0}, keys = dict_keys(['val']) kw = 'val', p = def set(self, *args, skip_dep_par_warning=False, **keywords): """ sets a parameter value checking for physical boundaries Parameters: keywords of the constructor """ keys = keywords.keys() if self.immutable is False or skip_dep_par_warning is True: pass else: #warnings.warn('\n\n *** you are trying to set a dependent parameter:%s *** \n'%self.name) raise RuntimeError('\n\n *** you are trying to set a dependent parameter:%s *** \n'%self.name) #return for kw in keys: if kw in self.allowed_keywords.keys() : if kw == 'val': if self.allowed_values is not None: if keywords[kw] not in self.allowed_values: raise RuntimeError('parameter %s' %(self.name), 'the value', keywords[kw] , 'is not in the allowed list',self.allowed_values) self._val.val = keywords[kw] if self._depending_pars is not []: for p in self._depending_pars: #print("==> setting dep par",p.name, 'to',p._func(),'p=',p,'master',self,'name',self.name) > p.set(val=p._func(),skip_dep_par_warning=True) ../anaconda3/envs/jetset_expanding_shell/lib/python3.9/site-packages/jetset/model_parameters.py:429: _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ self = def _eval_par_func(self): #transform par name and value into a local var #print('working on ',self.name) #TODO:THIS HOLDS ONLY FOR NUMPY <1.22, should be removed warnings.filterwarnings('ignore', message='invalid value encountered in reciprocal*') if type(self._depending_par_expr) == str: _par_values= [None]*len(self._master_pars) for ID, _user_par_ in enumerate(self._master_pars): _par_values[ID] = _user_par_.val_lin*_user_par_.units #print('==> _eval_par_func',_user_par_.name,_par_values[ID]) exec(_user_par_.name + '=_par_values[ID]') res = eval(self._depending_par_expr) elif callable(self._depending_par_expr) is True: _par_values={} for ID, _user_par_ in enumerate(self._master_pars): _par_values[_user_par_.name] = _user_par_.val_lin*_user_par_.units #print(_user_par_.name,_user_par_.val_lin) res=self._depending_par_expr(**_par_values) #TODO:THIS HOLDS ONLY FOR NUMPY <1.22, should be removed > res=self._depending_par_expr(**_par_values) E TypeError: s1+1 argument after ** must be a mapping, not list ../anaconda3/envs/jetset_expanding_shell/lib/python3.9/site-packages/jetset/model_parameters.py:378: TypeError During handling of the above exception, another exception occurred: self = plot = def test_dep_pars(self,plot=plot): t=TestDependingParameters() > t.integration_suite() ../anaconda3/envs/jetset_expanding_shell/lib/python3.9/site-packages/jetset/tests/test_users.py:37: _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ ../anaconda3/envs/jetset_expanding_shell/lib/python3.9/site-packages/jetset/tests/test_depending_parameters.py:7: in integration_suite self._all(plot=plot) ../anaconda3/envs/jetset_expanding_shell/lib/python3.9/site-packages/jetset/tests/base_class.py:9: in _all m[1](kwargs) ../anaconda3/envs/jetset_expanding_shell/lib/python3.9/site-packages/jetset/tests/test_depending_parameters.py:33: in test_dep_par j.make_dependent_par(par='s2', depends_on=['s1'], par_expr='s1+1') ../anaconda3/envs/jetset_expanding_shell/lib/python3.9/site-packages/jetset/base_model.py:385: in make_dependent_par m.val=m.val ../anaconda3/envs/jetset_expanding_shell/lib/python3.9/site-packages/jetset/model_parameters.py:257: in val self.set(val=val) _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ args = (,), kwargs = {'val': 2.0} message = 'the jetkernel failed\n\n exception message: \ns1+1 argument after ** must be a mapping, not list' def func_wrapper(*args, **kwargs): try: return func(*args, **kwargs) except Exception as e: message = 'the jetkernel failed\n' message += '\n exception message: \n' message += '%s'%str(e) new_version_warning() > raise JetkerneltException(message=message) E jetset.utils.JetkerneltException: the jetkernel failed E E exception message: E s1+1 argument after ** must be a mapping, not list ../anaconda3/envs/jetset_expanding_shell/lib/python3.9/site-packages/jetset/utils.py:100: JetkerneltException ----------------------------------- Captured stdout call ------------------------------------ name par type units val phys. bound. min phys. bound. max log frozen ----------- ------------------- --------------- ------------ ---------------- ---------------- ----- ------ gmin low-energy-cut-off lorentz-factor* 2.000000e+00 1.000000e+00 1.000000e+09 False False gmax high-energy-cut-off lorentz-factor* 1.000000e+06 1.000000e+00 1.000000e+15 False False N emitters_density 1 / cm3 1.000000e+02 0.000000e+00 -- False False gamma_break turn-over-energy lorentz-factor* 1.000000e+03 1.000000e+00 -- False False s1 LE_spectral_slope 2.500000e+00 -1.000000e+01 1.000000e+01 False False s2 LE_spectral_slope 3.200000e+00 -1.000000e+01 1.000000e+01 False False name par type units val phys. bound. min phys. bound. max log frozen ----------- ------------------- --------------- ------------ ---------------- ---------------- ----- ------ gmin low-energy-cut-off lorentz-factor* 2.000000e+00 1.000000e+00 1.000000e+09 False False gmax high-energy-cut-off lorentz-factor* 1.000000e+06 1.000000e+00 1.000000e+15 False False N emitters_density 1 / cm3 1.000000e+02 0.000000e+00 -- False False gamma_break turn-over-energy lorentz-factor* 1.000000e+03 1.000000e+00 -- False False s1 LE_spectral_slope 2.000000e+00 -1.000000e+01 1.000000e+01 False False s2 LE_spectral_slope 3.500000e+00 -1.000000e+01 1.000000e+01 False False ===================================== warnings summary ====================================== : 104 warnings /Users/sonal/anaconda3/envs/jetset_expanding_shell/lib/python3.9/site-packages/jetset/model_parameters.py:1121: DeprecationWarning: `np.float` is a deprecated alias for the builtin `float`. To silence this warning, use `float` by itself. Doing this will not modify any behavior and is safe. If you specifically wanted the numpy scalar type, use `np.float64` here. Deprecated in NumPy 1.20; for more details and guidance: https://numpy.org/devdocs/release/1.20.0-notes.html#deprecations t[n] = MaskedColumn(t[n].data, name=n, dtype=np.float) : 67 warnings /Users/sonal/anaconda3/envs/jetset_expanding_shell/lib/python3.9/site-packages/jetset/model_parameters.py:1119: DeprecationWarning: `np.float` is a deprecated alias for the builtin `float`. To silence this warning, use `float` by itself. Doing this will not modify any behavior and is safe. If you specifically wanted the numpy scalar type, use `np.float64` here. Deprecated in NumPy 1.20; for more details and guidance: https://numpy.org/devdocs/release/1.20.0-notes.html#deprecations t[n] = MaskedColumn(t[n].data, name=n, dtype=np.float, mask=t[n].data==None) ::TestUser::test_jet ::TestUser::test_jet ::TestUser::test_emitters ::TestUser::test_emitters /Users/sonal/anaconda3/envs/jetset_expanding_shell/lib/python3.9/site-packages/jetset/jet_emitters.py:213: DeprecationWarning: `np.bool` is a deprecated alias for the builtin `bool`. To silence this warning, use `bool` by itself. Doing this will not modify any behavior and is safe. If you specifically wanted the numpy scalar type, use `np.bool_` here. Deprecated in NumPy 1.20; for more details and guidance: https://numpy.org/devdocs/release/1.20.0-notes.html#deprecations msk = np.ones(x.shape, dtype=np.bool) ::TestUser::test_dep_pars /Users/sonal/anaconda3/envs/jetset_expanding_shell/lib/python3.9/site-packages/jetset/utils.py:54: UserWarning: ******************************************************************************** Something wrong has happened. Please, look at the exception message. ******************************************************************************** warnings.warn(m) ::TestUser::test_composit_model ::TestUser::test_composit_ebl ::TestUser::test_composit_ebl /Users/sonal/anaconda3/envs/jetset_expanding_shell/lib/python3.9/site-packages/jetset/template_2Dmodel.py:295: DeprecationWarning: `np.float` is a deprecated alias for the builtin `float`. To silence this warning, use `float` by itself. Doing this will not modify any behavior and is safe. If you specifically wanted the numpy scalar type, use `np.float64` here. Deprecated in NumPy 1.20; for more details and guidance: https://numpy.org/devdocs/release/1.20.0-notes.html#deprecations x_values = np.array(data.meta['redshift'], dtype=np.float) ::TestUser::test_composit_model ::TestUser::test_composit_ebl /Users/sonal/anaconda3/envs/jetset_expanding_shell/lib/python3.9/site-packages/jetset/model_manager.py:147: UserWarning: no cosmology defined, using default FlatLambdaCDM(name="Planck13", H0=67.8 km / (Mpc s), Om0=0.307, Tcmb0=2.725 K, Neff=3.05, m_nu=[0. 0. 0.06] eV, Ob0=0.0483) warnings.warn('no cosmology defined, using default %s'%self.cosmo) ::TestUser::test_composit_ebl ::TestUser::test_composit_ebl /Users/sonal/anaconda3/envs/jetset_expanding_shell/lib/python3.9/site-packages/jetset/template_2Dmodel.py:293: DeprecationWarning: `np.float` is a deprecated alias for the builtin `float`. To silence this warning, use `float` by itself. Doing this will not modify any behavior and is safe. If you specifically wanted the numpy scalar type, use `np.float64` here. Deprecated in NumPy 1.20; for more details and guidance: https://numpy.org/devdocs/release/1.20.0-notes.html#deprecations x_values = np.array(data.meta['REDSHIFT'], dtype=np.float) -- Docs: https://docs.pytest.org/en/stable/how-to/capture-warnings.html ================================== short test summary info ================================== FAILED ::TestUser::test_jet - RuntimeError: ("if you provide a proton_distribution, you can't set emitters_type=", 'el... FAILED ::TestUser::test_dep_pars - jetset.utils.JetkerneltException: the jetkernel failed ======================== 2 failed, 3 passed, 183 warnings in 35.01s =========================

Thanks!

[andreatramacere]

@sona-patel, Hi, I did not finalise all the tests. I have pushed a change, now the test should work. You can test also multithreading, and it will be faster, for example, if you have a jet model you can:

jet._blob.N_THREADS=8

This will build 8 threads in the C code. On my laptop, with N_THREADS=8, I get ~ 3 times faster. This is still experimental, the thumb rule is to use N_THREADS~2*(number of your CPU cores)

These days I am super busy, but I should have a release during the next two months