[Getfem-users] Outsourcing calculation of stiffness matrix

[Roman Putanowicz]

Dear All, 

I am working on the following problem:

I would like to be able to "outsource" calculation of
element matrices and vectors outside GetFEM++.
In  other words I would like to keep all GetFEM++ machinery
but during the assembly force GetFEM++ to call external 
function that will provide the matrices to be assembled.
The whole point is to be able to write the function in 
a scripting language for an interpreter embedded in the
application (I aim at Python).

I __do not__ want to pre compute element matrices and 
pass them as data to assembly routines.

The solution I see is to write a special class
that inherits getfem::nonlinear_elem_term and use
the ability to perform arbitrary calculations during
assembly via NonLin term.

In fact I envisage writing a special brick that 
will initialize the external routines (either
load shared library or load the scripts)
and during the assembly will use mentioned 
special nonlinear term.

The first attempts are promising but I would like
to ask for your opinions on this idea. Have anybody
tried this? Are there any subtle points to watch out?
I will appreciate any comments.

Regards

Roman
-- 
Roman Putanowicz, PhD  < address@hidden  >
Institute for Computational Civil Engng (L-5)
Dept. of Civil Engng, Cracow Univ. of Technology
www.l5.pk.edu.pl, tel. +48 12 628 2569, fax 2034