@@ -13,7 +13,7 @@ See License.txt for details.
1313
1414// ----------------------------------------------------------------------------
1515
16- vtkStandardNewMacro (vtkPlusPolydataForce)
16+ vtkStandardNewMacro ( vtkPlusPolydataForce )
1717
1818// ----------------------------------------------------------------------------
1919vtkPlusPolydataForce::vtkPlusPolydataForce()
@@ -27,9 +27,9 @@ vtkPlusPolydataForce::vtkPlusPolydataForce()
2727}
2828
2929// ----------------------------------------------------------------------------
30- void vtkPlusPolydataForce::PrintSelf (ostream& os, vtkIndent indent)
30+ void vtkPlusPolydataForce::PrintSelf ( ostream& os, vtkIndent indent )
3131{
32- this ->Superclass ::PrintSelf (os, indent.GetNextIndent ());
32+ this ->Superclass ::PrintSelf ( os, indent.GetNextIndent () );
3333 os << indent.GetNextIndent () << " Gamma Sigmoid: " << this ->gammaSigmoid << endl;
3434}
3535
@@ -40,86 +40,86 @@ vtkPlusPolydataForce::~vtkPlusPolydataForce()
4040}
4141
4242// ----------------------------------------------------------------------------
43- void vtkPlusPolydataForce::SetInput (vtkPolyData * poly)
43+ void vtkPlusPolydataForce::SetInput ( vtkPolyData * poly )
4444{
4545 this ->poly = poly;
4646}
4747
4848// ----------------------------------------------------------------------------
49- int vtkPlusPolydataForce::GenerateForce (vtkMatrix4x4 * transformMatrix, double force[3 ])
49+ int vtkPlusPolydataForce::GenerateForce ( vtkMatrix4x4 * transformMatrix, double force[3 ] )
5050{
5151 double distance;
52- distance = CalculateDistance (transformMatrix->GetElement (0 , 3 ), transformMatrix->GetElement (1 , 3 ), transformMatrix->GetElement (2 , 3 ) );
52+ distance = CalculateDistance ( transformMatrix->GetElement ( 0 , 3 ), transformMatrix->GetElement ( 1 , 3 ), transformMatrix->GetElement ( 2 , 3 ) );
5353
54- if (distance <= 5 )
54+ if ( distance <= 5 )
5555 {
56- CalculateForce (transformMatrix->GetElement (0 , 3 ), transformMatrix->GetElement (1 , 3 ), transformMatrix->GetElement (2 , 3 ), force);
56+ CalculateForce ( transformMatrix->GetElement ( 0 , 3 ), transformMatrix->GetElement ( 1 , 3 ), transformMatrix->GetElement ( 2 , 3 ), force );
5757 }
5858 else
5959 {
60- force[0 ] = (0 );
61- force[1 ] = (0 );
62- force[2 ] = (0 );
60+ force[0 ] = ( 0 );
61+ force[1 ] = ( 0 );
62+ force[2 ] = ( 0 );
6363 }
6464 cout << " FORCE: " << force[0 ] << " , " << force[1 ] << " , " << force[2 ] << endl;
6565 return 1 ;
6666}
6767
6868// ----------------------------------------------------------------------------
69- double vtkPlusPolydataForce::CalculateDistance (double x, double y, double z)
69+ double vtkPlusPolydataForce::CalculateDistance ( double x, double y, double z )
7070{
7171 int pointID;
7272 double distance;
73- pointID = this ->poly ->FindPoint (x,y,z );
74- this ->poly ->GetPoint (pointID, this ->lastPos );
73+ pointID = this ->poly ->FindPoint ( x, y, z );
74+ this ->poly ->GetPoint ( pointID, this ->lastPos );
7575
76- double tmp = pow ((x- lastPos[0 ]), 2 ) + pow ((y- lastPos[1 ]), 2 ) + pow ((z- lastPos[2 ]), 2 );
77- distance = sqrt (tmp);
76+ double tmp = pow ( ( x - lastPos[0 ] ), 2 ) + pow ( ( y - lastPos[1 ] ), 2 ) + pow ( ( z - lastPos[2 ] ), 2 );
77+ distance = sqrt ( tmp );
7878 return distance;
7979}
8080
8181// ----------------------------------------------------------------------------
82- int vtkPlusPolydataForce::SetGamma (double gamma)
82+ int vtkPlusPolydataForce::SetGamma ( double gamma )
8383{
8484 gammaSigmoid = gamma;
8585
8686 return 0 ;
8787}
8888
8989// ----------------------------------------------------------------------------
90- void vtkPlusPolydataForce::CalculateForce (double x, double y, double z, double force[3 ])
90+ void vtkPlusPolydataForce::CalculateForce ( double x, double y, double z, double force[3 ] )
9191{
9292 double vector[3 ];
9393
94- vector[0 ] = fabs (x - this ->lastPos [0 ]);
95- vector[1 ] = fabs (y - this ->lastPos [1 ]);
96- vector[2 ] = fabs (z - this ->lastPos [2 ]);
94+ vector[0 ] = fabs ( x - this ->lastPos [0 ] );
95+ vector[1 ] = fabs ( y - this ->lastPos [1 ] );
96+ vector[2 ] = fabs ( z - this ->lastPos [2 ] );
9797
9898 cout << " vector[0]: " << vector[0 ] << endl;
9999 cout << " vector[1]: " << vector[1 ] << endl;
100100 cout << " vector[2]: " << vector[2 ] << endl;
101101
102- for (int i= 0 ; i<= 3 ; i++)
102+ for ( int i = 0 ; i < 3 ; i++ )
103103 {
104- if (vector[i] > 0 )
104+ if ( vector[i] > 0 )
105105 {
106- force[i] = (0.1 /( vector[i]* vector[i]))* .6 ;
106+ force[i] = ( 0.1 / ( vector[i] * vector[i] ) ) * .6 ;
107107 }
108108 }
109109 cout << " X: " << force[0 ] << endl;
110110 cout << " Y: " << force[1 ] << endl;
111111 cout << " Z: " << force[2 ] << endl;
112112
113- if (force[0 ] > 1 )
113+ if ( force[0 ] > 1 )
114114 {
115- force[0 ]= .6 ;
115+ force[0 ] = .6 ;
116116 }
117- if (force[1 ] > 1 )
117+ if ( force[1 ] > 1 )
118118 {
119- force[1 ]= .6 ;
119+ force[1 ] = .6 ;
120120 }
121- if (force[2 ] > 1 )
121+ if ( force[2 ] > 1 )
122122 {
123- force[2 ]= .6 ;
123+ force[2 ] = .6 ;
124124 }
125125}
0 commit comments