Did you mean "linear_elastic_probes.html" instead of "linear_elastic_tractions.html" example? I found the "linear_elastic_probes.html" very useful example for my purposes to probe a solution in the given (x, y, z) points. Also the documentation "src/sfepy/fem/probes.html" gives all necessary information to help me implement what I want to do. Thank you!

But, despite this, could you tell me more about low-level way to evaluate a variable in the given (x, y, z) point?

Sincerely,

Alexander

On Sunday, August 26, 2012 12:47:50 AM UTC+4, Robert Cimrman wrote:

Hi Alec,

On 08/25/2012 05:45 PM, Alec Kalinin wrote:

> Dear SfePy users,

>

> Is it possible to evaluate a solution not only in the FEM mesh node, but in

> any arbitrary point in the domain with the given (x, y, z) coordinates?

Yes, it is possible. Either, you could use a probe as described in the

Primer [1] - the available probes are described in [2]. Or, you could

directly evaluate a variable in given points - this is a bit low-level

operation, but I could provide you instructions, if the probes are not

enough for you.

Cheers,

r.

> For example, consider Dirichlet problem for Poisson equation. We apply

> essential boundary conditions on the surface nodes and after the problem

> has been solved we have the solution vector, i.e. vector of values in the

> FEM mesh nodes. But I want to know the solution in point v(x, y, z) that is

> not FEM mesh node. What is the best way to obtain solution in this point v?

>

> Sincerely,

> Alec Kalinin

[1] doc-devel/examples/linear_elasticity/linear_elastic_ tractions.html

[2] http://sfepy.org/doc-devel/src/sfepy/fem/probes.html