Re: [Getfem-users] nonlinear brick

[Roman Putanowicz]

Dear Yves,
> 
> Your term
>  double a = vecl[0][0];
> cannot be correct. vecl is to be filled by  asm_real_tangent_terms. It is 
> initialized with 0.
Thank you a lot.
This is the point I was looking for. With your remark I was able to 
write  a code which solves the equation a*x^2+b*x+c = 0.
When I clean it a bit I will post it to the list.

Writing such program may seem strange at first, but for me such program
gives the opportunity to play with some of GetFEM++ concepts without the 
whole baggage of FEM. What is really appealing to me is explicit or
almost explicit representation of mathematical model in GetFEM.

This "explicitness" manifests itself  for instance in the possibility
to enhance bricks and the model classes so they could generate a LaTeX 
math expressions for the variational form. 

Regards,

Roman