Robert, Ondrej,
how about f2c to turn dft.f into dft.c you just need f2c.h.
I just edited the .f file, moved the integer declaration to above the
"data" statement and ran f2c.
gcc -o dft.so --shared -fPIC dft.c created the .so file.
shroedinger.py still seems to work.
br.
-osman
Hi Robert, Ondrej,
Finally I was able to re-build numpy/scipy and tried the new version of
sfepy.
Even during building of scipy, whenever a fortran or f2py based shared
lib to be built, I need to manually insert "-shared" to the command
(that is after it fails to build an executable! . Same with sfepy
physics module. But tests did run and two failed :
...
tests/test_lcbc_2d.py
--- test_linear_rigid_body_bc: failed!
tests/test_lcbc_3d.py
--- test_linear_rigid_body_bc: failed!
...
expected?
Regards,
-osman
I am pleased announce the release of SfePy 00.50.00.
SfePy is a finite element analysis software in Python, based primarily
on Numpy and SciPy.
Mailing lists, issue tracking, mercurial repository: http://sfepy.org
Home page: http://sfepy.kme.zcu.cz
People who contributed to this release: Ondrej Certik, Ryan Krauss,
Vladimir Lukes.
Major improvements:
- finite strain elasticity: neo-Hookean, Mooney-Rivlin materials in the
total Lagrangian (TL) formulation
- solving problems in complex numbers
- generalized equations to allow linear combination of terms
- run-time type of state term arguments
- refactoring to follow Python coding style guidelines
- new terms
For more information on this release, see
http://sfepy.googlecode.com/svn/web/releases/005000_RELEASE_NOTES.txt
Best regards,
Robert Cimrman
Hi all, and Ryan in particular,
I have started implementing large deformations into the code. Nothing is
finished yet, but I have already some proof-of-concept.
I will commit relatively soon (i.e. after my holidays) a patch adding
the hyperelastic neo-Hookean material.
Ultimately, I would like to have in SfePy all the functionality I had in
its predecessor (a matlab code called mafest1), but it is a long term
goal, as it was more a one-purpose tool than a general FE system.
All that the old code can do is nicely summarized in my thesis, that you
can find at
http://ui505p06-mbs.ntc.zcu.cz/sfe/RobertCimrman
at the bottom (thesis.pdf).
Ryan, try to look in particular at Section 2.1.2 (TL formulation),
Chapter 3 (FE discretization in TL context) and Chapter 6 (FE
discretization in a linear problem context), as it may provide you some
insight. Appendix A might be also interesting, as the usual constitutive
(stress-strain) relations are summarized there, e.g. the neo-Hookean one.
r.
Hi,
I'd like to be able to do these things with sfepy+sympy together:
http://reference.wolfram.com/mathematica/tutorial/NDSolvePDE.html
I.e. really just specify the equation in sympy and sfepy should be
clever enough to solve it correctly, or say why it cannot solve it.
Ondrej
