Issue loading CBLAS library

Dear all,

My name is Antonio Garcia-Uceda, and I’m a researcher on Flux Reconstruction methods. I’m using pyFR as part of my project.

I’m trying to run pyFR with backend “opemp”, but I’m coming across some issues when using a CBLAS library. I get the following error message on the screen:

Let me know if you know the source of the problem. Thanks a lotin advance.

Best regards,

Antonio

File "/develop/home/antonio/Software_Local/pyFR/pyfr/scripts/pyfr-sim", line 112, in <module>
main()
File "/usr/lib/pymodules/python2.7/mpmath/ctx_mp.py", line 1248, in g
return f(*args, **kwargs)
File "/develop/home/antonio/Software_Local/pyFR/pyfr/scripts/pyfr-sim", line 88, in main
solver = get_solver(backend, rallocs, mesh, soln, cfg)
File "/develop/home/antonio/Software_Local/pyFR/pyfr/solvers/__init__.py", line 14, in get_solver
return get_integrator(backend, systemcls, rallocs, mesh, initsoln, cfg)
File "/develop/home/antonio/Software_Local/pyFR/pyfr/integrators/__init__.py", line 29, in get_integrator
return integrator(backend, systemcls, rallocs, mesh, initsoln, cfg)
File "/develop/home/antonio/Software_Local/pyFR/pyfr/integrators/controllers.py", line 11, in __init__
super(BaseController, self).__init__(*args, **kwargs)
File "/develop/home/antonio/Software_Local/pyFR/pyfr/integrators/steppers.py", line 11, in __init__
super(BaseStepper, self).__init__(*args, **kwargs)
File "/develop/home/antonio/Software_Local/pyFR/pyfr/integrators/writers.py", line 68, in __init__
super(FileWriter, self).__init__(*args, **kwargs)
File "/develop/home/antonio/Software_Local/pyFR/pyfr/integrators/writers.py", line 16, in __init__
super(BaseWriter, self).__init__(*args, **kwargs)
File "/develop/home/antonio/Software_Local/pyFR/pyfr/integrators/base.py", line 50, in __init__
self._system = systemcls(backend, rallocs, mesh, initsoln, nreg, cfg)
File "/develop/home/antonio/Software_Local/pyFR/pyfr/solvers/base/system.py", line 31, in __init__
eles, elemap = self._load_eles(rallocs, mesh, initsoln)
File "/develop/home/antonio/Software_Local/pyFR/pyfr/solvers/base/system.py", line 66, in _load_eles
elemap[t] = self.elementscls(basismap[t], mesh[f], self._cfg)
File "/develop/home/antonio/Software_Local/pyFR/pyfr/solvers/base/elements.py", line 53, in __init__
self._gen_pnorm_fpts()
File "/develop/home/antonio/Software_Local/pyFR/pyfr/solvers/base/elements.py", line 196, in _gen_pnorm_fpts
smats = self._get_smats(self._basis.fpts).transpose(1, 3, 0, 2)
File "/develop/home/antonio/Software_Local/pyFR/pyfr/solvers/base/elements.py", line 228, in _get_smats
jac = self._get_jac_eles_at(pts)
File "/develop/home/antonio/Software_Local/pyFR/pyfr/solvers/base/elements.py", line 216, in _get_jac_eles_at
jacop = np.rollaxis(self._basis.sbasis.jac_nodal_basis_at(pts), 2)
File "/develop/home/antonio/Software_Local/pyFR/pyfr/nputil.py", line 28, in newfn
arr = fn(*args, **kwargs)
File "/develop/home/antonio/Software_Local/pyFR/pyfr/polys.py", line 50, in jac_nodal_basis_at
return np.linalg.solve(self.vdm, self.jac_ortho_basis_at(epts))
File "/usr/lib/python2.7/dist-packages/numpy/linalg/linalg.py", line 311, in solve
_assertRank2(a, b)
File "/usr/lib/python2.7/dist-packages/numpy/linalg/linalg.py", line 155, in _assertRank2
'two-dimensional' % len(a.shape))
numpy.linalg.linalg.LinAlgError: 3-dimensional array given. Array must be two-dimensional

Hi Antonio,

My psychic debugging skills suggest that the version of NumPy you are
using is older than 1.8. PyFR requires NumPy 1.8 or later to function.

Regards, Freddie.

Hi Freddie,

I confirm your phyco debugging skills, :), thanks a lot once more.

Btw, have you ever run the Euler case Ringleb Flow with pyFR?

Best regards,
Antonio

Hi Antonio,

I have no experience with this particular test case.

Regards, Freddie.