Files
EgtGeomKernel/SbzStandard.cpp
Daniele Bariletti 5d62db3e73 EgtGeomKernel :
- aggiunta funzione per creare una sfera come sup Bezier.
2024-02-14 14:40:57 +01:00

81 lines
2.5 KiB
C++

//----------------------------------------------------------------------------
// EgalTech 2024
//----------------------------------------------------------------------------
// File : SbzSphere.cpp Data : 14.02.2024 Versione : 2.6b2
// Contenuto : Implementazione di funzioni per creazione di superfici Sbz
// standard : Box, Pyramid, Cylinder, Sphere, Cone.
//
//
// Modifiche : 14.02.2024 DB Creazione modulo.
//
//
//----------------------------------------------------------------------------
//--------------------------- Include ----------------------------------------
#include "stdafx.h"
#include "CurveArc.h"
#include "SurfTriMesh.h"
#include "SurfBezier.h"
#include "/EgtDev/Include/EGkSbzStandard.h"
using namespace std ;
//-------------------------------------------------------------------------------
ISurfBezier*
CreateBezierSphere( const Point3d& ptCenter, double dR)
{
// creo una superficie di Bezier di grado 2 con 45 punti di controllo
PtrOwner<ISurfBezier> pSrfBez( CreateSurfBezier()) ;
int nDegU = 2 ;
int nDegV = 2 ;
int nSpanU = 4 ; // i poli della sfera sono a coordinate ( 0,0,-R) e ( 0,0,R)
int nSpanV = 2 ;
bool bRat = true ;
double dW = SQRT2 ;
pSrfBez->Init(nDegU, nDegV, nSpanU, nSpanV, bRat) ;
// polo inferiore // dW = 1, dW = SQRT2 / 2
for ( int i = 0 ; i < 9; ++i) {
if ( i % 2 == 0)
dW = 1 ;
else
dW = SQRT2 / 2 ;
pSrfBez->SetControlPoint( i, 0, ptCenter + Point3d( 0, 0, -dR), dW) ;
}
// definisco la gabbia esterna // dW = SQRT2 / 2, dW = 1 / 2
// parto dal punto ( 0,-dR, -dR) e completo riga per riga
double dH = -dR ;
for ( int j = 1 ; j < 4 ; ++j ) {
Vector3d vtDir ( 0, -dR, 0) ;
Point3d pt(-dR,-dR,dH) ;
for ( int i = 0; i < 9 ; ++i ) {
// ogni due punti ruoto di 90 gradi a sinistra
if ( i%2 == 0) {
vtDir.Rotate( Z_AX, 90) ;
if ( j%2 == 1)
dW = SQRT2 / 2 ;
else
dW = 1 ;
}
else {
if ( j%2 == 1)
dW = 1. / 2. ;
else
dW = SQRT2 / 2 ;
}
pt += vtDir ;
pSrfBez->SetControlPoint( i, j, ptCenter + pt, dW) ;
}
dH += dR ;
}
// polo superiore // dW = 1, dW = SQRT2 / 2
for ( int i = 0 ; i < 9; ++i) {
if ( i % 2 == 0)
dW = 1 ;
else
dW = SQRT2 / 2 ;
pSrfBez->SetControlPoint( i, 4, ptCenter + Point3d( 0, 0, dR), dW) ;
}
return Release( pSrfBez) ;
}