Revision 106932597feedfc2ece945dc45c4b0f6122a80d8 authored by Leif Strand on 22 June 2007, 04:36:08 UTC, committed by Leif Strand on 22 June 2007, 04:36:08 UTC
1 parent fcdc999
stokes_solver.c
/*
//~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
//
//<LicenseText>
//
// CitcomS.py by Eh Tan, Eun-seo Choi, and Pururav Thoutireddy.
// Copyright (C) 2002-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>
//
//~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
*/
#include <Python.h>
#include "exceptions.h"
#include "stokes_solver.h"
#include "global_defs.h"
#include "drive_solvers.h"
void assemble_forces(struct All_variables*, int);
void construct_stiffness_B_matrix(struct All_variables*);
void general_stokes_solver(struct All_variables *);
void general_stokes_solver_setup(struct All_variables*);
void get_system_viscosity(struct All_variables*, int, float**, float**);
void set_cg_defaults(struct All_variables*);
void set_mg_defaults(struct All_variables*);
void solve_constrained_flow_iterative(struct All_variables*);
void assemble_forces_pseudo_surf(struct All_variables*, int);
void general_stokes_solver_pseudo_surf(struct All_variables *);
void solve_constrained_flow_iterative_pseudo_surf(struct All_variables*);
char pyCitcom_assemble_forces__doc__[] = "";
char pyCitcom_assemble_forces__name__[] = "assemble_forces";
PyObject * pyCitcom_assemble_forces(PyObject *self, PyObject *args)
{
PyObject *obj;
struct All_variables* E;
if (!PyArg_ParseTuple(args, "O:assemble_forces", &obj))
return NULL;
E = (struct All_variables*)(PyCObject_AsVoidPtr(obj));
assemble_forces(E,0);
Py_INCREF(Py_None);
return Py_None;
}
char pyCitcom_assemble_forces_pseudo_surf__doc__[] = "";
char pyCitcom_assemble_forces_pseudo_surf__name__[] = "assemble_forces_pseudo_surf";
PyObject * pyCitcom_assemble_forces_pseudo_surf(PyObject *self, PyObject *args)
{
PyObject *obj;
struct All_variables* E;
if (!PyArg_ParseTuple(args, "O:assemble_forces_pseudo_surf", &obj))
return NULL;
E = (struct All_variables*)(PyCObject_AsVoidPtr(obj));
assemble_forces_pseudo_surf(E,0);
Py_INCREF(Py_None);
return Py_None;
}
char pyCitcom_construct_stiffness_B_matrix__doc__[] = "";
char pyCitcom_construct_stiffness_B_matrix__name__[] = "construct_stiffness_B_matrix";
PyObject * pyCitcom_construct_stiffness_B_matrix(PyObject *self, PyObject *args)
{
PyObject *obj;
struct All_variables* E;
if (!PyArg_ParseTuple(args, "O:construct_stiffness_B_matrix", &obj))
return NULL;
E = (struct All_variables*)(PyCObject_AsVoidPtr(obj));
construct_stiffness_B_matrix(E);
Py_INCREF(Py_None);
return Py_None;
}
char pyCitcom_general_stokes_solver__doc__[] = "";
char pyCitcom_general_stokes_solver__name__[] = "general_stokes_solver";
PyObject * pyCitcom_general_stokes_solver(PyObject *self, PyObject *args)
{
PyObject *obj;
struct All_variables* E;
if (!PyArg_ParseTuple(args, "O:general_stokes_solver", &obj))
return NULL;
E = (struct All_variables*)(PyCObject_AsVoidPtr(obj));
if(E->control.pseudo_free_surf)
if(E->mesh.topvbc==2)
general_stokes_solver_pseudo_surf(E);
else
assert(0);
else
general_stokes_solver(E);
Py_INCREF(Py_None);
return Py_None;
}
char pyCitcom_general_stokes_solver_setup__doc__[] = "";
char pyCitcom_general_stokes_solver_setup__name__[] = "general_stokes_solver_setup";
PyObject * pyCitcom_general_stokes_solver_setup(PyObject *self, PyObject *args)
{
PyObject *obj;
struct All_variables* E;
if (!PyArg_ParseTuple(args, "O:general_stokes_solver_setup", &obj))
return NULL;
E = (struct All_variables*)(PyCObject_AsVoidPtr(obj));
general_stokes_solver_setup(E);
Py_INCREF(Py_None);
return Py_None;
}
char pyCitcom_get_system_viscosity__doc__[] = "";
char pyCitcom_get_system_viscosity__name__[] = "get_system_viscosity";
PyObject * pyCitcom_get_system_viscosity(PyObject *self, PyObject *args)
{
PyObject *obj;
struct All_variables* E;
if (!PyArg_ParseTuple(args, "O:get_system_viscosity", &obj))
return NULL;
E = (struct All_variables*)(PyCObject_AsVoidPtr(obj));
get_system_viscosity(E,1,E->EVI[E->mesh.levmax],E->VI[E->mesh.levmax]);
Py_INCREF(Py_None);
return Py_None;
}
char pyCitcom_set_cg_defaults__doc__[] = "";
char pyCitcom_set_cg_defaults__name__[] = "set_cg_defaults";
PyObject * pyCitcom_set_cg_defaults(PyObject *self, PyObject *args)
{
PyObject *obj;
struct All_variables* E;
if (!PyArg_ParseTuple(args, "O:set_cg_defaults", &obj))
return NULL;
E = (struct All_variables*)(PyCObject_AsVoidPtr(obj));
E->control.CONJ_GRAD = 1;
set_cg_defaults(E);
Py_INCREF(Py_None);
return Py_None;
}
char pyCitcom_set_mg_defaults__doc__[] = "";
char pyCitcom_set_mg_defaults__name__[] = "set_mg_defaults";
PyObject * pyCitcom_set_mg_defaults(PyObject *self, PyObject *args)
{
PyObject *obj;
struct All_variables* E;
if (!PyArg_ParseTuple(args, "O:set_mg_defaults", &obj))
return NULL;
E = (struct All_variables*)(PyCObject_AsVoidPtr(obj));
E->control.NMULTIGRID = 1;
set_mg_defaults(E);
Py_INCREF(Py_None);
return Py_None;
}
char pyCitcom_set_mg_el_defaults__doc__[] = "";
char pyCitcom_set_mg_el_defaults__name__[] = "set_mg_el_defaults";
PyObject * pyCitcom_set_mg_el_defaults(PyObject *self, PyObject *args)
{
PyObject *obj;
struct All_variables* E;
if (!PyArg_ParseTuple(args, "O:set_mg_el_defaults", &obj))
return NULL;
E = (struct All_variables*)(PyCObject_AsVoidPtr(obj));
E->control.EMULTIGRID = 1;
set_mg_defaults(E);
Py_INCREF(Py_None);
return Py_None;
}
char pyCitcom_solve_constrained_flow_iterative__doc__[] = "";
char pyCitcom_solve_constrained_flow_iterative__name__[] = "solve_constrained_flow_iterative";
PyObject * pyCitcom_solve_constrained_flow_iterative(PyObject *self, PyObject *args)
{
PyObject *obj;
struct All_variables* E;
if (!PyArg_ParseTuple(args, "O:solve_constrained_flow_iterative", &obj))
return NULL;
E = (struct All_variables*)(PyCObject_AsVoidPtr(obj));
solve_constrained_flow_iterative(E);
return Py_BuildValue("d", E->viscosity.sdepv_misfit);
}
char pyCitcom_solve_constrained_flow_iterative_pseudo_surf__doc__[] = "";
char pyCitcom_solve_constrained_flow_iterative_pseudo_surf__name__[] = "solve_constrained_flow_iterative_pseudo_surf";
PyObject * pyCitcom_solve_constrained_flow_iterative_pseudo_surf(PyObject *self, PyObject *args)
{
PyObject *obj;
struct All_variables* E;
if (!PyArg_ParseTuple(args, "O:solve_constrained_flow_iterative_pseudo_surf", &obj))
return NULL;
E = (struct All_variables*)(PyCObject_AsVoidPtr(obj));
solve_constrained_flow_iterative_pseudo_surf(E);
return Py_BuildValue("d", E->viscosity.sdepv_misfit);
}
/* $Id$ */
/* End of file */
![swh spinner](/static/img/swh-spinner.gif)
Computing file changes ...