[Getdp] Dual Formulation in the wiki

Amit.Itagi at seagate.com Amit.Itagi at seagate.com
Mon Jan 29 20:04:18 CET 2007


Hi,

I have a question regarding the dual formulation example (
http://www.geuz.org/getdp/wiki/DualFormulations ). I find that the first
order and the higher order elements satisfy different constraints (See
excerpt below). Could someone kindly explain the motivation for defining
separate constraints ?

The higher order elements are being used as a result of P1 adaptation. I
would think that since the same problem is being solved, the constraints
would remain the same.

Thanks

Rgds,
Amit


//----------------------------------------------------------------------------------------------------------------------------
FunctionSpace {
  { Name Hgrad_phi_MagSta; Type Form0;
    BasisFunction {
      { Name sn1n ; NameOfCoef wn1n ; Function BF_Node ; Support Tot_Mag;
Entity NodesOf[All]; }
      { Name sn2e ; NameOfCoef wn2e ; Function BF_Node_2E ; Support Tot_Mag
; Entity EdgesOf[ All ] ; }
      .
       .
       .
    }
    Constraint {
      { NameOfCoef wn1n; EntityType NodesOf ; NameOfConstraint Sta_phi; }
      { NameOfCoef wn2e ; EntityType EdgesOf ; NameOfConstraint
Sta_phi_Zero ; }
      .
      .
      .
    }
  }

Constraint {
  { Name Sta_phi ;
    Case {
      { Region left ; Value -100. ; }
      { Region right ; Value 100. ; }
    }
  }
  { Name Sta_phi_Zero ;
    Case {
      { Region left ; Value 0. ; }
      { Region right ; Value 0. ; }
    }
  }
//---------------------------------------------------------------------------------