GradientModel.h File Reference

#include "Model.h"
#include "Gradient.h"
#include "GradientMatrix.h"
#include "GeomFields.h"
#include "NumType.h"
#include "SquareTensor.h"
#include "Mesh.h"

Include dependency graph for GradientModel.h:

This graph shows which files directly or indirectly include this file:

Go to the source code of this file.

Classes
class	GradientModelBase
class	GradientModel< X >

Functions
template<class T >
void	reflectGradient (Gradient< T > &gr, const Gradient< T > &g0, const Vector< T, 3 > &en)
template<class T , int N>
void	reflectGradient (Gradient< Vector< T, N > > &gr, const Gradient< Vector< T, N > > &g0, const Vector< T, 3 > &en)
template<class T >
void	reflectGradient (Gradient< Vector< T, 2 > > &gr, const Gradient< Vector< T, 2 > > &g0, const Vector< T, 3 > &en)
template<class T >
void	reflectGradient (Gradient< Vector< T, 3 > > &gr, const Gradient< Vector< T, 3 > > &g0, const Vector< T, 3 > &en)

Function Documentation

template<class T >

void reflectGradient	(	Gradient< T > &	gr,
		const Gradient< T > &	g0,
		const Vector< T, 3 > &	en
	)

Definition at line 23 of file GradientModel.h.

Referenced by GradientModel< Vector< T, 3 > >::compute(), and MomentumPressureGradientDiscretization< X >::discretize().

{
  const T g0_dot_en_x2 = T(2.0)*(g0*en);
  for(int i=0; i<3; i++)
    gr[i] = g0[i] - g0_dot_en_x2*en[i];
}

template<class T , int N>

void reflectGradient	(	Gradient< Vector< T, N > > &	gr,
		const Gradient< Vector< T, N > > &	g0,
		const Vector< T, 3 > &	en
	)

Definition at line 34 of file GradientModel.h.

{
  for(int k=0; k<N; k++)
  {
      const T g0_dot_en_x2 = T(2.0)*(g0[0][k]*en[0] + g0[1][k]*en[1] + g0[2][k]*en[2]);
      for(int i=0; i<3; i++)
        gr[i][k] = g0[i][k] - g0_dot_en_x2*en[i];
  }
}

template<class T >

void reflectGradient	(	Gradient< Vector< T, 2 > > &	gr,
		const Gradient< Vector< T, 2 > > &	g0,
		const Vector< T, 3 > &	en
	)

Definition at line 51 of file GradientModel.h.

{
  const Vector<T,2> g0_dot_en_x2 = T(2.0)*(g0*en);
  for(int i=0; i<3; i++)
  {
      gr[i][0] = g0[i][0] - g0_dot_en_x2[0]*en[i];
      gr[i][1] = g0[i][1] - g0_dot_en_x2[1]*en[i];
  }
}

template<class T >

void reflectGradient	(	Gradient< Vector< T, 3 > > &	gr,
		const Gradient< Vector< T, 3 > > &	g0,
		const Vector< T, 3 > &	en
	)

Definition at line 64 of file GradientModel.h.

{
  SquareTensor<T,3> R, GT0;
  
  for(int i=0;i<3;i++)
    for(int j=0; j<3; j++)
    {
        if (i==j)
          R(i,j) = 1.0 - 2*en[i]*en[j];
        else
          R(i,j) = - 2*en[i]*en[j];

        GT0(i,j) = g0[j][i];
    }

  SquareTensor<T,3> GTR(R*GT0*R);
  
  for(int i=0;i<3;i++)
    for(int j=0; j<3; j++)
    {
        gr[j][i]  = GTR(i,j);
    }
}