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
Citcom.c -rw-r--r-- 5.9 KB
CitcomSFull.c -rw-r--r-- 1.2 KB
CitcomSRegional.c -rw-r--r-- 1.2 KB
Makefile.am -rw-r--r-- 4.1 KB
citcoms.in -rw-r--r-- 1.5 KB
pycitcoms.c -rw-r--r-- 2.7 KB

back to top