From 85e081cab1f60fcb91d3c6dcd4bf2ad3a0b97f36 Mon Sep 17 00:00:00 2001
From: Eduardo DAVILA <davila@localhost.localdomain>
Date: Wed, 8 Mar 2017 10:08:47 +0100
Subject: [PATCH] #3070 creaMaracasVisu Feature New Normal  - 
 bbmaracasvisuAxeVolume with fopenmp

---
 bbtk/src/bbmaracasvisuAxeVolume.cxx | 23 ++++++++++++++---------
 1 file changed, 14 insertions(+), 9 deletions(-)

diff --git a/bbtk/src/bbmaracasvisuAxeVolume.cxx b/bbtk/src/bbmaracasvisuAxeVolume.cxx
index 5cb228e..24d2ef0 100644
--- a/bbtk/src/bbmaracasvisuAxeVolume.cxx
+++ b/bbtk/src/bbmaracasvisuAxeVolume.cxx
@@ -93,27 +93,32 @@ void AxeVolume::Process()
 	mimage->SetOrigin(bbGetInputIn()->GetOrigin());
 	mimage->SetScalarType( outputformat );
 	mimage->AllocateScalars();
-	int i,j,k;
+
     int sizeLstPointR 	= bbGetInputlstPointR().size();
 	int iAxe,sizeAxe	= bbGetInputlstPointX().size();
-	double rx,ry,rz;
-	double r,rr;
-	unsigned short *p;
+	int ii;
 	int sizeImage 		= sizeX*sizeY*sizeZ;
-	double px,py,pz;
-	double px1,py1,pz1;
-	double px2,py2,pz2;
+	unsigned short *p;
+
 	// Clean image
 	p = (unsigned short*)mimage->GetScalarPointer (0, 0, 0);
-	for ( i=0 ; i<sizeImage ; i++)
+	for ( ii=0 ; ii<sizeImage ; ii++)
 	{
 		*p = 0;
 		p++;
 	}
 
+
+#pragma omp parallel for 
 	for (iAxe=0 ; iAxe<sizeAxe; iAxe++)
 	 {
-printf("AxeVolume %d/%d\n ",iAxe,sizeAxe);
+		int i,j,k;
+		double rx,ry,rz;
+		double r,rr;
+		double px,py,pz;
+		double px1,py1,pz1;
+		double px2,py2,pz2;
+		printf("AxeVolume %d/%d\n ",iAxe,sizeAxe);
 		if (iAxe<sizeLstPointR)
 		{
 			r = bbGetInputlstPointR()[ iAxe ]* invSpc[0];
-- 
2.45.0