https://github.com/geodynamics/citcoms
Revision ef7299ce82f140ef3b4db2c6e34396c6fef2e814 authored by Eh Tan on 29 October 2008, 23:22:10 UTC, committed by Eh Tan on 29 October 2008, 23:22:10 UTC
* New functions global_v_norm2(), global_p_norm2(), global_div_norm2() compute the L2 norm of velocity, pressure, and divergence. * The L2 norm of the vector f is defined as: \sum_{all elements} ( \int_{element} (f*f dA) ) / A where A is the volume of the domain. * The convergence criterion is controlled by parameter "accuracy" under "[CitcomS.solver.vsolver]". The iteration stops if (norm(div(rho*V))/norm(V) < accuracy) or ( (norm(dV)/norm(V) < accuracy) and (norm(dP)/norm(P) < accuracy) for the last two iterations ) * The default value of "accuracy" becomes 1e-4. * The input parameters "tole_compressibility" and "relative_err_accuracy" under "[CitcomS.solver.vsolver]" are gone. * The mass matrix E->MASS becomes double precision (was single precision).
1 parent 2ba507e
Tip revision: ef7299ce82f140ef3b4db2c6e34396c6fef2e814 authored by Eh Tan on 29 October 2008, 23:22:10 UTC
(for r13196) Fixing convergence criterion for Stokes solver.
(for r13196) Fixing convergence criterion for Stokes solver.
Tip revision: ef7299c
File | Mode | Size |
---|---|---|
CitcomS | ||
bin | ||
etc | ||
examples | ||
lib | ||
module | ||
tests | ||
visual | ||
.gitignore | -rw-r--r-- | 345 bytes |
AUTHORS | -rw-r--r-- | 242 bytes |
COPYING | -rw-r--r-- | 17.6 KB |
ChangeLog | -rw-r--r-- | 408.6 KB |
INSTALL | -rw-r--r-- | 900 bytes |
Makefile.am | -rw-r--r-- | 1.4 KB |
NEWS | -rw-r--r-- | 6.7 KB |
README | -rw-r--r-- | 366 bytes |
configure.ac | -rw-r--r-- | 5.7 KB |
setup.py | -rw-r--r-- | 1.2 KB |
Computing file changes ...