/export/home/gkanscha/ANS/submissions/Salgado/include/Material.h

Go to the documentation of this file.

#ifndef _MATERIAL_PARAMETERS_H_
#define _MATERIAL_PARAMETERS_H_

#include <cmath>
#include <utility>

enum TestCase{
  TEST_CASE_MOVE, 
  TEST_CASE_SPLIT, 
  TEST_CASE_MERGE 
};

struct material_function{
  public:
    material_function( const double _m1 = 1., const double _m2 = 1., const double _delta = 1. );
    material_function( const std::pair<double, double> &params, const double _delta );
    double deriv( const double x ) const;
    double max_deriv() const;
    double operator()( const double x ) const;
    double average( const double x, const double y ) const;
  private:
    const double avg, 
                 diff, 
                 delta; 

    double value ( const double x ) const;
};

struct GLPotential{
  public:
    GLPotential();
    double operator()( const double x ) const;
  private:
    const double max_second_deriv;
  public:
    const double A;
};

struct InterfaceEnergy{
  public:
    InterfaceEnergy( const double gamma, const double theta_s );
    double operator()( const double x ) const;
  private:
    const double coeff, 
                 max_second_deriv; 
  public:
    const double B;
};

struct Material_Parameters{
  public:
    Material_Parameters( const double epsilon1, const double epsilon2, const double epsilon3,
                         const double K1, const double K2, const double M1, const double M2,
                         const double rho1, const double rho2, const double beta1, const double beta2,
                         const double eta1, const double eta2, const double _lambda, const double _gamma,
                         const double theta_s, const double _delta, const double _alpha );
    Material_Parameters( std::pair<double, double> _epsilon, const double epsilon3,
                         std::pair<double, double> _K, 
                         std::pair<double, double> _M, std::pair<double, double> _rho,
                         std::pair<double, double> _beta, 
                         std::pair<double, double> _eta, const double _lambda, const double _gamma,
                         const double theta_s, const double _delta, const double _alpha );
    material_function epsilon, K, M, rho, beta, eta;
    const double lambda, gamma, delta, alpha, epsilon_plate;
    GLPotential W;
    InterfaceEnergy gamma_fs;
};

#endif //_MATERIAL_PARAMETERS_H_