[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