[Gmsh] gmsh <--> getdp

Michael Grant m.grant at ee.wits.ac.za
Mon Jun 27 11:22:21 CEST 2011


Dear All,

I have an equivalence problem: in gmsh I can generate a field through the
inspection interface.  I simply use MathEval over the area of interest, and my
field is completely described by the co-ordinates of the problem.

Is there an equivalent way of creating the same field in getdp?  I wish to
integrate this field within a region, and have attached my getdp problem
definition for reference. 

Kind regards,

Mike.

-----------------

Group {
    Space = Region[101]; Building = Region[102];     
    Problem = Region[ {Space, Building} ];

}

Function {
    prob[] = Exp(-1*($X + $Y);
}

FunctionSpace {
{ Name gspace; Type Scalar;
    BasisFunction {
{Name sn; NameOfCoef vn; Function BF_Node;
Support Problem; Entity NodesOf[All]; }
    }
    GlobalQuantity {
    {Name gquantity; Type AssociatedWith; NameOfCoef vn;}
    }
}
}

Formulation {
{ Name testf; Type FemEquation;
    Quantity {
    { Name gquantity; Type Global; NameOfSpace gspace; }
    }
    Equation {
         GlobalTerm {  [ prob[], {gquantity} ]  ; In Problem; }
    }
}
}

Resolution {
{ Name testr;
    System {
    {Name gsystem; NameOfFormulation testf;}
    }
    Operation {
    Generate [gsystem]; Solve[gsystem]; SaveSolution[gsystem];
    }
}
}

PostProcessing {
{ Name test; NameOfFormulation testf;
    Quantity {
    {Name gquantity; Value {Local{ [ {gquantity} ] ; In Problem; }}}
    }
}
}

PostOperation {
{Name post_out; NameOfPostProcessing test;
    Operation {
         Print [gquantity, OnElementsOf Problem, File  "map.pos"];
    }
}
} 
<html><p><font face = "verdana" size = "0.8" color = "navy">This communication is intended for the addressee only. It is confidential. If you have received this communication in error, please notify us immediately and destroy the original message. You may not copy or disseminate this communication without the permission of the University. Only authorized signatories are competent to enter into agreements on behalf of the University and recipients are thus advised that the content of this message may not be legally binding on the University and may contain the personal views and opinions of the author, which are not necessarily the views and opinions of The University of the Witwatersrand, Johannesburg. All agreements between the University and outsiders are subject to South African Law unless the University agrees in writing to the contrary.</font></p></html>