working on kickdrift #94, fixed update pos, still pendant stresses an…

…d strain

src/cuda/Mechanical.cu

@@ -166,6 +166,7 @@ __device__ inline void Domain_d::UpdatePos(double dt){
 
 	if ( i < particle_count ) {
     	x[i] += dt*v[i];   
+      u[i] += dt*v[i];
   }
 }
 

src/main_mech.cu

@@ -262,7 +262,11 @@ int main(int argc, char **argv) //try
 	dom_d->auto_ts = true;
   dom_d->Alpha = 1.0;
 	//dom_d->MechSolve(0.0101,1.0e-4);
-  dom_d->MechKickDriftSolve(0.0101,0.0101);
+
+  //New solver
+  dom_d->auto_ts = false;
+  timestep = (0.4*h/(Cs));
+  dom_d->MechKickDriftSolve(0.0101,1.0e-4);
 
   //First example
   // dom_d->deltat = 1.0e-7;