+ mimage->Delete();
+ }
+
+ int ext[6];
+ bbGetInputIn()->GetExtent(ext);
+ int sizeX=ext[1]-ext[0]+1;
+ int sizeY=ext[3]-ext[2]+1;
+ int sizeZ=ext[5]-ext[4]+1;
+
+ double spc[3];
+ bbGetInputIn()->GetSpacing(spc);
+ double invSpc[3];
+ invSpc[0] = 1/spc[0];
+ invSpc[1] = 1/spc[1];
+ invSpc[2] = 1/spc[2];
+
+ int outputformat = GetTypeFormat( bbGetInputOutputFormat() , bbGetInputIn() );
+
+ mimage = vtkImageData::New();
+ mimage->SetSpacing(bbGetInputIn()->GetSpacing());
+ mimage->SetDimensions(bbGetInputIn()->GetDimensions());
+ mimage->SetExtent(bbGetInputIn()->GetExtent());
+ mimage->SetOrigin(bbGetInputIn()->GetOrigin());
+ //EED 2017-01-01 Migration VTK7
+ #if VTK_MAJOR_VERSION <= 5
+ mimage->SetScalarType( outputformat );
+ mimage->AllocateScalars();
+ #else
+ mimage->AllocateScalars( outputformat , 1);
+ #endif
+
+ int sizeLstPointR = bbGetInputlstPointR().size();
+ int iAxe,sizeAxe = bbGetInputlstPointX().size();
+ int ii;
+ int sizeImage = sizeX*sizeY*sizeZ;
+ // unsigned short *p;
+
+ // Clean image
+ memset( (char*)mimage->GetScalarPointer(0,0,0) , 0, sizeImage*(mimage->GetScalarSize()) );
+
+ // p = (unsigned short*)mimage->GetScalarPointer (0, 0, 0);
+ // for ( ii=0 ; ii<sizeImage ; ii++)
+ // {
+ // *p = 0;
+ // p++;
+ // }
+
+
+ DEF_POINTER_IMAGE_VTK_CREA(vI,ssI,pI,stI,mimage)
+
+
+ int k1omp=(double)(sizeAxe-1)*0.0;
+ int k2omp=(double)(sizeAxe-1)*0.2;
+ int k3omp=(double)(sizeAxe-1)*0.4;
+ int k4omp=(double)(sizeAxe-1)*0.6;
+ int k5omp=(double)(sizeAxe-1)*0.8;
+ int k6omp=(double)(sizeAxe-1)*1.0;
+
+ printf("EED AxeVolume::Process (with openmp)\n");
+
+ #pragma omp parallel for
+ for (iAxe=0 ; iAxe<sizeAxe; iAxe++)
+ {
+ if ( (k1omp==iAxe) || (k2omp==iAxe) || (k3omp==iAxe) ||
+ (k4omp==iAxe) || (k5omp==iAxe) || (k6omp==iAxe) ) { printf(" %d%\n", (int)(((double)iAxe/(double)(sizeAxe-1))*100 )); }
+ int i,j,k;
+ double rx,ry,rz;
+ double r,rr;
+ double px,py,pz;
+ double px1,py1,pz1;
+ double px2,py2,pz2;
+ double vItmpOMP;
+
+ if (iAxe<sizeLstPointR)