Flatten mesh to single Object due to access errors

[luchete80]

There are not neighbors

[luchete80]

All IDs are zero

[luchete80]

Problem iis in CalculateSurface

[luchete80]

• It will be ommited SetMeshVelKernel (which sets Mesh_d::m_v)

• now will be used node_v (which sets the above function)

• For setting node_v will be used SetNodesVelKernel<<>> which calls to

global inline void SetNodesVelKernel(TriMesh_d mesh_d,const double3 &v){
//int i = threadIdx.x + blockDim.xblockIdx.x;
for (int i=0;mesh_d->nodecount;i++)
//if (i<mesh_d->nodecount)
mesh_d->node_v[i] = v;
}

         cudaDeviceSynchronize();  
            cout << "Mesh ID " << id << ", "<< "Velocity set to : "
                                        <<dom_d->bConds[bc].value.x 
                                        <<", " <<dom_d->bConds[bc].value.y 
                                        <<", " << dom_d->bConds[bc].value.z<<endl;

AND Mesh::Move is changed by ommitting node_v calculation from m_v
'''

inline device void TriMesh_d::Move(double dt){
//printf("node count %d\n", nodecount);
int n = threadIdx.x + blockDim.x*blockIdx.x; //Parallelize by node
if ( n < nodecount ){

double3 vr 	= cross(m_w, node[n]);
//// FOR THE NEW FLATTENED MESH IS NT ANYMORE AN SCALAR M_V
///node_v[n] = m_v + vr; //MUST BE DELETED, SET PREVIOUSLY

//printf("NODE V: %f %f %f \n",node_v[n].x, node_v[n].y, node_v[n].z);

'''

[luchete80]

Inside Solver are Functions:

MeshUpdateKernel (which calls Move() )

UpdateContactParticlesKernel() which calls mesh_id

//TODO: CHANGE TO SEVERAL CONTACT SURFACES 
inline void __device__ Domain_d::UpdateContactParticles(int mesh_id){
  
  int e = threadIdx.x + blockDim.x*blockIdx.x;	
  //printf("elemn %d\n", e);
  
  if (e < trimesh[mesh_id]->elemcount) {
    //printf("UPDATING e %d\n",e);
    //int e = element[i];
    double3 vv = make_double3(0.);
    for (int en = 0; en<3; en++){
      //printf("particle %d \n",i);
      //printf ("node %d \n",trimesh->elnode[3*e+en]);
      // if (trimesh->elnode[3*e+en] < trimesh->nodecount )
      vv += trimesh[mesh_id] -> node_v[trimesh[mesh_id]->elnode[3*e+en]];
      // else 
        // printf("error \n");
    }
    

    v [first_fem_particle_idx[mesh_id] + e] = vv/3.;
    a [first_fem_particle_idx[mesh_id] + e] = make_double3(0.);
    //printf(" particle %d , v %f %f %f \n", e, vv.x, vv.y, vv.z);
    // if (length(normal[e])<1.0e-3)
      // printf("UPDATING ERROR ZERO mesh normal, %f %f %f\n", trimesh -> normal[e].x , trimesh -> normal[e].y, trimesh -> normal[e].z);
    normal[first_fem_particle_idx[mesh_id] + e] = trimesh[mesh_id] -> normal[e];
    //printf("mesh normal, %f %f %f\n", trimesh -> normal[e].x , trimesh -> normal[e].y, trimesh -> normal[e].z);
  }
}