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
File Mode Size
Makefile.am -rw-r--r-- 1.1 KB
__vault__.odb -rw-r--r-- 1.1 KB
full-sphere.odb -rw-r--r-- 1.2 KB
full.odb -rw-r--r-- 1.2 KB
incomp-newtonian.odb -rw-r--r-- 1.2 KB
incomp-non-newtonian.odb -rw-r--r-- 1.2 KB
regional-sphere.odb -rw-r--r-- 1.2 KB
regional.odb -rw-r--r-- 1.2 KB
temp.odb -rw-r--r-- 1.2 KB

back to top