Changeset 4797
- Timestamp:
- 07/26/10 07:48:41 (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
issm/trunk/src/c/objects/Elements/Penta.cpp
r4780 r4797 2375 2375 double vy[numvertices]; 2376 2376 double vz[numvertices]; 2377 double vel[numvertices]; 2377 2378 double pressure[numvertices]; 2378 2379 int dummy; … … 2387 2388 } 2388 2389 2389 /*Ok, we have vx and vy in values, fill in vx and vyarrays: */2390 /*Ok, we have vx and vy in values, fill in all arrays: */ 2390 2391 for(i=0;i<numvertices;i++){ 2391 2392 vx[i]=values[i*numdofpervertex+0]; … … 2400 2401 pressure[i]=pressure[i]*stokesreconditioning; 2401 2402 } 2403 2404 /*Now Compute vel*/ 2405 for(i=0;i<numvertices;i++) vel[i]=pow( pow(vx[i],2.0) + pow(vy[i],2.0) + pow(vz[i],2.0) , 0.5); 2402 2406 2403 2407 /*Now, we have to move the previous inputs to old … … 2412 2416 this->inputs->AddInput(new PentaVertexInput(VyEnum,vy)); 2413 2417 this->inputs->AddInput(new PentaVertexInput(VzEnum,vz)); 2418 this->inputs->AddInput(new PentaVertexInput(VelEnum,vel)); 2414 2419 this->inputs->AddInput(new PentaVertexInput(PressureEnum,pressure)); 2415 2420 }
Note:
See TracChangeset
for help on using the changeset viewer.