[Getdp] PML example

Christophe Geuzaine Christophe.Geuzaine at ulg.ac.be
Thu Dec 20 13:07:09 CET 2001


Christophe Geuzaine wrote:
> 
> c[] = 340 ;
> cx[] = Complex[1,-3] ;
> detJ [ pmlx ] = cx[] ;
> invJ [ pmlx ] = 1/TensorDiag[cx[],1,1] ;
> 
> FunctionSpace {
>   { Name Hgrad_p ; Type Form0 ;
>     BasisFunction {
>       { Name sn ; NameOfCoef pn ; Function BF_Node ;
>         Support Region[{Vol, PML}] ; Entity NodesOf[ All ] ; }

'PML' should of course coincide with 'pmlx' in this example.

-- 
Christophe Geuzaine

Tel: 32 (0) 4 366 37 10    http://www.geuz.org
Fax: 32 (0) 4 366 29 10    mailto:Christophe.Geuzaine at ulg.ac.be