---

markers.h File Reference

#include "interpolation.h"

Include dependency graph for markers.h:

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

Go to the source code of this file.

Data Structures
struct	interface
struct	intersec
Defines
#define	NOTENSION   0
#define	TENSION_U   1
#define	TENSION_V   2
#define	NSEC   10
#define	MASSPREC   1.e-3
#define	TRUE_PERIODIC   1
#define	APPROX_PERIODIC   2
#define	AXITRUE   3
#define	AXINATURAL   4
#define	AXIBOTH   5
#define	INTERPREC   1.e-6
#define	MARK_TYPE   0
#define	HOR_TYPE   1
#define	VER_TYPE   2
Functions
int	interface_arc (interface in, real t)
void	interface_write_markers (interface in, FILE *fptr)
void	interface_write_markers_speed (real time, interface in, real2D u, real2D v, real tau, real du, int nx, int ny, FILE *fptr)
void	interface_write_splines (interface in, FILE *fptr, int nb)
void	curvature_write (interface in, FILE *fptr)
real	interfaces_area (interface *in, int n)
real	interfaces_axivolume (interface *in, int n)
real	interfaces_xmoment (interface *in, int n)
real	interfaces_ymoment (interface *in, int n)
real	interfaces_x_amp (interface *in, int n)
real	interfaces_y_amp (interface *in, int n)
void	interface_init (interface *dst)
void	interface_splines (interface in, int nx, int ny)
void	interface_redis (interface *in, real space)
void	interface_redis_scale (interface *in, interface src, real space, real scale)
void	interface_axisymmetry (interface *in, real space, int nx, int ny)
int	reconnect_axis (interface *in, int nx, int ny, real length)
void	interface_copy (interface *dst, interface in)
void	interface_free (interface in)
void	interface_advect (interface in, real2D u, real2D v, int nx, int ny)
void	interface_to_c (interface *in, int nint, real2D work, real2D c, real xo, real yo, int nx, int ny, int sflag, real2D rhox, real2D p, real sigma, real tau)
real	interface_inside (interface in, real x, real y)
void	interface_fill (interface *in, int nin, real2D ap, int nx, int ny)
int	interface_split_axis (interface **inlist, real length)
int	interface_merge_axis (interface **inlist, real length)
real	bilinu (real2D u, real x1, real y1)
real	bilinv (real2D v, real x1, real y1)

---

Define Documentation

#define APPROX_PERIODIC   2

Definition at line 10 of file markers.h. Referenced by interface_splines().

#define AXIBOTH   5

Definition at line 13 of file markers.h. Referenced by ini_impact(), ini_impact2(), inifile(), and interface_splines().

#define AXINATURAL   4

Definition at line 12 of file markers.h. Referenced by interface_splines(), and repulse().

#define AXITRUE   3

Definition at line 11 of file markers.h. Referenced by inifusion(), inijet(), interface_splines(), and repulse().

#define HOR_TYPE   1

Definition at line 17 of file markers.h. Referenced by cut_interface(), interface_cut_write(), interface_fluxes(), interface_surfaces(), and interface_tensionv().

#define INTERPREC   1.e-6

Definition at line 15 of file markers.h. Referenced by cut_interface(), interface_fluxes(), interface_hfrac(), interface_splines(), and interface_vfrac().

#define MARK_TYPE   0

Definition at line 16 of file markers.h. Referenced by cut_interface(), interface_cut_write(), and interface_fluxes().

#define MASSPREC   1.e-3

Definition at line 7 of file markers.h.

#define NOTENSION   0

Definition at line 3 of file markers.h. Referenced by timestep().

#define NSEC   10

Definition at line 6 of file markers.h.

#define TENSION_U   1

Definition at line 4 of file markers.h.

#define TENSION_V   2

Definition at line 5 of file markers.h.

#define TRUE_PERIODIC   1

