---

utilf.h File Reference

#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include "options.h"

Include dependency graph for utilf.h:

This graph shows which files directly or indirectly include this file:

Go to the source code of this file.

Defines
#define	PI   3.14159265359
#define	sq(x)   ((x) * (x))
#define	cube(x)   ((x) * (x) * (x))
#define	MAX(x, y)   ((x) > (y) ? (x) : (y))
#define	MIN(x, y)   ((x) < (y) ? (x) : (y))
#define	POWER(x, y)   exp((y)*log(x))
#define	realreadSI   doublereadSI
Typedefs
typedef double	real
typedef real *	real1D
typedef real **	real2D
typedef int *	int1D
typedef int **	int2D
Functions
real	findmax (real2D u, int nx, int ny, int *imax, int *jmax)
real	fmodmax (real2D u, int nx, int ny, int *imax, int *jmax)
real	fmean (real2D u, real umean, int xs, int nx, int ny)
real	findmin (real2D u, int nx, int ny, int *imin, int *jmin)
real	sumfield (real2D u, int nx, int ny)
real	findninf (real2D u, real2D v, int nx, int ny, int *imax, int *jmax)
real	dp (real2D p, real2D cc, int nx, int ny)
real	rdivmax (real2D div, int nx, int ny, int *imax, int *jmax)
void	copy (real2D out, real2D in, int nx, int ny)
void	add (real2D a, real2D b, int nx, int ny)
void	fill0 (real2D u, int nx, int ny)

---

Define Documentation

#define cube (  x    )     ((x) * (x) * (x))

Definition at line 17 of file utilf.h. Referenced by bc_vector_bound(), and readglob().

#define MAX (  x, y    )     ((x) > (y) ? (x) : (y))

Definition at line 18 of file utilf.h. Referenced by adaptative().

#define MIN (  x, y    )     ((x) < (y) ? (x) : (y))

Definition at line 19 of file utilf.h. Referenced by computeng().

#define PI   3.14159265359

Definition at line 13 of file utilf.h.

#define POWER (  x, y    )     exp((y)*log(x))

Definition at line 20 of file utilf.h.

#define realreadSI   doublereadSI

Definition at line 29 of file utilf.h. Referenced by ini_impact(), ini_impact2(), inibubble(), inidroplet(), inidroplets(), inifusion(), inirealwave(), initorus(), iniwave(), and readglob().

#define sq (  x    )     ((x) * (x))

Definition at line 16 of file utilf.h. Referenced by advect(), avg_kappa_normals(), bc_vector_bound(), closest_normal(), closest_u(), closest_v(), CubicPolynom(), curvature_write(), distmin(), extra_poly(), extra_velocity_normals(), extra_velocity_normals2(), ini_impact(), ini_impact2(), inijet(), interface_fluxes(), interface_inside(), interface_rrdz(), interface_splines(), interface_surfaces(), interface_vfrac(), interfacial_pressure(), kappa_normals(), momentum(), polyroots(), printnumbers(), readglob(), repulse(), rescale(), and timestep().

---

Typedef Documentation

typedef int* int1D

Definition at line 26 of file utilf.h.

typedef int** int2D

Definition at line 27 of file utilf.h.

typedef double real

Definition at line 22 of file utilf.h. Referenced by adaptative(), adaptative_tau(), advect(), aligned(), avg_kappa_normals(), bc_pressure(), bc_vector(), bc_vector_bound(), bc_vector_div(), bcucof(), bcuint(), bilinu(), bilinv(), closest_normal(), closest_u(), closest_v(), coli(), CubicPolynom(), curvature_write(), cut_interface(), datbvector(), distmin(), dp(), extra_poly(), extra_velocity(), extra_velocity_normals(), extra_velocity_normals2(), findmax(), findmin(), fmean(), fmodmax(), four_aligned(), gauss13(), gauss2(), gauss3(), gnuvectors(), ini_impact(), ini_impact2(), inibubble(), inidroplet(), inidroplets(), inifile(), inifusion(), iniget(), inijet(), inirealwave(), initglobal(), initorus(), initpressure(), iniwave(), interface_arc(), interface_cut_write(), interface_fill(), interface_filter(), interface_fluxes(), interface_hfrac(), interface_inside(), interface_redis(), interface_redis_scale(), interface_rrdz(), interface_rzdr(), interface_splines(), interface_surfaces(), interface_tensionu(), interface_tensionv(), interface_vfrac(), interface_write_markers_speed(), interface_write_splines(), interfaces_area(), interfaces_axivolume(), interfaces_x_amp(), interfaces_xmoment(), interfaces_y_amp(), interfaces_ymoment(), interfacial_pressure(), kappa_normals(), lop(), mgsolve(), minau(), minav(), momentum(), polyderiv(), polyeval(), polyint(), polymul(), polynomial(), polyprint(), polyroot(), polyroots(), pressure(), printfree(), printnumbers(), printxpl(), Pspline(), qtrap(), rdivmax(), readglob(), relax(), repulse(), rescale(), save(), spline(), splint_find(), stresses(), sumfield(), sx_fill(), sx_write(), sy_fill(), sy_write(), timestep(), trapzd(), tridag(), and viscous_tensor().

typedef real* real1D

Definition at line 23 of file utilf.h. Referenced by Pspline(), spline(), splint_find(), and tridag().

typedef real** real2D

