#include <weight_parameter_function.h>

Public Member Functions
	weight_parameter_function (libMesh::Mesh &alpha_mesh, double alpha_eps, double alpha_coupling_BIG, int subdomain_idx_BIG, int subdomain_idx_micro, int subdomain_idx_both, int subdomain_idx_coupling)

	weight_parameter_function (libMesh::Mesh &alpha_mesh)

	~weight_parameter_function ()

void	set_parameters (double alpha_eps, double alpha_coupling_BIG, int subdomain_idx_BIG, int subdomain_idx_micro, int subdomain_idx_both, int subdomain_idx_coupling)

void	set_parameters (std::string &filename)

double	get_alpha (const libMesh::Point &qpoint, WeightFunctionSystemType system_type)

double	get_alpha_BIG (const libMesh::Point &qpoint)

double	get_alpha_micro (const libMesh::Point &qpoint)

void	clear ()

Protected Attributes
libMesh::Mesh &	m_alpha_mesh

std::unique_ptr< libMesh::PointLocatorBase >	m_locator_unique_ptr

double	m_alpha_eps

double	m_alpha_coupling_BIG

double	m_alpha_coupling_micro

int	m_subdomain_idx_BIG

int	m_subdomain_idx_micro

int	m_subdomain_idx_both

int	m_subdomain_idx_coupling

Detailed Description

Definition at line 13 of file weight_parameter_function.h.

Constructor & Destructor Documentation

weight_parameter_function::weight_parameter_function	(	libMesh::Mesh &	alpha_mesh,
		double	alpha_eps,
		double	alpha_coupling_BIG,
		int	subdomain_idx_BIG,
		int	subdomain_idx_micro,
		int	subdomain_idx_both,
		int	subdomain_idx_coupling
	)

inline

Definition at line 33 of file weight_parameter_function.h.

                                                                                                                             :
         m_alpha_mesh{ alpha_mesh },
         m_locator_unique_ptr { alpha_mesh.sub_point_locator() },
         m_alpha_eps { alpha_eps },
         m_alpha_coupling_BIG { alpha_coupling_BIG },
         m_alpha_coupling_micro { 1. - alpha_coupling_BIG },
         m_subdomain_idx_BIG { subdomain_idx_BIG },
         m_subdomain_idx_micro { subdomain_idx_micro },
         m_subdomain_idx_both { subdomain_idx_both },
         m_subdomain_idx_coupling { subdomain_idx_coupling }
     {
     };

weight_parameter_function::weight_parameter_function ( libMesh::Mesh & alpha_mesh )

inline

Definition at line 48 of file weight_parameter_function.h.

                                                           :
         m_alpha_mesh{ alpha_mesh },
         m_locator_unique_ptr { alpha_mesh.sub_point_locator() },
         m_alpha_eps { 10E-2 },
         m_alpha_coupling_BIG { 0.5 },
         m_alpha_coupling_micro { 0.5 },
         m_subdomain_idx_BIG { -1 },
         m_subdomain_idx_micro { -1 },
         m_subdomain_idx_both { -1 },
         m_subdomain_idx_coupling { -1 }
     {
     };

weight_parameter_function::~weight_parameter_function ( )

inline

Definition at line 62 of file weight_parameter_function.h.

     {
         clear();
     }

Member Function Documentation

void weight_parameter_function::clear ( )

inline

Definition at line 207 of file weight_parameter_function.h.

     {
         m_locator_unique_ptr.reset(NULL);
     };

double weight_parameter_function::get_alpha	(	const libMesh::Point &	qpoint,
		WeightFunctionSystemType	system_type
	)

inline

Definition at line 137 of file weight_parameter_function.h.

     {
         switch (system_type)
         {
             case WeightFunctionSystemType::MICRO :  
                     return this->get_alpha_micro(qpoint);
                     break;
             case WeightFunctionSystemType::MACRO :  
                     return this->get_alpha_BIG(qpoint);
                     break;
             case WeightFunctionSystemType::NO_WEIGHT :  
                     return 1.;
                     break;
         }
         homemade_error_msg("Why are you here!?");
     }

double weight_parameter_function::get_alpha_BIG ( const libMesh::Point & qpoint )

inline

