https://github.com/geodynamics/citcoms
Revision 38bdaf9dd9be244384bb4889b7dd6ce0834035d0 authored by Eh Tan on 25 January 2007, 23:30:49 UTC, committed by Eh Tan on 25 January 2007, 23:30:49 UTC
The previous code computed and stored the Cartesian velocity in E->temp in
v_from_vector() and accessed E->temp in velo_from_element. This is bug prone.
E->temp is a temporary array and should be used for storage.

Also, in previous code, if E->sphere.cap[m].V is modified, one has to recompute
E->temp. Otherwise, the Cartesian velocity will be incorrect.

1 parent 8734680
History
Tip revision: 38bdaf9dd9be244384bb4889b7dd6ce0834035d0 authored by Eh Tan on 25 January 2007, 23:30:49 UTC
Compute the velocity in Cartesian coord on the fly.
Tip revision: 38bdaf9

back to top