Definition at line 24 of file utilf.h. Referenced by adaptative(), adaptative_tau(), add(), advect(), bc_pressure(), bc_pressure_inside(), bc_scalar(), bc_scalar_periodic(), bc_tensor(), bc_tensor_periodic(), bc_vector(), bc_vector_bound(), bc_vector_div(), bc_vector_periodic(), bcucof(), bilinu(), bilinv(), closest_u(), closest_v(), coli(), copy(), datbarray(), datbvector(), divergence(), dp(), extra_poly(), extra_velocity(), extra_velocity_normals(), extra_velocity_normals2(), fill0(), findmax(), findmin(), fmean(), fmodmax(), gnuvectors(), ini_impact(), ini_impact2(), inibubble(), inidroplet(), inidroplets(), iniget(), inijet(), inirealwave(), initorus(), initpressure(), iniwave(), interface_fill(), interface_fluxes(), interface_hfrac(), interface_rrdz(), interface_rzdr(), interface_surfaces(), interface_tensionu(), interface_tensionv(), interface_vfrac(), interface_write_markers_speed(), interfacial_pressure(), lop(), mgfas(), mginit(), mgsolve(), minau(), minav(), momentum(), pressure(), printfree(), printxpl(), rdivmax(), relax(), rescale(), resid(), save(), stresses(), sumfield(), sx_fill(), sx_write(), sy_fill(), sy_write(), and viscous_tensor().

---

Function Documentation

void add (  real2D    a, real2D    b, int    nx, int    ny )

Definition at line 182 of file utilf.c. References real2D. Referenced by mgsolve(), and mgvcycle(). { int i, j; for (i = 1; i <= nx; i++) for (j = 1; j <= ny; j++) a[i][j] += b[i][j]; }

void copy (  real2D    out, real2D    in, int    nx, int    ny )

Definition at line 173 of file utilf.c. References real2D. Referenced by mgfas(), mgsolve(), and pressure(). { int i, j; for (i = 1; i <= nx; i++) for (j = 1; j <= ny; j++) out[i][j] = in[i][j]; }

real dp (  real2D    p, real2D    cc, int    nx, int    ny )

Definition at line 133 of file utilf.c. References real, and real2D. Referenced by timestep(). { int i, j, nin = 0, nout = 0; real pin = 0, pout = 0; for (i = 2; i <= nx - 1; i++) for (j = 2; j <= ny - 1; j++) if (cc[i][j] == 1.0) { pin += p[i][j]; nin++; } else if (cc[i][j] == 0.0) { pout += p[i][j]; nout++; } if (nin > 0 && nout > 0) return pin / (real) nin - pout / (real) nout; else return 0.0; }

void fill0 (  real2D    u, int    nx, int    ny )

Definition at line 192 of file utilf.c. References real2D. Referenced by interface_fill(), interface_fluxes(), interface_hfrac(), interface_rrdz(), interface_rzdr(), interface_surfaces(), interface_tensionu(), interface_tensionv(), interface_vfrac(), mgsolve(), mgvcycle(), and pressure(). { int i, j; for (i = 1; i <= nx; i++) for (j = 1; j <= ny; j++) u[i][j] = 0.0; }

real findmax (  real2D    u, int    nx, int    ny, int *    imax, int *    jmax )

Definition at line 10 of file utilf.c. References real, and real2D. Referenced by timestep(). { int i, j; real res; real fmax; fmax = -1.0; for (i = 2; i <= nx - 1; i++) for (j = 2; j <= ny - 1; j++) { res = u[i][j]; if (res > fmax) { fmax = res; *imax = i; *jmax = j; } } return fmax; }

real findmin (  real2D    u, int    nx, int    ny, int *    imin, int *    jmin )

Definition at line 97 of file utilf.c. References real, and real2D. Referenced by timestep(). { int i, j; real res, fmin; fmin = 1.0; for (i = 2; i <= nx - 1; i++) for (j = 2; j <= ny - 1; j++) { res = u[i][j]; if(res < fmin) { fmin = res; *imin = i; *jmin = j; } } return fmin; }

real findninf (  real2D    u, real2D    v, int    nx, int    ny, int *    imax, int *    jmax )

Referenced by adaptative_tau(), and timestep().

real fmean (  real2D    u, real    umean, int    xs, int    nx, int    ny )

Definition at line 33 of file utilf.c. References real, and real2D. { int i, j; real mean = 0; for (i = xs; i <= nx - 1; i++) for (j = 2; j <= ny - 1; j++) mean += (u[i][j] - umean) * (u[i][j] - umean); return sqrt(mean / (nx - xs) / (ny - 2)); }

real fmodmax (  real2D    u, int    nx, int    ny, int *    imax, int *    jmax )

Definition at line 46 of file utilf.c. References real, real2D, and UNDEFINED. Referenced by adaptative(), printfree(), and timestep(). { int i, j; real res, fmax; fmax = 0.0; for (i = 2; i <= nx - 1; i++) for (j = 2; j <= ny - 1; j++) if (u[i][j] != UNDEFINED) { res = u[i][j]; if(fabs(res) > fmax) { fmax = fabs(res); *imax = i; *jmax = j; } } return fmax; }

real rdivmax (  real2D    div, int    nx, int    ny, int *    imax, int *    jmax )

Definition at line 156 of file utilf.c. References real, and real2D. Referenced by adaptative(), adaptative_tau(), mgsolve(), and timestep(). { int i, j; real val, maxdiv = 0.0; for (i = 2; i < nx; i++) for (j = 2; j < ny; j++) { /* 2D-axi: the divergence is div, the volume is (r + 1/2) */ val = fabs(div[i][j]/((real)j - 1.5)); if (val > maxdiv) { maxdiv = val; *imax = i; *jmax = j; } } return maxdiv; }

real sumfield (  real2D    u, int    nx, int    ny )

Definition at line 119 of file utilf.c. References real, and real2D. Referenced by timestep(). { int i, j; real sum; sum = 0.0; for (i = 2; i <= nx - 1; i++) for (j = 2; j <= ny - 1; j++) sum += u[i][j]; return sum; }

---