https://github.com/geodynamics/citcoms
Revision adf8b195539123ce36fff2eb604485eb82c6ecf5 authored by Thorsten Becker on 21 November 2009, 00:15:57 UTC, committed by Thorsten Becker on 21 November 2009, 00:15:57 UTC
convergence at finite compressibility > accuracy. Changed criterion to loop while (convergence < 2) OR (compressibility < imp) Please review Also modified the sloppy iteration branch, where pressure isn't checked.
1 parent 5f38ecb
Tip revision: adf8b195539123ce36fff2eb604485eb82c6ecf5 authored by Thorsten Becker on 21 November 2009, 00:15:57 UTC
Exit criterion for incompressibility iteration allowed for velocity and pressure
Exit criterion for incompressibility iteration allowed for velocity and pressure
Tip revision: adf8b19
h5util.h
/*
* h5util.h by Luis Armendariz and Eh Tan.
* Copyright (C) 2006, 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
*/
#include "hdf5.h"
typedef struct field_t
{
const char *name;
int rank;
hsize_t *dims;
hsize_t *maxdims;
hsize_t *offset;
hsize_t *count;
int n;
float *data;
} field_t;
field_t *open_field(hid_t group, const char *name);
herr_t read_field(hid_t group, field_t *field, int cap);
herr_t close_field(field_t *field);
herr_t get_attribute_str(hid_t obj_id, const char *attr_name, char **data);
herr_t get_attribute_int(hid_t input, const char *name, int *val);
herr_t get_attribute_float(hid_t input, const char *name, float *val);
Computing file changes ...