Face.h

Go to the documentation of this file.

//
//  This file is part of the FFEA simulation package
//
//  Copyright (c) by the Theory and Development FFEA teams,
//  as they appear in the README.md file.
//
//  FFEA is free software: you can redistribute it and/or modify
//  it under the terms of the GNU General Public License as published by
//  the Free Software Foundation, either version 3 of the License, or
//  (at your option) any later version.
//
//  FFEA is distributed in the hope that it will be useful,
//  but WITHOUT ANY WARRANTY; without even the implied warranty of
//  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
//  GNU General Public License for more details.
//
//  You should have received a copy of the GNU General Public License
//  along with FFEA.  If not, see <http://www.gnu.org/licenses/>.
//
//  To help us fund FFEA development, we humbly ask that you cite
//  the research papers on the package.
//

#ifndef FACE_H_INCLUDED
#define FACE_H_INCLUDED

#include <math.h>

#include "mesh_node.h"
#include "tetra_element_linear.h"
#include "SecondOrderFunctions.h"
#include "SimulationParams.h"
#include "CheckTetrahedraOverlap.h"
#include "VolumeIntersection.h"
#include "mat_vec_types.h"
#include "mat_vec_fns.h"

class Face {
public:
    Face();

    ~Face();

    mesh_node *n[4];

    int index;

    tetra_element_linear *e;

    int ssint_interaction_type;

    scalar area_0;

    scalar area;

    vector3 normal;

    vector3 centroid;

    vector3 force[4];


    SecondOrderFunctions::stu centroid_stu;

    bool ssint_xz_interaction_flag;
    // bool *ssint_bb_interaction_flag; // DEPRECATED
    bool kinetically_active;


    int num_blobs;
    // vector3 *vdw_bb_force; // DEPRECATED
    // scalar *vdw_bb_energy; // DEPRECATED
    // vector3 *vdw_xz_force; // DEPRECATED
    // scalar vdw_xz_energy; // DEPRECATED

    bool checkTetraIntersection(Face *f2);
    bool checkTetraIntersection(Face *f2, scalar *blob_corr,int f1_daddy_blob_index,int f2_daddy_blob_index);

    scalar getTetraIntersectionVolume(Face *f2);
    scalar getTetraIntersectionVolume(Face *f2, scalar *blob_corr,int f1_daddy_blob_index,int f2_daddy_blob_index);

    scalar checkTetraIntersectionAndGetVolume(Face *f2);
    scalar checkTetraIntersectionAndGetVolume(Face *f2,scalar *blob_corr,int f1_daddy_blob_index,int f2_daddy_blob_index);


    void getTetraIntersectionVolumeAndArea(Face *f2, geoscalar &vol, geoscalar &area);
    void getTetraIntersectionVolumeAndArea(Face *f2, geoscalar &vol, geoscalar &area, scalar *blob_corr,int f1_daddy_blob_index,int f2_daddy_blob_index);


    bool getTetraIntersectionVolumeTotalGradientAndShapeFunctions(Face *f2, geoscalar dr, grr3 (&dVdr), geoscalar &vol, grr4 (&phi1), grr4 (&phi2));
    bool getTetraIntersectionVolumeTotalGradientAndShapeFunctions(Face *f2, geoscalar dr, grr3 (&dVdr), geoscalar &vol, grr4 (&phi1), grr4 (&phi2), scalar *blob_corr, int f1_daddy_blob_index, int f2_daddy_blob_index);


    Blob *daddy_blob;

    int init(int index, tetra_element_linear *e, mesh_node *n0, mesh_node *n1, mesh_node *n2, mesh_node *oposite, SecondOrderFunctions::stu centroid_stu, Blob *daddy_blob, SimulationParams *params);

    int init(int index, mesh_node *n0, mesh_node *n1, mesh_node *n2, mesh_node *opposite, Blob *daddy_blob, SimulationParams *params);

    void set_ssint_interaction_type(int ssint_interaction_type);

    int build_opposite_node();

    void calc_area_normal_centroid();
    void set_kinetic_state(bool state);

    vector3 * get_centroid();
    void print_centroid();
    void print_nodes();
    scalar get_area();

    void barycentric_calc_point(scalar b1, scalar b2, scalar b3, vector3 *p);
    void barycentric_calc_point_f2(scalar b1, scalar b2, scalar b3, vector3 *p,scalar *blob_corr,int f1_daddy_blob_index,int f2_daddy_blob_index);

    scalar average_phi();

    scalar get_normal_flux();

    void add_force_to_node(int i, vector3 *f);
    template <class brr3> void add_force_to_node(int i, brr3 (&f));

    void add_force_to_node_atomic(int i, vector3 *f);

    // void add_bb_vdw_force_to_record(vector3 *f, int other_blob_index); // DEPRECATED
    // template <class brr3> void add_bb_vdw_force_to_record(brr3 &f, int other_blob_index); // DEPRECATED

    // void add_bb_vdw_energy_to_record(scalar energy, int other_blob_index); // DEPRECATED

    // void add_xz_vdw_force_to_record(vector3 *f); // DEPRECATED

    // void add_xz_vdw_energy_to_record(scalar energy); // DEPRECATED

    void zero_force();

    // void zero_vdw_bb_measurement_data(); // DEPRECATED

    // void zero_vdw_xz_measurement_data(); // DEPRECATED

    void set_ssint_xz_interaction_flag(bool state);

    // void set_vdw_bb_interaction_flag(bool state, int other_blob_index); // DEPRECATED

    template <class brr3> void vec3Vec3SubsToArr3Mod(Face *f2, brr3 (&w), scalar *blob_corr,int f1_daddy_blob_index,int f2_daddy_blob_index);

    bool is_ssint_active();
    bool is_kinetic_active();

    scalar length_of_longest_edge(); 

private:
    int stuff;
    bool dealloc_n3;

};

#endif