pzhyperplane.h

Go to the documentation of this file.

/*****************************************************************************
 * O contedo desse arquivo �de propriedade do LabMeC-DES-FEC-UNICAMP e do
 * CENPES-Petrobras. 
 * O uso de qualquer parte ou do todo est�condicionado �expressa autoriza�o
 * dos propriet�ios.
 *****************************************************************************/

#ifndef PZREADMESHHR_H
#define PZREADMESHHR_H

#include "pzmanvector.h"
#include <string>

class TPZGeoMesh;
class TPZGeoEl;
class TPZCompMesh;
class TPZGeoElSide;

class TPZHyperPlane
{
    
    TPZManVector<REAL,3> fCenter;
    
    TPZManVector<REAL,3> fNormal;
    
public:

    TPZHyperPlane(TPZVec<REAL> &center, TPZVec<REAL> &normal) : fCenter(center), fNormal(normal)
    {
        
    }

    bool IsLeft(TPZVec<REAL> &point) const
    {
        REAL val=0.;
        for (int i=0; i<3; i++) {
            val += (point[i]-fCenter[i])*fNormal[i];
        }
        return val <= 0;
    }
    
    REAL Distance(const TPZVec<REAL> &point, TPZVec<REAL> &jacobian) const
    {
        REAL distance = 0;
        for (int i=0; i<3; i++) {
            distance += (point[i]-fCenter[i])*fNormal[i];
            jacobian[i] = fNormal[i];
        }
        return distance;
    }
};

class TPZHyperPlaneIntersect
{
public:
    void Intersect(TPZGeoMesh &input, const TPZHyperPlane &plane, TPZGeoMesh &target);
    
    REAL EdgeIntersect(const TPZGeoElSide &gelside, const TPZHyperPlane &plane);
    
private:
    void Reorder(TPZGeoEl *gel, TPZGeoMesh &target, TPZVec<std::pair<int64_t,int64_t> > &sidenodepair);
    int ReorderGeneral(TPZGeoMesh &target, TPZVec<std::pair<int64_t,int64_t> > &sidenodepair);
};

#endif