smesh/src/StdMeshers/StdMeshers_Quadrangle_2D.hxx

//  SMESH SMESH : implementaion of SMESH idl descriptions
//
//  Copyright (C) 2003  OPEN CASCADE, EADS/CCR, LIP6, CEA/DEN,
//  CEDRAT, EDF R&D, LEG, PRINCIPIA R&D, BUREAU VERITAS 
// 
//  This library is free software; you can redistribute it and/or 
//  modify it under the terms of the GNU Lesser General Public 
//  License as published by the Free Software Foundation; either 
//  version 2.1 of the License. 
// 
//  This library 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 
//  Lesser General Public License for more details. 
// 
//  You should have received a copy of the GNU Lesser General Public 
//  License along with this library; if not, write to the Free Software 
//  Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307 USA 
// 
// See http://www.salome-platform.org/ or email : webmaster.salome@opencascade.com
//
//
//
//  File   : StdMeshers_Quadrangle_2D.hxx
//           Moved here from SMESH_Quadrangle_2D.hxx
//  Author : Paul RASCLE, EDF
//  Module : SMESH
//  $Header$

#ifndef _SMESH_QUADRANGLE_2D_HXX_
#define _SMESH_QUADRANGLE_2D_HXX_

#include "SMESH_2D_Algo.hxx"
#include "SMESH_Mesh.hxx"
#include "Utils_SALOME_Exception.hxx"

#include "gp_XY.hxx"

#include "SMESH_MesherHelper.hxx"

//class SMDS_MeshNode;

typedef struct uvPtStruct
{
  double param;
  double normParam;
  double u; // original 2d parameter
  double v;
  double x; // 2d parameter, normalized [0,1]
  double y; 
  const SMDS_MeshNode * node;
} UVPtStruct;

typedef struct faceQuadStruct
{
  int nbPts[4];
  TopoDS_Edge edge[4];
  double first[4];
  double last[4];
  bool isEdgeForward[4];
  bool isEdgeOut[4]; // true, if an edge has more nodes, than the opposite
  UVPtStruct* uv_edges[4];
  UVPtStruct* uv_grid;
} FaceQuadStruct;

class StdMeshers_Quadrangle_2D:
  public SMESH_2D_Algo
{
public:
  StdMeshers_Quadrangle_2D(int hypId, int studyId, SMESH_Gen* gen);
  virtual ~StdMeshers_Quadrangle_2D();

  virtual bool CheckHypothesis(SMESH_Mesh& aMesh,
                               const TopoDS_Shape& aShape,
                               SMESH_Hypothesis::Hypothesis_Status& aStatus);

  virtual bool Compute(SMESH_Mesh& aMesh,
		       const TopoDS_Shape& aShape)
    throw (SALOME_Exception);

  FaceQuadStruct* CheckAnd2Dcompute(SMESH_Mesh& aMesh,
				    const TopoDS_Shape& aShape,
                                    const bool CreateQuadratic)
    throw (SALOME_Exception);

  static void QuadDelete(FaceQuadStruct* quad);

  /**
   * Returns NLinkNodeMap from myTool
   */
  const NLinkNodeMap& GetNLinkNodeMap() { return myTool->GetNLinkNodeMap(); }

  ostream & SaveTo(ostream & save);
  istream & LoadFrom(istream & load);
  friend ostream & operator << (ostream & save, StdMeshers_Quadrangle_2D & hyp);
  friend istream & operator >> (istream & load, StdMeshers_Quadrangle_2D & hyp);

protected:

  FaceQuadStruct* CheckNbEdges(SMESH_Mesh& aMesh,
                               const TopoDS_Shape& aShape)
    throw (SALOME_Exception);

  void SetNormalizedGrid(SMESH_Mesh& aMesh,
			 const TopoDS_Shape& aShape,
			 FaceQuadStruct*& quad)
    throw (SALOME_Exception);

  /**
   * Special function for creation only quandrangle faces
   */
  bool ComputeQuadPref(SMESH_Mesh& aMesh,
                       const TopoDS_Shape& aShape,
                       FaceQuadStruct* quad)
    throw (SALOME_Exception);

  UVPtStruct* LoadEdgePoints2(SMESH_Mesh& aMesh,
			      const TopoDS_Face& F, const TopoDS_Edge& E,
                              bool IsReverse);

  UVPtStruct* LoadEdgePoints(SMESH_Mesh& aMesh,
			     const TopoDS_Face& F, const TopoDS_Edge& E,
			     double first, double last);

  UVPtStruct* MakeEdgePoints(SMESH_Mesh& aMesh,
			     const TopoDS_Face& F, const TopoDS_Edge& E,
			     double first, double last, int nb_segm);

  // true if QuadranglePreference hypothesis is assigned that forces
  // construction of quadrangles if the number of nodes on opposite edges
  // is not the same in the case where the global number of nodes on edges is even
  bool myQuadranglePreference;

  SMESH_MesherHelper* myTool; // toll for working with quadratic elements
};

#endif