mceq-project / MCEq

Matrix cascade equation core code
BSD 3-Clause "New" or "Revised" License
30 stars 33 forks source link

import problem with the newest version #13

Closed pkalaczynski closed 4 years ago

pkalaczynski commented 4 years ago

Hi, with the newest version (1.1.1.) I get an error. When I execute this code:

from MCEq.core import config, MCEqRun
import crflux.models as crf
import matplotlib.pyplot as plt

# Initalize MCEq by creating the user interface object MCEqRun

mceq = MCEqRun(
    interaction_model='SIBYLL23C',
    primary_model = (crf.GaisserStanevTilav,'3-gen'),
    theta_deg = 0.
)

# Solve the equation system
mceq.solve()

The output is:

MCEqRun::set_interaction_model(): SIBYLL23C
ParticleManager::_init_default_tracking(): Initializing default tracking categories (pi, K, mu)
MCEqRun::set_density_model(): Setting density profile to CORSIKA ('BK_USStd', None)
MCEqRun::set_primary_model(): GaisserStanevTilav 3-gen
---------------------------------------------------------------------------
ImportError                               Traceback (most recent call last)
<ipython-input-4-8652a7e1c444> in <module>()
     13 
     14 # Solve the equation system
---> 15 mceq.solve()

D:\Anaconda\lib\site-packages\MCEq\core.py in solve(self, int_grid, grid_var, **kwargs)
    747                     config.kernel_config))
    748 
--> 749         self._solution, self.grid_sol = kernel(*args)
    750 
    751         info(

D:\Anaconda\lib\site-packages\MCEq\solvers.py in solv_MKL_sparse(nsteps, dX, rho_inv, int_m, dec_m, phi, grid_idcs)
    193 
    194     from ctypes import c_int, c_char, POINTER, byref
--> 195     from mceq_config import mkl
    196 
    197     gemv = None

ImportError: cannot import name 'mkl'
afedynitch commented 4 years ago

Sorry, the tests are neither great nor very detailed and mkl was unfortunately not part of standard tests. I fixed this in version 1.1.2. Please update through pip. Thanks.