From 7a21ac205e7e00b99caa09ecde5dd55c1b6bb711 Mon Sep 17 00:00:00 2001
From: Julien Loiseau <jloiseau@lanl.gov>
Date: Tue, 20 Jun 2017 16:15:08 -0600
Subject: [PATCH] Fixed sodtube

---
 app/sodtube/generator/main.cc | 25 ++++---------------------
 1 file changed, 4 insertions(+), 21 deletions(-)

diff --git a/app/sodtube/generator/main.cc b/app/sodtube/generator/main.cc
index 5b394514..81eddda6 100644
--- a/app/sodtube/generator/main.cc
+++ b/app/sodtube/generator/main.cc
@@ -99,14 +99,16 @@ int main(int argc, char * argv[]){
       P[part] = pressure_2;
       rho[part] = rho_2; 
       u[part] = u_2;
-      m[part] = m_2;
+      //m[part] = m_2;
     }else{
       P[part] = pressure_1;
       rho[part] = rho_1;
       u[part] = u_1;
-      m[part] = m_1;
+      //m[part] = m_1;
     }
 
+    m[part] = rho[part]*middle/(nparticles/2.);
+
     //m[part] = 0.;
     // Y and Z not used 
     // VX, VY, VZ and AX, AY, AZ stay to 0
@@ -119,25 +121,6 @@ int main(int argc, char * argv[]){
     //std::cout<<x[part]<<": "<<h[part]<<std::endl;
   }
 
-  // For each particle, search the neighbors to find the mass
-  //for(int64_t p1=0;p1<nparticlesproc; ++p1){
-  //  for(int64_t p2=0;p2<nparticlesproc;++p2){
-  //    if( fabs(x[p1]-x[p2]) < 2*smoothing_length ){
-         // Then consider it for my mass 
-  //      m[p1] += rho[p2]/kernel::cubic_spline_kernel
-  //        (fabs(x[p1]-x[p2]),smoothing_length);
-  //    }
-  //  }
-  //}
-
-  //double totalmass = 0.;
-  //for(int64_t p1=0;p1<nparticlesproc; ++p1){
-  //  totalmass += m[p1];  
-  //}
-  //for(int64_t p1=0;p1<nparticlesproc; ++p1){
-  //  m[p1] = m[p1]/totalmass;
-  //}
-
   char filename[128];
   //sprintf(filename,"%s_%d.h5part",fileprefix,nparticles);
   sprintf(filename,"%s.h5part",fileprefix,nparticles);