Definition at line 154 of file weight_parameter_function.h.

     {
         libMesh::PointLocatorBase& locator = *m_locator_unique_ptr.get();
         const libMesh::Elem* elem = locator(qpoint);
         double output = 0;
 
         if(elem->subdomain_id() == m_subdomain_idx_BIG)
         {
             output = 1;
         }
         else if(elem->subdomain_id() == m_subdomain_idx_micro)
         {
             output = 0;
         }
         else if(elem->subdomain_id() == m_subdomain_idx_both)
         {
             output = m_alpha_eps;
         }
         else if(elem->subdomain_id() == m_subdomain_idx_coupling)
         {
             output = m_alpha_coupling_BIG;
         }
 
         return output;
     }

double weight_parameter_function::get_alpha_micro ( const libMesh::Point & qpoint )

inline

Definition at line 180 of file weight_parameter_function.h.

     {
         libMesh::PointLocatorBase& locator = *m_locator_unique_ptr.get();
         const libMesh::Elem* elem = locator(qpoint);
         double output = 0;
 
         if(elem->subdomain_id() == m_subdomain_idx_BIG)
         {
             output = 0;
         }
         else if(elem->subdomain_id() == m_subdomain_idx_micro)
         {
             output = 1;
         }
         else if(elem->subdomain_id() == m_subdomain_idx_both)
         {
             output = 1 - m_alpha_eps;
         }
         else if(elem->subdomain_id() == m_subdomain_idx_coupling)
         {
             output = m_alpha_coupling_micro;
         }
 
         return output;
     }

void weight_parameter_function::set_parameters	(	double	alpha_eps,
		double	alpha_coupling_BIG,
		int	subdomain_idx_BIG,
		int	subdomain_idx_micro,
		int	subdomain_idx_both,
		int	subdomain_idx_coupling
	)

inline

Definition at line 67 of file weight_parameter_function.h.

     {
         m_alpha_eps = alpha_eps;
 
         m_alpha_coupling_BIG = alpha_coupling_BIG;
         m_alpha_coupling_micro = 1. - alpha_coupling_BIG;
 
         m_subdomain_idx_BIG  = subdomain_idx_BIG;
         m_subdomain_idx_micro = subdomain_idx_micro;
         m_subdomain_idx_both = subdomain_idx_both;
         m_subdomain_idx_coupling = subdomain_idx_coupling;
     };

void weight_parameter_function::set_parameters ( std::string & filename )

inline

Definition at line 81 of file weight_parameter_function.h.

     {
         GetPot field_parser;
         field_parser.parse_input_file(filename, "#", "\n", " \t\n");
 
         if (field_parser.search(1, "MacroDomainOnlyIdx")) {
             m_subdomain_idx_BIG = field_parser.next(
                     m_subdomain_idx_BIG);
         } else {
             homemade_error_msg("Missing the macro domain only index!");
         }
 
         if (field_parser.search(1, "MicroDomainOnlyIdx")) {
             m_subdomain_idx_micro = field_parser.next(
                     m_subdomain_idx_micro);
         } else {
             homemade_error_msg("Missing the micro domain only index!");
         }
 
         if (field_parser.search(1, "CouplingDomainIdx")) {
             m_subdomain_idx_coupling = field_parser.next(
                     m_subdomain_idx_coupling);
         } else {
             homemade_error_msg("Missing the coupling domain index!");
         }
 
         if (field_parser.search(1, "BothDomainsIdx")) {
             m_subdomain_idx_both = field_parser.next(
                     m_subdomain_idx_both);
         } else {
             m_subdomain_idx_both = -1;
             std::cout << " >> Warning: Both system's domain index not defined!" << std::endl;
         }
 
         if (field_parser.search(1, "AlphaCouplingMacro")) {
             m_alpha_coupling_BIG = field_parser.next(
                     m_alpha_coupling_BIG);
             m_alpha_coupling_micro = 1. - m_alpha_coupling_BIG;
         } else if (field_parser.search(1, "AlphaCouplingMicro")) {
             m_alpha_coupling_micro = field_parser.next(
                     m_alpha_coupling_micro);
             m_alpha_coupling_BIG = 1. - m_alpha_coupling_micro;
         } else {
             m_alpha_coupling_BIG = 0.5;
             m_alpha_coupling_micro = 0.5;
             std::cout << " >> Warning: Coupling weights not defined, using 0.5 for both systems!" << std::endl;
         }
 
         if (field_parser.search(1, "AlphaEps")) {
             m_alpha_eps = field_parser.next(m_alpha_eps);
         } else {
             m_alpha_eps = 1e-2;
             std::cout << " >> Warning: Alpha epsilon not defined, using 1e-2!" << std::endl;
         }
     }