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
CitcomSmodule.c -rw-r--r-- 1.9 KB
CitcomSmodule.h -rw-r--r-- 1.3 KB
Makefile.am -rw-r--r-- 2.0 KB
advdiffu.c -rw-r--r-- 3.5 KB
advdiffu.h -rw-r--r-- 1.8 KB
bindings.c -rw-r--r-- 9.5 KB
bindings.h -rw-r--r-- 1.2 KB
exceptions.c -rw-r--r-- 1.1 KB
exceptions.h -rw-r--r-- 1.2 KB
getProperty.h -rw-r--r-- 2.9 KB
initial_conditions.c -rw-r--r-- 3.8 KB
initial_conditions.h -rw-r--r-- 2.0 KB
mesher.c -rw-r--r-- 2.2 KB
mesher.h -rw-r--r-- 1.5 KB
misc.c -rw-r--r-- 10.0 KB
misc.h -rw-r--r-- 3.9 KB
outputs.c -rw-r--r-- 2.5 KB
outputs.h -rw-r--r-- 1.6 KB
setProperties.c -rw-r--r-- 26.1 KB
setProperties.h -rw-r--r-- 3.3 KB
stokes_solver.c -rw-r--r-- 7.4 KB
stokes_solver.h -rw-r--r-- 3.1 KB

back to top