[Getdp] [Gmsh] ¿how to express a weak formulation for an integral operator?

Luis Saavedra luis94855510 at gmail.com
Thu Mar 13 09:59:09 CET 2008


Dear Community,

i'm attempting to express this week formulation (in LaTeX notation):

\Int_{\Omega} T(x,y) v(x,y) dxdy = 
\Int_{\Omega}\Int_{\Omega}(1/r).(grad' u(x',y'))(grad v(x,y)) dxdydx'dy'

where:

u = u(x',y') is the unknown function,
v = v(x,y) is a test function,
T is a known function,
r = || (x-x')^2 + (y-y')^2 ||, and
grad' is the gradient with respect to de x', y' variables.

it would be something as this:

Formulation {
  {Name WeekIEM; Type FemEquation;
    Quantity {
      { Name  v; Type Local;  NameOfSpace  TestSpace; }
    }
    Equation {
      Galerkin { [[ (1/r) , Dof{Grad v} ], {Grad v}] ; In Omega ; 
Jacobian Vol ; Integration Int ; }
      Galerkin { [ -T[],{v}] ; In Omega ; Jacobian Vol ; Integration Int ; }
    }
}

but I know that this will not work.

I like very much the free software GetDP, is my made dream reality, i'm 
very grateful with you.

Thanks for all and best regards,
Luis Saavedra.

Santiago Chile.