Namespace for exact solution and problem parameters. More...

Functions
void	source (const Vector< double > &x, Vector< double > &f)
	Source function. More...

void	mass_source (const Vector< double > &x, double &m)
	Mass source. More...

void	boundary_pressure (const double &time, const Vector< double > &x, const Vector< double > &n, double &result)
	Pressure around the boundary of the domain. More...

void	exact_soln (const Vector< double > &x, Vector< double > &soln)
	Exact solution: q1,q2,div_q,p. More...

template<class ELEMENT >
void	edge_sign_setup (Mesh *mesh_pt)

Variables
bool	Use_point_source_solution =false

double	Mass_source_strength =1000.0
	Mass source strength. More...

double	Mass_source_decay_exponent =50.0
	Mass source strength. More...

double	Element_area =0.01
	Target element area for Triangle. More...

std::string	Directory ="RESLT"
	The directory in which the solution is output. More...

Detailed Description

Namespace for exact solution and problem parameters.

Function Documentation

◆ boundary_pressure()

void TestProblem::boundary_pressure	(	const double &	time,
		const Vector< double > &	x,
		const Vector< double > &	n,
		double &	result
	)

Pressure around the boundary of the domain.

   {
    if (Use_point_source_solution)
     {
      result=0.0;
     }
    else
     {
      result=2*x[0]+3*x[1]-3.0/2.0;
     }
   }

References Use_point_source_solution, and plotDoE::x.

Referenced by DarcyProblem< ELEMENT >::create_pressure_elements().

◆ edge_sign_setup()

template<class ELEMENT >

void TestProblem::edge_sign_setup ( Mesh * mesh_pt )

Global function that completes the edge sign setup – has to be called before projection in unstructured adaptation

  {
   // The dictionary keeping track of edge signs
   std::map<Edge,unsigned> assignments;
   
   // Loop over all elements
   unsigned n_element = mesh_pt->nelement();
   for(unsigned e=0;e<n_element;e++)
    {
     ELEMENT* el_pt = dynamic_cast<ELEMENT*>(mesh_pt->element_pt(e));
   
     // Assign edge signs: Loop over the vertex nodes (always 
     // first 3 nodes for triangles)
     for(unsigned i=0;i<3;i++)
      {
       Node *nod_pt1, *nod_pt2;
       nod_pt1 = el_pt->node_pt(i);
       nod_pt2 = el_pt->node_pt((i+1)%3);
       Edge edge(nod_pt1,nod_pt2);
       unsigned status = assignments[edge];
       
       // This relies on the default value for an int being 0 in a map
       switch(status)
        {
         // If not assigned on either side, give it a + on current side
        case 0:
         assignments[edge]=1;
         break;
         // If assigned + on other side, give it a - on current side
        case 1:
         assignments[edge]=2;
         el_pt->sign_edge(i)=-1;
         break;
         // If assigned - on other side, give it a + on current side
        case 2:
         assignments[edge]=1;
         break;
        }
      } // end of loop over vertex nodes
    } // end of loop over elements
  
  
  }

References anonymous_namespace{skew_symmetric_matrix3.cpp}::assignments(), e(), oomph::Mesh::element_pt(), i, and oomph::Mesh::nelement().

◆ exact_soln()

void TestProblem::exact_soln	(	const Vector< double > &	x,
		Vector< double > &	soln
	)

Exact solution: q1,q2,div_q,p.

   {
  
    if (Use_point_source_solution)
     {
      // q[0]
      soln[0]=0.0;
      // q[1]
      soln[1]=0.0;
      // div q
      soln[2]=0.0;
      // p
      soln[3]=0.0;
     }
    else
     {
      // q[0]
      soln[0]=x[0]*x[1]-x[1]*x[1];
      // q[1]
      soln[1]=x[0]+x[0]*x[0]-0.5*x[1]*x[1];
      // div q
      soln[2]=0;
      // p
      soln[3]=2*x[0]+3*x[1]-3.0/2.0;
     }
   }

References Use_point_source_solution, and plotDoE::x.

Referenced by DarcyProblem< ELEMENT >::complete_problem_setup(), and DarcyProblem< ELEMENT >::doc_solution().

◆ mass_source()

void TestProblem::mass_source	(	const Vector< double > &	x,
		double &	m
	)

Mass source.

   {
    if (Use_point_source_solution)
     {
      double r=sqrt(x[0]*x[0]+x[1]*x[1]);
      m=Mass_source_strength*exp(-r*Mass_source_decay_exponent);
     }
    else
     {
      m=0.0;
     }
   }

References Eigen::bfloat16_impl::exp(), m, Mass_source_decay_exponent, Mass_source_strength, UniformPSDSelfTest::r, sqrt(), Use_point_source_solution, and plotDoE::x.

Referenced by DarcyProblem< ELEMENT >::complete_problem_setup().

◆ source()

void TestProblem::source	(	const Vector< double > &	x,
		Vector< double > &	f
	)

Source function.

   {
    if (Use_point_source_solution)
     {
      f[0]=0.0;
      f[1]=0.0;
     }
    else
     {
      f[0]=x[0]*x[1]-x[1]*x[1]+2;
      f[1]=x[0]+x[0]*x[0]-0.5*x[1]*x[1]+3;
     }
   }

References f(), Use_point_source_solution, and plotDoE::x.

Variable Documentation

◆ Directory

std::string TestProblem::Directory ="RESLT"

The directory in which the solution is output.

Referenced by DarcyProblem< ELEMENT >::DarcyProblem(), DarcyProblem< ELEMENT >::doc_shape_functions(), DarcyProblem< ELEMENT >::doc_solution(), and main().

◆ Element_area

double TestProblem::Element_area =0.01

Target element area for Triangle.

Referenced by DarcyProblem< ELEMENT >::DarcyProblem(), DarcyProblem< ELEMENT >::doc_solution(), and main().

◆ Mass_source_decay_exponent

double TestProblem::Mass_source_decay_exponent =50.0

Mass source strength.

Referenced by mass_source().

◆ Mass_source_strength

double TestProblem::Mass_source_strength =1000.0