SfePy April 2016

sfepy@python.org

5 participants
10 discussions

by Robert Cimrman

In examples/large_deformation/hyperelastic.py a rotation by displacements is applied. By using a similar function the vectors defining the force couples could be defined for dw_surface_ltr (IMHO). Does it make sense? r. ----- Reply message ----- From: "Andre Smit" <freev...(a)gmail.com> To: <sfepy...(a)googlegroups.com> Subject: Torque Date: Sat, Dec 18, 2010 05:10 What is the best way to apply a torque load to a model? -- Andre -- You received this message because you are subscribed to the Google Groups "sfepy-devel" group. To post to this group, send email to sfepy...(a)googlegroups.com. To unsubscribe from this group, send email to sfepy-devel...(a)googlegroups.com. For more options, visit this group at http://groups.google.com/group/sfepy-devel?hl=en.

2 years, 12 months

new interactive example with time stepping

by Robert Cimrman

Hi, I have just updated the time stepping solvers in sfepy for interactive use, as demonstrated in the new example [1]. For basic use, ignore the probing code - the time stepper can be used as simply as: tss = SimpleTimeSteppingSolver({'t0' : 0.0, 't1' : 100.0, 'n_step' : 11}, problem=problem) tss.init_time() for step, time, state in tss(): pass r. [1] http://sfepy.org/doc-devel/examples/diffusion/time_poisson_interactive.html

3 years, 11 months

Need help to start with sfepy

by Nicolas DELAN

Hi everyone, I am a new user of sfepy, i need it to compute the diffusions equations on a turbine blade. But i am a little lost with the boundary conditions. I have tried to impose some dirichlet boundary conditions on some points of the edges of the airfoil but it keeps crashing with exceptions. For example in my .mesh file i have this at the line where my point is defined 0.1 0.5 25 (x) (y) So I have written this in the region definition regions = { 'Omega' : 'all', 'point25' : ('vertices of group 25', 'edge'), } I have tried before with : refions = { 'Omega' : 'all', 'extrados' : ('vertices in (y>0.1) ', 'facet'), 'intrados' : ('vertices in (y<0)', 'facet'), } and it was giving resuls (it was only for test purpose) I don't understand very well how to define the boundary conditions (and to be honest the hole configuration file despite the documentation) So if somebody could help me, i would be grateful ! cordially, Nicolas DELAN

4 years, 9 months

postprocessing

by Nina Halabi

Dear Robert I could overcome my previous problems of not having root privileges, so it all worked up till the command: ./postproc.py cylinder.vtk This echoed the below message: Exception In /home/gmh/.local/lib/python2.7/site-packages/mayavi-4.4.4.dev0-py2.7-linux-x86_64.egg/mayavi/sources/vtk_xml_file_reader.py:86 AttributeError: 'AlgorithmOutput' object has no attribute 'point_data' (in get_all_attributes) Traceback (most recent call last): File "/home/gmh/.local/lib/python2.7/site-packages/traits-4.6.0.dev354-py2.7-linux-x86_64.egg/traits/trait_notifiers.py", line 340, in __call__ self.handler( *args ) File "/home/gmh/sfepy/sfepy/postprocess/dataset_manager.py", line 177, in _dataset_changed self._assign_attribute.input = value File "/home/gmh/.local/lib/python2.7/site-packages/traits-4.6.0.dev354-py2.7-linux-x86_64.egg/traits/trait_handlers.py", line 104, in _read_only name, class_of( object ) ) TraitError: The 'input' trait of an AssignAttribute instance is 'read only'. sfepy: point scalars t at [-0.05 -0.02 0. ] sfepy: range: -2.00e+00 2.00e+00 l2 norm range: 1.06e-02 2.00e+00 No handlers could be found for logger "mayavi.core.common" I'm using Python 2.7.8 if that's relevant to the current issue. Thanks again Cheers

4 years, 10 months

Region definition

by Dimitri Balasoiu

Hello everyone, I have a question about sfepy's definition of regions. I want to define a region Omega_1, wich is a subregion of Omega. So I type : regions = { 'Omega' : 'all', 'Omega_1' : ('vertices of group 1', 'cell', Omega),} , where group 1 is a group of vertices. Question : wich field element P1 will be created over the region Omega_1 ? More precisely, I have a triangle (A,B,C), with node A in group 1, but the other ones outside of group 1. Will the finite element wich is worth 1 at A, and 0 at B and C be in my field element basis ? In wich file are the basis functions for the fe fields defined ? Dimitri

4 years, 10 months

Re: [sfepy-devel] install sfepy vcvarshell error

by Robert Cimrman

