// initialize the Neumann vector void setNeumann(FVMesh2D &m,FVVect<double> &Vn,Parameter ¶) { FVEdge2D *ptr_e; m.beginEdge(); while((ptr_e=m.nextEdge())) { Vn[ptr_e->label-1]=ptr_e->getMeanValue(Neumann,para); } }
// initialize the Dirichelt vector void setDirichlet(FVMesh2D &m,FVVect<double> &Vd,Parameter ¶) { FVEdge2D *ptr_e; m.beginEdge(); while((ptr_e=m.nextEdge())) { Vd[ptr_e->label-1]=ptr_e->getMeanValue(Dirichlet,para); } }
// initialize the velocity void setVelocity(FVMesh2D &m,FVVect<FVPoint2D<double> > &u,Parameter ¶) { FVEdge2D *ptr_e; m.beginEdge(); while((ptr_e=m.nextEdge())) { u[ptr_e->label-1].x=Xvelocity(ptr_e->centroid,para); u[ptr_e->label-1].y=Yvelocity(ptr_e->centroid,para); } }