PRISMS-PF  v2.1
initialConditions.cc
Go to the documentation of this file.
1 /*
2 //initial condition
3 template <int dim>
4 double InitialCondition<dim>::value (const dealii::Point<dim> &p, const unsigned int component) const
5 {
6  double scalar_IC=0;
7  // =====================================================================
8  // ENTER THE INITIAL CONDITIONS HERE FOR SCALAR FIELDS
9  // =====================================================================
10  // Enter the function describing conditions for the fields at point "p".
11  // Use "if" statements to set the initial condition for each variable
12  // according to its variable index.
13 
14  // =====================================================================
15  return scalar_IC;
16 }
17 
18 template <int dim>
19 void InitialConditionVec<dim>::vector_value (const dealii::Point<dim> &p, dealii::Vector<double> &vector_IC) const
20 {
21  // =====================================================================
22  // ENTER THE INITIAL CONDITIONS HERE FOR VECTOR FIELDS
23  // =====================================================================
24  // Enter the function describing conditions for the fields at point "p".
25  // Use "if" statements to set the initial condition for each variable
26  // according to its variable index.
27 
28 
29  // =====================================================================
30 }
31 */