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