Hi Mark, yes this error is annoying, and I am not aware of a solution that actually works for sfepy (failed myself to install sfepy in a W10 virtual box). I assume you tried the solutions proposed for example in [1, 2, 3]? I do not think that I have tried everything, so in case there is something you did not try, could you give it a shot? r. [1] https://support.enthought.com/hc/en-us/articles/204469210-Windows-Unable-to… [2] https://stackoverflow.com/questions/24380442/getting-error-unable-to-find-v… [3] https://github.com/cython/cython/wiki/CythonExtensionsOnWindows On 04/07/2016 05:36 PM, 'Mark Shart' via sfepy-devel wrote: > Hi everyone, > > I try to install sfepy for window10 following the instruction with > python(xy). > I can accomplish everything till I need to compile sfepy. The I get the > error:Unable to find vcvarsall.bat. > I installed the MS C++ Compiler for Python 2.7 but its not working. I also > tried other C++ compiler... but still this error drives me crazy. > How can I fix this? Please help :\ >

4 years, 10 months

issues #317, #318

by Robert Cimrman

Hi, any windows users around that could help with the issues #317, #318? r. https://github.com/sfepy/sfepy/issues/317 https://github.com/sfepy/sfepy/issues/318

4 years, 10 months

Install error Windows10

by Mark Shart

Hi everyone, I'm trying to install sfepy on windows10 with python(xy) as suggested, but get the following error when compiling python setup.py build_ext --inplace --compiler=mingw32 $ python setup.py build_ext --inplace igakit (optional) is missing; you may get run-time errors petsc4py (optional) is missing; you may get run-time errors mpi4py (optional) is missing; you may get run-time errors pymetis (optional) is missing; you may get run-time errors pysparse (optional) is missing; you may get run-time errors running build_ext running build_src build_src building library "sfepy_common" sources building library "sfepy_terms" sources building extension "sfepy.discrete.common.extmods._fmfield" sources building extension "sfepy.discrete.common.extmods.mappings" sources building extension "sfepy.discrete.common.extmods.assemble" sources building extension "sfepy.discrete.common.extmods.cmesh" sources building extension "sfepy.discrete.common.extmods.crefcoors" sources building extension "sfepy.discrete.common.extmods._geommech" sources building extension "sfepy.discrete.fem.extmods.bases" sources building extension "sfepy.discrete.fem.extmods.lobatto_bases" sources building extension "sfepy.discrete.iga.extmods.igac" sources building extension "sfepy.linalg.extmods.crcm" sources building extension "sfepy.terms.extmods.terms" sources building data_files sources build_src: building npy-pkg config files customize Mingw32CCompiler customize Mingw32CCompiler using build_clib customize Mingw32CCompiler customize Mingw32CCompiler using build_ext extending extension 'sfepy.discrete.common.extmods._fmfield' defined_macros with [('__SDIR__', '\'"sfepy\\discrete\\common\\extmods"\''), ('SFEPY_PLATFORM', 1)] extending extension 'sfepy.discrete.common.extmods.mappings' defined_macros with [('__SDIR__', '\'"sfepy\\discrete\\common\\extmods"\''), ('SFEPY_PLATFORM', 1)] extending extension 'sfepy.discrete.common.extmods.crefcoors' defined_macros with [('__SDIR__', '\'"sfepy\\discrete\\common\\extmods"\''), ('SFEPY_PLATFORM', 1)] extending extension 'sfepy.discrete.common.extmods._geommech' defined_macros with [('__SDIR__', '\'"sfepy\\discrete\\common\\extmods"\''), ('SFEPY_PLATFORM', 1)] extending extension 'sfepy.discrete.fem.extmods.bases' defined_macros with [('__SDIR__', '\'"sfepy\\discrete\\common\\extmods"\''), ('SFEPY_PLATFORM', 1)] extending extension 'sfepy.discrete.fem.extmods.lobatto_bases' defined_macros with [('__SDIR__', '\'"sfepy\\discrete\\common\\extmods"\''), ('SFEPY_PLATFORM', 1)] extending extension 'sfepy.terms.extmods.terms' defined_macros with [('__SDIR__', '\'"sfepy\\terms\\extmods"\''), ('SFEPY_PLATFORM', 1), ('__SDIR__', '\'"sfepy\\discrete\\common\\extmods"\'')] building 'sfepy.discrete.common.extmods._fmfield' extension compiling C sources C compiler: gcc -O2 -Wall -Wstrict-prototypes compile options: '-D__SDIR__='"sfepy\discrete\common\extmods"' -DSFEPY_PLATFORM=1 -D__SDIR__='"sfepy\discrete\common\extmods"' -DSFEPY_PLATFORM=1 -DNPY_MINGW_USE_CUSTOM_MSVCR -D__MSVCRT_VERSION__=0x0900 -Isfepy\discrete\common\extmods -IC:\Python27\lib\site-packages\numpy\core\include -IC:\Python27\include -IC:\Python27\PC -c' extra options: '-g -O2' gcc -O2 -Wall -Wstrict-prototypes -D__SDIR__='"sfepy\discrete\common\extmods"' -DSFEPY_PLATFORM=1 -D__SDIR__='"sfepy\discrete\common\extmods"' -DSFEPY_PLATFORM=1 -DNPY_MINGW_USE_CUSTOM_MSVCR -D__MSVCRT_VERSION__=0x0900 -Isfepy\discrete\common\extmods -IC:\Python27\lib\site-packages\numpy\core\include -IC:\Python27\include -IC:\Python27\PC -c sfepy\discrete\common\extmods\_fmfield.c -o build\temp.win32-2.7\Release\sfepy\discrete\common\extmods\_fmfield.o -g -O2 Found executable C:\MinGW\bin\gcc.exe g++ -shared build\temp.win32-2.7\Release\sfepy\discrete\common\extmods\_fmfield.o -LC:\Python27\libs -LC:\Python27\PCbuild -Lbuild\temp.win32-2.7 -lsfepy_common -lpython27 -lmsvcr90 -o sfepy\discrete\common\extmods\_fmfield.pyd Found executable C:\MinGW\bin\g++.exe C:\Python27\libs/libpython27.a: error adding symbols: File format not recognized collect2.exe: error: ld returned 1 exit status error: Command "g++ -shared build\temp.win32-2.7\Release\sfepy\discrete\common\extmods\_fmfield.o -LC:\Python27\libs -LC:\Python27\PCbuild -Lbuild\temp.win32-2.7 -lsfepy_common -lpython27 -lmsvcr90 -o sfepy\discrete\common\extmods\_fmfield.pyd" failed with exit status 1 Can anyone help me please?

