My complete problem is actually a dynamic problem, at every time-step run sfepy as stress-strain solver. However, when the mesh is large such as 513*513, sfepy runs slowly, especially during solving equlibrium equation from iter 0 to iter 1. The solving process takes roughly 1 min/time-step. I am wondering is there any good idea to improve the running speed?

在 2015年3月5日星期四 UTC+1下午5:37:00，Robert Cimrman写道：

On 03/05/2015 05:01 PM, Ronghai Wu wrote:
>
> 在 2015年3月3日星期二 UTC+1下午11:35:05，Robert Cimrman写道：

<snip>

>>> Stress is plotted in this way. A line along x-axis(or y-axis) across the
>>> center of hole, plot from the center of hole to boundary, distance(r) is
>>> normalized by hole radius(R).
>>
>> BTW. you can get the stress by using the ev_cauchy_stress term, just like
>> the
>> strain... also check the primer for data probes.
>>
>>
> I tried to do this, but got error which I cannot fix (please run
> "elasticity-periodicBC3.py"). Additionally, I have some questions regarding
> the strain and stress calculated by pb.evalate(please see "questions about
> strain and stress.pdf").

You just need to use the probing code from
examples/linear_elasticity/its2D_interactive.py, because your script is also
interactive. Attached is the fixed script (+ some additional small fixes).

As for the questions, the strain and stress coming from ev_cauchy_*() terms are
without the prestrain. You need to subtract that yourself.

r.