std::vector<state> flic(const std::vector<state> & X, int BCL, int BCR, double dx, double C, double t, int l) { double t0 = 0; int iter = 0; std::vector<state> ret = X; while (t0 < t) { std::vector<state> X0 = add_boundary_conditions(ret, BCL, BCR, 2); double dt = timestep(X0, dx, C, t0, t, iter); flic_stepper(ret, X0, dx, dt, C, l); t0 += dt; iter++; } return ret; }
void EssentialBCs<Scalar>::add_boundary_condition(EssentialBoundaryCondition<Scalar> * boundary_condition) { Hermes::vector<EssentialBoundaryCondition<Scalar> *> boundary_conditions; boundary_conditions.push_back(boundary_condition); add_boundary_conditions(boundary_conditions); }
EssentialBCs<Scalar>::EssentialBCs(EssentialBoundaryCondition<Scalar> * boundary_condition) : HermesAnyBC(NULL) { Hermes::vector<EssentialBoundaryCondition<Scalar> *> boundary_conditions; boundary_conditions.push_back(boundary_condition); add_boundary_conditions(boundary_conditions); }
EssentialBCs<Scalar>::EssentialBCs(Hermes::vector<EssentialBoundaryCondition<Scalar> *> essential_bcs) : HermesAnyBC(NULL) { add_boundary_conditions(essential_bcs); }
EssentialBCs::EssentialBCs(EssentialBoundaryCondition * boundary_condition) { Hermes::vector<EssentialBoundaryCondition *> boundary_conditions; boundary_conditions.push_back(boundary_condition); add_boundary_conditions(boundary_conditions); };
EssentialBCs::EssentialBCs(Hermes::vector<EssentialBoundaryCondition *> essential_bcs) { add_boundary_conditions(essential_bcs); };