Definition at line 9 of file markers.h. Referenced by initorus(), and interface_splines().

#define VER_TYPE   2

Definition at line 18 of file markers.h. Referenced by cut_interface(), interface_cut_write(), interface_fluxes(), interface_surfaces(), and interface_tensionu().

---

Function Documentation

real bilinu (  real2D    u, real    x1, real    y1 )

Definition at line 13 of file advect.c. References OUTU, real, real2D, and UNDEFINED. Referenced by advect(). { int i, j; i = x; j = (int) ((real)y - 0.5); x -= (real)i; y -= 0.5 + (real)j; if (OUTU(i,j)) { x = 1. - x; y = 1. - y; if (x + y > 1.0 || OUTU(i,j+1) || OUTU(i+1,j) || OUTU(i+1,j+1)) return UNDEFINED; return x*u[i][j+1] + y*u[i+1][j] + (1. - x - y)*u[i+1][j+1]; } if (OUTU(i+1,j+1)) { if (x + y > 1.0 || OUTU(i+1,j) || OUTU(i,j+1)) return UNDEFINED; return x*u[i+1][j] + y*u[i][j+1] + (1. - x - y)*u[i][j]; } if (OUTU(i+1,j)) { y = 1. - y; if (x + y > 1.0 || OUTU(i,j+1)) return UNDEFINED; return x*u[i+1][j+1] + y*u[i][j] + (1. - x - y)*u[i][j+1]; } if (OUTU(i,j+1)) { x = 1. - x; if (x + y > 1.0) return UNDEFINED; return x*u[i][j] + y*u[i+1][j+1] + (1. - x - y)*u[i+1][j]; } return u[i][j]*(1. - x - y + x*y) + u[i+1][j]*x*(1. - y) + u[i][j+1]*y*(1. - x) + u[i+1][j+1]*x*y; }

real bilinv (  real2D    v, real    x1, real    y1 )

Definition at line 47 of file advect.c. References OUTV, real, real2D, and UNDEFINED. Referenced by advect(). { int i, j; i = (int)((real) x - 0.5); j = y; x -= 0.5 + (real)i; y -= (real) j; if (OUTV(i,j)) { x = 1. - x; y = 1. - y; if (x + y > 1.0 || OUTV(i,j+1) || OUTV(i+1,j) || OUTV(i+1,j+1)) return UNDEFINED; return x*v[i][j+1] + y*v[i+1][j] + (1. - x - y)*v[i+1][j+1]; } if (OUTV(i+1,j+1)) { if (x + y > 1.0 || OUTV(i,j+1) || OUTV(i+1,j)) return UNDEFINED; return x*v[i+1][j] + y*v[i][j+1] + (1. - x - y)*v[i][j]; } if (OUTV(i+1,j)) { y = 1. - y; if (x + y > 1.0 || OUTV(i,j+1)) return UNDEFINED; return x*v[i+1][j+1] + y*v[i][j] + (1. - x - y)*v[i][j+1]; } if (OUTV(i,j+1)) { x = 1. - x; if (x + y > 1.0) return UNDEFINED; return x*v[i][j] + y*v[i+1][j+1] + (1. - x - y)*v[i+1][j]; } return v[i][j]*(1. - x - y + x*y) + v[i+1][j]*x*(1. - y) + v[i][j+1]*y*(1. - x) + v[i+1][j+1]*x*y; }

void curvature_write (  interface    in, FILE *    fptr )

Definition at line 68 of file markers.c. References interface::n, real, splint1, splint2, interface::spx, interface::spy, sq, and interface::t. Referenced by initialize(), and timestep(). { int i; real x1, y1, dl, kappa; fprintf(fptr, "\n"); for (i = 0; i < in.n - 1; i++) { x1 = splint1(in.spx[i], in.t[i]); y1 = splint1(in.spy[i], in.t[i]); dl = sqrt(sq(x1) + sq(y1)); kappa = (x1*splint2(in.spy[i], in.t[i]) - y1*splint2(in.spx[i], in.t[i])) / (dl*dl*dl); fprintf(fptr, "%g %g\n", in.t[i]/in.t[in.n-1], kappa); } }

