BiotSavart::BiotSavart(Mesh &ms, const Vect<complex_t> &J, Vect<complex_t> &B, int code)
{
_mu = 4*OFELI_PI*1.e-7;
_theMesh = &ms;
_code = code;
setSurface();
_BC = &B;
_C = true;
_bound = false;
_type = J.getDOFType();
_theEdgeList = NULL;
_theSideList = NULL;
_JC = &J;
if (_type==SIDE_FIELD)
_theSideList = new SideList(*_theMesh);
if (_type==EDGE_FIELD)
_theEdgeList = new EdgeList(*_theMesh);
}
void BiotSavart::setCurrentDensity(const Vect<complex_t> &J)
{
_JC = &J;
_C = true;
_type = J.getDOFType();
}
void BiotSavart::setCurrentDensity(const Vect<real_t> &J)
{
_J = &J;
_C = false;
_type = J.getDOFType();
}