4 years, 11 months

Re: [sfepy-devel] Obtain the required data from a *.vtk file?

by Robert Cimrman

On 04/02/2016 02:50 PM, Kid Guo wrote: > Thank you for your reply. > Last time I did not have a good description of my problem,but I made some > changes on your script and got some useful results for me. > > Now I have a new question,I have got all points(such as,x=0.1) > and corresponding cells,if you run the modified script,you can see it, > there are 11 cells on the face where x=0.1, OK, so you want cells containing faces that are in x=0.1 plane, right? I have modified the script to find those - there are 16 of them - there was a mistake in your script, as your cells have 4 vertices, not only three. > I want to get the values(p and u) of the four points of each unit instead > of all cells. See the attachment. r. > Any ideas? > > Thank you advance. > > 在 2016年4月1日星期五 UTC+8下午11:07:34，Robert Cimrman写道： >> >> On 04/01/2016 05:08 AM, Kid Guo wrote: >>> I am trying to Obtain the required data from a *.vtk file. >>> For example, >>> 1. I want to output the all points where{x=0.1},such as (0.1,0.1,0.1) >>> 2. I want to obtain the CELLS including points in step1;such as >>> (12,13,56,23) >>> 3. I want to obtain the values(p,u) in the *.vtk file >> >> The attached script should show you how to do what you want. >> >>> < >> https://lh3.googleusercontent.com/-edlklG7xlrM/Vv3mDxDgzDI/AAAAAAAAACo/BajI…> >> >> >> This seems to correspond to an old version of sfepy - the script is w.r.t. >> the >> latest version. >> >>> I am using ipython,but I have no idea. >> >> Yeah, these low-level features are not extensively documented - one has to >> know >> what to look for. >> >> r. >> >>> Here is the attachment >>> >>> >>> What should I do? >>> >>> Thank you advance. >> >> >

4 years, 11 months

Obtain the required data from a *.vtk file?

by Kid Guo

I am trying to Obtain the required data from a *.vtk file. For example, 1. I want to output the all points where{x=0.1},such as (0.1,0.1,0.1) 2. I want to obtain the CELLS including points in step1;such as (12,13,56,23) 3. I want to obtain the values(p,u) in the *.vtk file <https://lh3.googleusercontent.com/-edlklG7xlrM/Vv3mDxDgzDI/AAAAAAAAACo/BajI…> I am using ipython,but I have no idea. Here is the attachment What should I do? Thank you advance.

4 years, 11 months

2021

2020

2019

2018

2017

2016

2015

2014

2013

2012

2011

2010

2009

2008

SfePy April 2016