[Libmesh-users] error in FE reinit

[Manav B. <manav@u.washington.edu>]

Hi
    I am getting an error of negative-jacobian from this code.
    Could anyone help me with this please?


                    Node point1(.2308846387,.4,.3773333333,0),
                      point2(0,0,.4,0);
                    Edge2 edge;
                    edge.set_node(0) = &point1;
                    edge.set_node(1) = &point2;

                    std::auto_ptr<FEBase> fe(new FE<1,LAGRANGE>(FEType 
()));
                    std::auto_ptr<QBase> qbase(new QGauss(1,FIFTH));
                    fe->attach_quadrature_rule(qbase.get());
                    fe->reinit(&edge);



Thanks
Manav