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
--
You received this message because you are subscribed to the Google Groups "PyFR
Mailing List" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to [email protected].
To post to this group, send an email to [email protected].
Visit this group at http://groups.google.com/group/pyfrmailinglist.
For more options, visit https://groups.google.com/d/optout.