void interface_advect (  interface    in, real2D    u, real2D    v, int    nx, int    ny )

Referenced by timestep().

int interface_arc (  interface    in, real    t )

Definition at line 10 of file markers.c. References interface::n, real, and interface::t. Referenced by extra_poly(). { int imin = 0, imax = in.n - 1, i; if (t < 0.0 || t > in.t[in.n - 1]) return -1; while(imax - imin > 1) { i = (imin + imax)/2; if (t == in.t[i]) return i; if (t > in.t[i]) imin = i; else imax = i; } return imin; }

void interface_axisymmetry (  interface *    in, real    space, int    nx, int    ny )

Referenced by timestep().

void interface_copy (  interface *    dst, interface    in )

void interface_fill (  interface *    in, int    nin, real2D    ap, int    nx, int    ny )

void interface_free (  interface    in )

Referenced by mgsolve().

void interface_init (  interface *    dst )

Referenced by mgsolve().

real interface_inside (  interface    in, real    x, real    y )

Definition at line 9 of file surfaces.poly.c. { int i; real a, d, dmin = 1e30; for (i = 0; i < in.n - 1; i++) { d = sqrt(sq(x - 0.5*(in.x[i] + in.x[i+1])) + sq(y - 0.5*(in.y[i] + in.y[i+1]))); a = (in.x[i+1] - in.x[i])*(y - in.y[i]) - (x - in.x[i])*(in.y[i+1] - in.y[i]); if (d < fabs(dmin)) dmin = a >= 0.0 ? d : -d; } if (dmin >= 0.0) return dmin; return 0.0; }

int interface_merge_axis (  interface **    inlist, real    length )

void interface_redis (  interface *    in, real    space )

