From 3ac2499bedf736e9c70eda8e1d68619e9937f1d0 Mon Sep 17 00:00:00 2001
From: Jean-Francois Remacle <jean-francois.remacle@uclouvain.be>
Date: Fri, 2 Mar 2007 09:33:36 +0000
Subject: [PATCH] *** empty log message ***
---
Geo/gmshSurface.cpp | 9 +++++----
1 file changed, 5 insertions(+), 4 deletions(-)
diff --git a/Geo/gmshSurface.cpp b/Geo/gmshSurface.cpp
index ec0a1c5347..258bc9e8a9 100644
--- a/Geo/gmshSurface.cpp
+++ b/Geo/gmshSurface.cpp
@@ -1,4 +1,4 @@
-// $Id: gmshSurface.cpp,v 1.5 2007-03-02 09:20:21 remacle Exp $
+// $Id: gmshSurface.cpp,v 1.6 2007-03-02 09:33:36 remacle Exp $
//
// Copyright (C) 1997-2007 C. Geuzaine, J.-F. Remacle
//
@@ -73,7 +73,7 @@ gmshSurface * gmshPolarSphere::NewPolarSphere(int iSphere, double x, double y, d
SPoint3 gmshPolarSphere::point (double parA, double parB) const
{
- double ra,phi;
+ /*double ra,phi;
ra=hypot(parA,parB);
phi=2*atan((parB/ra)/(1+parA/ra));
double par1=-phi;
@@ -83,8 +83,9 @@ SPoint3 gmshPolarSphere::point (double parA, double parB) const
const double x = xc + r * sin(par2) * cos(par1);
const double y = yc + r * sin(par2) * sin(par1);
const double z = zc - r * cos(par2);
- // printf("%g %g - %g %g %g\n",par1,par2,x,y,z);
- return SPoint3(x, y, z);
+ // printf("%g %g - %g %g %g\n",par1,par2,x,y,z);*/
+ double f=2*r/(parA*parA+parB*parB+4*r*r);
+ return SPoint3(f*2*parA*r, f*2*parB*r, f*(parA*parA+parB*parB));
}
--
GitLab