Skip to content

Conversation

@EamonnGlynn
Copy link

set and get strain rate functions have been added, as well as code to test it.

Copy link

@onkarsahni onkarsahni left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Looks good, some changes needed.

h_strainRate(i,3) = xyShear[i];
h_strainRate(i,4) = xzShear[i];
h_strainRate(i,5) = yzShear[i];
}

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Need to discuss this numbering, i.e., for 6 components of symmetric tensor per vertex (i.e., 0 through 5).

Copy link
Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

function now uses one 2d array instead of six 1d arrays.

src/pmpo_c.cpp Outdated
Kokkos::deep_copy(arrayHost, array_d);
}

void polympo_setMeshVtxStrainRate_f(MPMesh_ptr p_mpmesh, const int nVertices, const double *xNormal, const double *yNormal, const double *zNormal, const double *xyShear, const double *xzShear, const double *yzShear){

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Need to discuss if each component (of the 6 components of symmetric tensor) will be given as a separate array.

Copy link
Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

It has been converted to take one 2d array instead

call assert((yArray(i) .eq. value2-i), "Assert MeshVel v-component Velocity Fail")
end do

!VtxStrainRate needs 6 arrays, using x,y,zArrays for all 6

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Only in this file, change 2d array Mesharray name to Mesharray_vel and use a new 2d array Mesharray_strainrate(numCompsStrainRate,nverts) for "VtxStrainRate" testing with 6 values per vertex

Copy link
Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Changes have been made.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

3 participants