Definition at line 327 of file markers.c. References interface::n, real, splint, interface::spx, interface::spy, interface::t, interface::x, and interface::y. Referenced by initialize(), and timestep(). { int i, j, nb, k = 0; real l, dl, xe, ye; nb = (int)(in->t[in->n-1] / space) + 2; dl = in->t[in->n-1] / (real)(nb - 1); if (dl >= 1.0) printf("WARNING interface_redis(): dl=%g\n", dl); xe = in->x[in->n-1]; /* saving the coordinates of the last point */ ye = in->y[in->n-1]; in->x = (real *)realloc(in->x, nb * sizeof(real)); in->y = (real *)realloc(in->y, nb * sizeof(real)); l = dl; /* Starting at i=1 the first point is unchanged */ for (i = 1; i < nb - 1; i++, l += dl) { /* find the next interval [j;j+1] such as t[j] <= l <= t[j+1] */ while(in->t[k+1] < l) k++; #ifdef DEBUG_REDIS if (k >= in->n-1) { printf("interface_redis: Error i=%d/%d k too large: %d/%d\n", i, nb - 1, k, in->n - 1); exit(1); } if (in->t[k] > l || in->t[k+1] < l) { printf("interface_redis: Error i=%d/%d k=%d, %g <= %g <= %g\n", i, nb - 1, k, in->t[k], l, in->t[k+1]); exit(1); } #endif in->x[i] = splint(in->spx[k], l); in->y[i] = splint(in->spy[k], l); } in->x[nb - 1] = xe; /* The last point is unchanged */ in->y[nb - 1] = ye; /* Copying the saved coordinates */ in->t = (real *)realloc(in->t, nb * sizeof(real)); in->spx = (polynom3 *)realloc(in->spx, (nb - 1) * sizeof(polynom3)); in->spy = (polynom3 *)realloc(in->spy, (nb - 1) * sizeof(polynom3)); #ifdef FREE in->p = (real *)realloc(in->p, nb * sizeof(real)); #endif in->n = nb; }

void interface_redis_scale (  interface *    in, interface    src, real    space, real    scale )

Definition at line 378 of file markers.c. References interface::n, real, splint, interface::spx, interface::spy, interface::t, interface::x, and interface::y. Referenced by mgsolve(). { int i, j, nb, k = 0; real l, dl, xe, ye; nb = (int)(src.t[src.n-1] / space) + 2; dl = src.t[src.n-1] / (real)(nb - 1); in->x = (real *)realloc(in->x, nb * sizeof(real)); in->y = (real *)realloc(in->y, nb * sizeof(real)); in->t = (real *)realloc(in->t, nb * sizeof(real)); in->spx = (polynom3 *)realloc(in->spx, (nb - 1) * sizeof(polynom3)); in->spy = (polynom3 *)realloc(in->spy, (nb - 1) * sizeof(polynom3)); #ifdef FREE in->p = (real *)realloc(in->p, nb * sizeof(real)); #endif in->n = nb; l = dl; /* Starting at i=1 the first point is unchanged */ for (i = 1; i < nb - 1; i++, l += dl) { /* find the next interval [j;j+1] such as t[j] <= l <= t[j+1] */ while(src.t[k+1] < l) k++; in->x[i] = scale*(splint(src.spx[k], l) - 2.) + 2.; in->y[i] = scale*(splint(src.spy[k], l) - 2.) + 2.; } in->x[0] = scale*(src.x[0] - 2.) + 2.; in->y[0] = scale*(src.y[0] - 2.) + 2.; in->x[nb - 1] = scale*(src.x[src.n - 1] - 2.) + 2.; in->y[nb - 1] = scale*(src.y[src.n - 1] - 2.) + 2.; }

void interface_splines (  interface    in, int    nx, int    ny )

Definition at line 246 of file markers.c. References APPROX_PERIODIC, AXIBOTH, AXINATURAL, AXITRUE, interface::bc, INTERPREC, interface::n, PERIODIC_DL, Pspline(), real, SP_DERIVATIVE, SP_NATURAL, spline(), interface::spx, interface::spy, sq, interface::t, TRUE_PERIODIC, interface::x, and interface::y. Referenced by initialize(), interface_filter(), mgsolve(), and timestep(). { real *d2y, l = 0.0, d0, d1, det, dy0, dy1, yp; real frac; int i; d2y = (real *)malloc(in.n * sizeof(real)); in.t[0] = 0.0; for (i = 1; i < in.n; i++) { /* Move slightly the points when they are near the grid lines */ frac = fabs(modf(in.x[i], &d0)); if (frac < INTERPREC || 1. - frac < INTERPREC || fabs(0.5 - frac) < INTERPREC) { printf("WARNING: interfaces_splines(): moving x[%d]: %.10f to %.10f\n", i, in.x[i], in.x[i] + 10.*INTERPREC); in.x[i] += 10.*INTERPREC; } frac = fabs(modf(in.y[i], &d0)); if (frac < INTERPREC || 1. - frac < INTERPREC || fabs(0.5 - frac) < INTERPREC) { printf("WARNING: interfaces_splines(): moving y[%d]: %.10f to %.10f\n", i, in.y[i], in.y[i] + 10.*INTERPREC); in.y[i] += 10.*INTERPREC; } d0 = sqrt(sq(in.x[i] - in.x[i-1]) + sq(in.y[i] - in.y[i-1])); /* if (d0 >= 1.0) printf("WARNING: interfaces_splines(): d0=%f >= 1.0\n", d0); */ l += d0; in.t[i] = l; } switch(in.bc) { case TRUE_PERIODIC: Pspline(in.t, in.x, d2y, in.t[in.n - 1], in.n - 1, in.spx); Pspline(in.t, in.y, d2y, in.t[in.n - 1], in.n - 1, in.spy); break; case APPROX_PERIODIC: /* We need an estimate of the derivative for the periodic ends We suppose a 2nd order interpolation for the points 1, 0, n-2 will do */ d0 = in.t[in.n-1] - in.t[in.n-2]; d1 = d0 + in.t[1] - in.t[0]; det = d1*d0*(d0 - d1); dy0 = in.x[0] - in.x[in.n-2]; dy1 = in.x[1] - in.x[in.n-2]; PERIODIC_DL(dy0, nx); PERIODIC_DL(dy1, nx); yp = (2.*(dy0*d1 - dy1*d0)*d0 + dy1*d0*d0 - dy0*d1*d1) / det; spline(in.t, in.x, in.n, yp, yp, d2y, in.spx, SP_DERIVATIVE); dy0 = in.y[0] - in.y[in.n-2]; dy1 = in.y[1] - in.y[in.n-2]; yp = (2.*(dy0*d1 - dy1*d0)*d0 + dy1*d0*d0 - dy0*d1*d1) / det; spline(in.t, in.y, in.n, yp, yp, d2y, in.spy, SP_DERIVATIVE); break; case AXITRUE: spline(in.t, in.x, in.n, 0.0, 0.0, d2y, in.spx, SP_DERIVATIVE); spline(in.t, in.y, in.n, 1.0, -1.0, d2y, in.spy, SP_DERIVATIVE); break; case AXINATURAL: spline(in.t, in.x, in.n, 0.0, 0.0, d2y, in.spx, SP_NATURAL); spline(in.t, in.y, in.n, 0.0, 0.0, d2y, in.spy, SP_NATURAL); break; case AXIBOTH: spline(in.t, in.x, in.n, 0.0, 0.0, d2y, in.spx, SP_DERIVATIVE); spline(in.t, in.y, in.n, 1.0, 1.0, d2y, in.spy, SP_DERIVATIVE); break; default: printf("interface_splines(): bad value for the bc parameter: bc = %d\n", in.bc); exit(1); } free(d2y); }

int interface_split_axis (  interface **    inlist, real    length )

void interface_to_c (  interface *    in, int    nint, real2D    work, real2D    c, real    xo, real    yo, int    nx, int    ny, int    sflag, real2D    rhox, real2D    p, real    sigma, real    tau )

Referenced by timestep().

void interface_write_markers (  interface    in, FILE *    fptr )

Definition at line 26 of file markers.c. References interface::n, interface::x, and interface::y. Referenced by initialize(), pressure(), and timestep(). { int i; fprintf(fptr, "\n"); for (i = 0; i <= in.n-1; i++) fprintf(fptr, "%g %g\n", in.x[i], in.y[i]); fprintf(fptr, "&\n"); }

void interface_write_markers_speed (  real    time, interface    in, real2D    u, real2D    v, real    tau, real    du, int    nx, int    ny, FILE *    fptr )

Definition at line 36 of file markers.c. References interface::n, real, real2D, interface::x, and interface::y. { /* Physical velocity*/ int i; real h=1./(nx-2); fprintf(fptr, "\n"); for (i = 0; i <= in.n-1; i++) fprintf(fptr, "%g %g %g %g %g\n", time, in.x[i], in.y[i], bilinu(u, in.x[i], in.y[i])*du*h/tau, bilinu(u, in.x[i], in.y[i])*du*h/tau); fprintf(fptr, "&\n"); }

void interface_write_splines (  interface    in, FILE *    fptr, int    nb )

Definition at line 52 of file markers.c. References interface::n, real, splint_find(), interface::spx, interface::spy, and interface::t. Referenced by closest_normal(), initialize(), and timestep(). { int i; real t; fprintf(fptr, "\n"); for (i = 0; i < nb; i++) { t = (real) i * in.t[in.n-1] / (real) (nb - 1); fprintf(fptr, "%g %g\n", splint_find(in.t, in.spx, in.n, t), splint_find(in.t, in.spy, in.n, t)); } fprintf(fptr, "&\n"); }

real interfaces_area (  interface *    in, int    n )

Definition at line 85 of file markers.c. References polynom3::a, polynom3::b, polynom3::c, polynom3::d, interface::n, real, interface::spx, interface::spy, and interface::t. { int i, j; real ar = 0.0, a1, a2, a3, a4, a5, a6, t2, t1; interface in1; for (j = 0; j < n; j++) { in1 = in[j]; for (i = 0; i < in1.n - 1; i++) { a1 = in1.spx[i].a*in1.spy[i].b; a2 = 0.5*(in1.spx[i].b*in1.spy[i].b + 2.*in1.spx[i].a*in1.spy[i].c); a3 = (in1.spx[i].c*in1.spy[i].b + 2.*in1.spx[i].b*in1.spy[i].c + 3.*in1.spx[i].a*in1.spy[i].d)/3.; a4 = 0.25*(in1.spx[i].d*in1.spy[i].b + 2.*in1.spx[i].c*in1.spy[i].c + 3.*in1.spx[i].b*in1.spy[i].d); a5 = (3.*in1.spx[i].c*in1.spy[i].d + 2.*in1.spx[i].d*in1.spy[i].c)/5.; a6 = 0.5*in1.spx[i].d*in1.spy[i].d; t1 = in1.t[i]; t2 = in1.t[i+1]; ar += t2*(a1 + t2*(a2 + t2*(a3 + t2*(a4 + t2*(a5 + t2*a6))))) - t1*(a1 + t1*(a2 + t1*(a3 + t1*(a4 + t1*(a5 + t1*a6))))); } } return ar; }

real interfaces_axivolume (  interface *    in, int    n )

Definition at line 112 of file markers.c. References polynom3::a, polynom3::b, polynom3::c, polynom3::d, interface::n, real, interface::spx, interface::spy, interface::t, interface::x, and interface::y. Referenced by timestep(). { int i, j; real ar = 0.0, a1, a2, a3, a4, a5, a6, a7, a8, a9, t2, t1; real ax, bx, cx, dx, ay, by, cy, dy; real x1, y1, x2, y2; interface in1; for (j = 0; j < n; j++) { in1 = in[j]; for (i = 0; i < in1.n - 1; i++) { #ifdef POLYVOLUME /* integration of xy dy along the polynomial segment */ ax = in1.spx[i].a; bx = in1.spx[i].b; cx = in1.spx[i].c; dx = in1.spx[i].d; ay = in1.spy[i].a - 2.0; by = in1.spy[i].b; cy = in1.spy[i].c; dy = in1.spy[i].d; a1 = ax*ay*by; a2 = ax*ay*cy+(ax*by+bx*ay)*by/2.; a3 = ax*ay*dy+2.0/3.0*(ax*by+bx*ay)*cy+(ax*cy+bx*by+cx*ay)*by/3.; a4 = 3.0/4.0*(ax*by+bx*ay)*dy+(ax*cy+bx*by+cx*ay)*cy/2. +(ax*dy+bx*cy+cx*by+dx*ay)*by/4.; a5 = 3.0/5.0*(ax*cy+bx*by+cx*ay)*dy +2.0/5.0*(ax*dy+bx*cy+cx*by+dx*ay)*cy+(bx*dy+cx*cy+dx*by)*by/5.; a6 = (ax*dy+bx*cy+cx*by+dx*ay)*dy/2+(bx*dy+cx*cy+dx*by)*cy/3 +(cx*dy+dx*cy)*by/6.; a7 = 3.0/7.0*(bx*dy+cx*cy+dx*by)*dy+2.0/7.0*(cx*dy+dx*cy)*cy +dx*dy*by/7.; a8 = 3.0/8.0*(cx*dy+dx*cy)*dy+dx*dy*cy/4.; a9 = dx*dy*dy/3.; t1 = in1.t[i]; t2 = in1.t[i+1]; ar += t2*(a1 + t2*(a2 + t2*(a3 + t2*(a4 + t2*(a5 + t2*(a6 + t2*(a7 + t2*(a8 + t2*a9)))))))) - t1*(a1 + t1*(a2 + t1*(a3 + t1*(a4 + t1*(a5 + t1*(a6 + t1*(a7 + t1*(a8 + t1*a9)))))))); #else /* integration of xy dy along the segment */ x1 = in1.x[i]; y1 = in1.y[i] - 2.; x2 = in1.x[i+1]; y2 = in1.y[i+1] - 2.; ar += ((x2 - x1)*(y2 - y1)/3. + (y1*(x2 - x1) + x1*(y2 - y1))/2. + x1*y1)*(y2 - y1); #endif } } return 2.*PI*ar; }

real interfaces_x_amp (  interface *    in, int    n )

Definition at line 218 of file markers.c. References interface::n, real, and interface::x. { int i, j; real max = -1e6, min = 1e6; for (j = 0; j < n; j++) for (i = 0; i < in[j].n; i++) { max = in[j].x[i] > max ? in[j].x[i] : max; min = in[j].x[i] < min ? in[j].x[i] : min; } return fabs(max - min); }

real interfaces_xmoment (  interface *    in, int    n )

Definition at line 160 of file markers.c. References interface::n, real, splint, splint1, interface::spx, interface::spy, and interface::t. { /* compute Sum{x^2.dx.dy} */ /* coefficients of Newton-Cotes formula */ static real coef[7] = {1., 3., 3., 2., 3., 3., 1.}; int i, j, k; real x, yp, h, t1, t2, t, sum, total = 0.0; interface in1; polynom3 px, py; for (k = 0; k < n; k++) { in1 = in[k]; for (j = 0; j < in1.n - 1; j++) { t1 = in1.t[j]; t2 = in1.t[j+1]; px = in1.spx[j]; py = in1.spy[j]; sum = 0.0; h = (t2 - t1)/6.; for (i = 0; i < 7; i++) { t = t1 + h*i; x = splint(px, t) - 2.; yp = splint1(py, t); sum += coef[i]*x*x*x*yp; } total += 3.*h*sum/8.; } } return total/3.; }

real interfaces_y_amp (  interface *    in, int    n )

Definition at line 232 of file markers.c. References interface::n, real, and interface::y. { int i, j; real max = -1e6, min = 1e6; for (j = 0; j < n; j++) for (i = 0; i < in[j].n; i++) { max = in[j].y[i] > max ? in[j].y[i] : max; min = in[j].y[i] < min ? in[j].y[i] : min; } return fabs(max - min); }

real interfaces_ymoment (  interface *    in, int    n )

Definition at line 189 of file markers.c. References interface::n, real, splint, splint1, interface::spx, interface::spy, and interface::t. { /* compute Sum{y^2.dx.dy} */ /* coefficients of Newton-Cotes formula */ static real coef[7] = {1., 3., 3., 2., 3., 3., 1.}; int i, j, k; real y, xp, h, t1, t2, t, sum, total = 0.0; interface in1; polynom3 px, py; for (k = 0; k < n; k++) { in1 = in[k]; for (j = 0; j < in1.n - 1; j++) { t1 = in1.t[j]; t2 = in1.t[j+1]; px = in1.spx[j]; py = in1.spy[j]; sum = 0.0; h = (t2 - t1)/6.; for (i = 0; i < 7; i++) { t = t1 + h*i; y = splint(py, t) - 2.; xp = splint1(px, t); sum += coef[i]*y*y*y*xp; } total += 3.*h*sum/8.; } } return total/3.; }

int reconnect_axis (  interface *    in, int    nx, int    ny, real    length )

Referenced by timestep().

---