https://github.com/geodynamics/citcoms
Revision dc9672a15e02a96f8eaca7205d7c1786c3212a62 authored by Eh Tan on 20 March 2009, 23:36:56 UTC, committed by Eh Tan on 20 March 2009, 23:36:56 UTC
A new input parameter solver.vsolver.max_mg_cycles (default to 50) controls the max. number of multigrid cycles when solving the velocity. Without this parameter, when the stiffness matrix is near singular, the the multigrid solver may become an infinite loop.
1 parent 25a7493
Tip revision: dc9672a15e02a96f8eaca7205d7c1786c3212a62 authored by Eh Tan on 20 March 2009, 23:36:56 UTC
limiting the number of multigrid cycles in solve_del2_u().
limiting the number of multigrid cycles in solve_del2_u().
Tip revision: dc9672a
advection_diffusion.h
/*
*~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
*
*<LicenseText>
*
* CitcomS by Louis Moresi, Shijie Zhong, Lijie Han, Eh Tan,
* Clint Conrad, Michael Gurnis, and Eun-seo Choi.
* Copyright (C) 1994-2005, California Institute of Technology.
*
* This program is free software; you can redistribute it and/or modify
* it under the terms of the GNU General Public License as published by
* the Free Software Foundation; either version 2 of the License, or
* (at your option) any later version.
*
* This program is distributed in the hope that it will be useful,
* but WITHOUT ANY WARRANTY; without even the implied warranty of
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
* GNU General Public License for more details.
*
* You should have received a copy of the GNU General Public License
* along with this program; if not, write to the Free Software
* Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
*
*</LicenseText>
*
*~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
*/
#if !defined(CitcomS_advection_diffusion_h)
#define CitcomS_advection_diffusion_h
#ifdef __cplusplus
extern "C" {
#endif
void PG_timestep_init(struct All_variables *);
void PG_timestep_solve(struct All_variables *);
void std_timestep(struct All_variables *);
#ifdef __cplusplus
}
#endif
#endif
Computing file changes ...
