Converting between nodal and elemental values

ryanwalker

Is there an easy way to get an elemental average of a variable defined on the vertices? (velocity, thickness, etc.) I can't figure out whether averaging() can do this.

mathieumorlighem

Hi Ryan,

to get the element average of a nodal vector, you can use the following command:


S=md.geometry.surface(md.mesh.elements)*repmat(1,size(md.mesh.elements,2),1)/size(md.mesh.elements,2);

where S is the element average of the surface topography md.geometry.surface.

To interpolate an element based vector onto a nodal vector, you can use averaging:


N=averaging(md,md.materials.rheology_n,0);

ryanwalker

Thanks! It might be nice to add the nodes-to-elements solution as a function.