StabilizationUtils.h

Go to the documentation of this file.

//==============================================================================
//==============================================================================
 
#ifndef STABILIZATIONUTILS_H_
#define STABILIZATIONUTILS_H_
 
#include <vector>
#include "Vec3.h"
#include "MatVec.h"
 
namespace StabilizationUtils {
  double getElementSize(const std::vector<Vec3>& XC, int nsd);
 
  double getTauPt(double dt, double mu, const Vector& U, const Matrix& G,
                  const double Ct = 2.0, const double Cl = 36.0);
 
  std::pair<double,double>
  getTauNSPt(double dt, double mu, const Vector& U, const Matrix& G,
             const double Ct = 2.0, const double Cl = 36.0);
 
  std::pair<double,double>
  getTauNSALEPt(double dt, double mu, const Vector& U, const Matrix& G,
                const double Ct = 2.0, const double Cl = 36.0);
 
  Vector getTauPtJac(const Vector& U, const Matrix& G,
                   const double tauM);
 
  std::pair<Vector,Vector>
  getTauNSPtJac(const Vector& U, const Matrix& G,
                const double tauM, const double& tauC);
}
 
#endif