Mailman 3 - SfePy - python.org

mailing list for issues
by Robert Cimrman 04 Feb '08

04 Feb '08

there is now sfepy-...(a)googlegroups.com r.

1 1

amd64 problems, scipy 0.6.0 patch
by Ondrej Certik 04 Feb '08

04 Feb '08

Hi, To make sfepy work with the latest released scipy, use this patch: diff --git a/sfe/solvers/ls.py b/sfe/solvers/ls.py --- a/sfe/solvers/ls.py +++ b/sfe/solvers/ls.py @@ -1,7 +1,7 @@ from sfe.base.base import * from sfe.base.base import * from sfe.solvers.solvers import LinearSolver -import scipy.splinalg.dsolve.umfpack as um +import scipy.linsolve.umfpack as um um.configure( assumeSortedIndices = True ) ## However, today I discovered problems when using sfepy on amd64 (Debian): It compiles fine, but: $ ./schroedinger.py input/schroed.py warning: other missing: ['functions', 'modules', 'epbc', 'lcbc'] warning: left over: ['funV', 'options'] mesh.read = 0.04 mesh.split = 0.0 t = 0.0 Traceback (most recent call last): File "./schroedinger.py", line 110, in ? main() File "./schroedinger.py", line 107, in main evp = solveEigenProblem( conf, options ) File "./schroedinger.py", line 24, in solveEigenProblem pb = ProblemDefinition.fromConf( conf ) File "/home/ondra/sfepy/sfe/fem/problemDef.py", line 40, in fromConf domain.setupNeighbourLists() File "/home/ondra/sfepy/sfe/fem/domain.py", line 349, in setupNeighbourLists mu.sortRows( aux, nm.array( sortCols[mode], nm.int32 ) ) File "/home/ondra/sfepy/sfe/fem/extmods/meshutils.py", line 102, in sortRows return _meshutils.sortRows(*args) TypeError: array cannot be safely cast to required type This is imho related to the fact, that the 64bit system uses integers of a different default length. Ondrej

2 16

paraview postprocessing and other stuff
by Ondrej Certik 04 Feb '08

04 Feb '08

Hi, I greatly recommend paraview for postprocessing. We should add this to our wiki too. I created a Debian package, accessible in my repository: http://debian.certik.cz/ However, more work is needed, before it can get to the main Debian. I have a question - do you know how to save the pipeline, so that when I fix my program, calculate the result again, I don't have to fiddle with mouse in paraview and simply see the same pipeline as before but with new data? Ideally, do you have some example script how to call paraview from python? Let's add all this information to the wiki. Our first issue: http://code.google.com/p/sfepy/issues/detail?id=1 And to make things moving, I created some more issues. Now I am working on incorporating my mesh conversion utilities + automatic mesh generation using tetgen, will send a patch soon. Ondrej

2 1

planet scipy
by Robert Cimrman 03 Jan '08

03 Jan '08

Hi all (=3), check http://planet.scipy.org! Having a sfepy blog would be nice. happy new year, r.

2 1