Example #1
0
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);
 }
Example #5
0
EssentialBCs::EssentialBCs(EssentialBoundaryCondition * boundary_condition) {
  Hermes::vector<EssentialBoundaryCondition *> boundary_conditions;
  boundary_conditions.push_back(boundary_condition);
  add_boundary_conditions(boundary_conditions);
};
Example #6
0
EssentialBCs::EssentialBCs(Hermes::vector<EssentialBoundaryCondition *> essential_bcs) {
  add_boundary_conditions(essential_bcs);
};