From 7a603f9bbbd0744155d5d287e733902c5c714c52 Mon Sep 17 00:00:00 2001
From: Jean-Francois Remacle <jean-francois.remacle@uclouvain.be>
Date: Wed, 15 Nov 2006 20:52:57 +0000
Subject: [PATCH] *** empty log message ***
---
Geo/OCCVertex.cpp | 58 +++++++++++++++++++++++++++++++++++++++++++++++
1 file changed, 58 insertions(+)
create mode 100644 Geo/OCCVertex.cpp
diff --git a/Geo/OCCVertex.cpp b/Geo/OCCVertex.cpp
new file mode 100644
index 0000000000..b0137cc6fe
--- /dev/null
+++ b/Geo/OCCVertex.cpp
@@ -0,0 +1,58 @@
+// $Id: OCCVertex.cpp,v 1.1 2006-11-15 20:52:57 remacle Exp $
+//
+// Copyright (C) 1997-2006 C. Geuzaine, J.-F. Remacle
+//
+// This program 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 2 of the License, or
+// (at your option) any later version.
+//
+// This program 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 this program; if not, write to the Free Software
+// Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307
+// USA.
+//
+// Please report all bugs and problems to <gmsh@geuz.org>.
+
+#if defined(HAVE_OCC)
+#include "GModel.h"
+#include "OCCVertex.h"
+#include "OCCEdge.h"
+#include "OCCFace.h"
+
+double max_surf_curvature ( double x, double y, double z , const GEdge *_myGEdge)
+{
+ std::list<GFace *> faces = _myGEdge->faces();
+ std::list<GFace *>::iterator it = faces.begin();
+ double curv = 0;
+ while (it != faces.end())
+ {
+ SPoint2 par = (*it)->parFromPoint(SPoint3 (x,y,z));
+ double cc = (*it)->curvature ( par );
+ if (cc < 1.e4) curv = std::max(curv, cc );
+ ++it;
+ }
+ return curv;
+}
+
+
+double OCCVertex::max_curvature_of_surfaces() const
+{
+ if (max_curvature <0)
+ {
+ for (std::list<GEdge*> :: const_iterator it = l_edges.begin() ; it != l_edges.end() ; ++it )
+ {
+ max_curvature = std::max ( max_surf_curvature (x(), y(), z(), *it) , max_curvature);
+ }
+ // printf("max curvature (%d) = %12.5E lc = %12.5E\n",tag(),max_curvature,prescribedMeshSizeAtVertex());
+
+ }
+ return max_curvature;
+}
+
+#endif
--
GitLab