#include <stddef.h>
#include <stdlib.h>
#include <string.h>
#include <math.h>
#include <float.h>
#include "c99.h"
#include "name.h"
#include "fail.h"
#include "types.h"
#include "mem.h"
#include "tensor.h"
#include "poly.h"

Include dependency graph for findpts_el_3.c:

Classes
struct	findpts_el_pt_3

struct	findpts_el_gface_3

struct	findpts_el_gedge_3

struct	findpts_el_gpt_3

struct	findpts_el_data_3

Macros
#define	findpts_el_setup_3 PREFIXED_NAME(findpts_el_setup_3)

#define	findpts_el_free_3 PREFIXED_NAME(findpts_el_free_3 )

#define	findpts_el_3 PREFIXED_NAME(findpts_el_3 )

#define	findpts_el_eval_3 PREFIXED_NAME(findpts_el_eval_3 )

#define	DIAGNOSTICS_ITERATIONS 0

#define	CONVERGED_FLAG (1u<<6)

#define	FLAG_MASK 0x7fu

#define	DO_MAX(x)

#define	SET_FACE(i, base, n)

#define	SET_EDGE1(j, k, d, rd, rn, base)

#define	SET_EDGE2(j, d, rd, rn, base)

#define	SET_EDGE(j, rd, rn, base)

#define	DXDN1(i, d) (fd->sides+base+(9(i) +(d))n)

#define	D2XDN1(i, d) (fd->sides+base+(9(i)+3+(d))n)

#define	D2XDN2(i, d) (fd->sides+base+(9(i)+6+(d))n)

#define	CHECK_CONSTRAINT(drcd, d3)

#define	SETDR(d)

#define	SETR(d)

#define	EVAL(r, s) -(y[0]r+y[1]s)+(rA[0]r+(2rA[1]+sA[2])s)/2

#define	SETR(d, d3)

#define	EVAL(dr) (drA-2y)*dr

Typedefs
typedef void	compute_face_data_fun(struct findpts_el_data_3 *fd)

typedef void	findpt_fun(struct findpts_el_pt_3 const out, struct findpts_el_data_3 const fd, const struct findpts_el_pt_3 *const p, const unsigned pn, const double tol)

Functions
static void	lin_solve_3 (double x[3], const double A[9], const double y[3])

static void	lin_solve_sym_2 (double x[2], const double A[3], const double y[2])

static unsigned	num_constrained (const unsigned flags)

static unsigned	pt_flags_to_bin_noC (const unsigned flags)

static unsigned	pt_flags_to_bin (const unsigned flags)

static unsigned	plus_1_mod_3 (const unsigned x)

static unsigned	plus_2_mod_3 (const unsigned x)

static unsigned	which_bit (const unsigned x)

static unsigned	face_index (const unsigned x)

static unsigned	edge_index (const unsigned x)

static unsigned	point_index (const unsigned x)

static unsigned	work_size (const unsigned nr, const unsigned ns, const unsigned nt, const unsigned npt_max)

void	findpts_el_setup_3 (struct findpts_el_data_3 *const fd, const unsigned n[3], const unsigned npt_max)

void	findpts_el_free_3 (struct findpts_el_data_3 *const fd)

static void	compute_face_data_rs (struct findpts_el_data_3 *fd)

static void	compute_face_data_st (struct findpts_el_data_3 *fd)

static void	compute_face_data_tr (struct findpts_el_data_3 *fd)

static const struct findpts_el_gface_3 *	get_face (struct findpts_el_data_3 *fd, unsigned fi)

static void	compute_edge_data (struct findpts_el_data_3 *fd, unsigned d)

static const struct findpts_el_gedge_3 *	get_edge (struct findpts_el_data_3 *fd, unsigned ei)

static void	compute_pt_data (struct findpts_el_data_3 *fd)

static const struct findpts_el_gpt_3 *	get_pt (struct findpts_el_data_3 *fd, unsigned pi)

static int	reject_prior_step_q (struct findpts_el_pt_3 const out, const double resid[3], const struct findpts_el_pt_3 const p, const double tol)

static void	newton_vol (struct findpts_el_pt_3 const out, const double jac[9], const double resid[3], const struct findpts_el_pt_3 const p, const double tol)

static void	newton_face (struct findpts_el_pt_3 const out, const double jac[9], const double rhes[3], const double resid[3], const unsigned d1, const unsigned d2, const unsigned dn, const unsigned flags, const struct findpts_el_pt_3 const p, const double tol)

static void	newton_edge (struct findpts_el_pt_3 const out, const double jac[9], const double rhes, const double resid[3], const unsigned de, const unsigned dn1, const unsigned dn2, unsigned flags, const struct findpts_el_pt_3 const p, const double tol)

static void	findpt_vol (struct findpts_el_pt_3 const out, struct findpts_el_data_3 const fd, const struct findpts_el_pt_3 *const p, const unsigned pn, const double tol)

static void	findpt_face (struct findpts_el_pt_3 const out, struct findpts_el_data_3 const fd, const struct findpts_el_pt_3 *const p, const unsigned pn, const double tol)

static void	findpt_edge (struct findpts_el_pt_3 const out, struct findpts_el_data_3 const fd, const struct findpts_el_pt_3 *const p, const unsigned pn, const double tol)

static void	findpt_pt (struct findpts_el_pt_3 const out, struct findpts_el_data_3 const fd, const struct findpts_el_pt_3 *const p, const unsigned pn, const double tol)

static void	seed (struct findpts_el_data_3 const fd, struct findpts_el_pt_3 const pt, const unsigned npt)

void	findpts_el_3 (struct findpts_el_data_3 *const fd, const unsigned npt, const double tol)

void	findpts_el_eval_3 (double const out_base, const unsigned out_stride, const double const r_base, const unsigned r_stride, const unsigned pn, const double const in, struct findpts_el_data_3 const fd)

Macro Definition Documentation

#define CHECK_CONSTRAINT	(	drcd,
		d3
	)

Value:

do { \
      const double rz = r0[d3]+dr[d3], lb=bnd[2*d3],ub=bnd[2*d3+1]; \
      const double delta=drcd, nr = r0[d3]+(dr[d3]+delta); \
      if((nr-lb)*(ub-nr)<0) { \
        if(nr<lb) { \
          double f = (lb-rz)/delta; \
          if(f<fac) fac=f, new_flags = 1u<<(2*d3); \
        } else { \
          double f = (ub-rz)/delta; \
          if(f<fac) fac=f, new_flags = 2u<<(2*d3); \
        } \
      } \
    } while(0)

Referenced by newton_vol().

#define CONVERGED_FLAG (1u<<6)

Definition at line 72 of file findpts_el_3.c.

Referenced by findpt_pt(), newton_edge(), newton_face(), newton_vol(), and reject_prior_step_q().

#define D2XDN1	(	i,
		d
	)	(fd->sides+base+(9(i)+3+(d))n)

Referenced by compute_edge_data().

#define D2XDN2	(	i,
		d
	)	(fd->sides+base+(9(i)+6+(d))n)

Referenced by compute_edge_data().

#define DIAGNOSTICS_ITERATIONS 0

Definition at line 22 of file findpts_el_3.c.

#define DO_MAX ( x )

Value:

do { const unsigned temp=(x); \

wsize=temp>wsize?temp:wsize; } while(0)

x

#define x

Referenced by work_size().

#define DXDN1	(	i,
		d
	)	(fd->sides+base+(9(i) +(d))n)

Referenced by compute_edge_data().

#define EVAL	(	r,
		s
	)	-(y[0]r+y[1]s)+(rA[0]r+(2rA[1]+sA[2])s)/2

Referenced by newton_edge(), and newton_face().

#define EVAL ( dr ) (dr*A-2*y)*dr

#define findpts_el_3 PREFIXED_NAME(findpts_el_3 )

Definition at line 16 of file findpts_el_3.c.

#define findpts_el_eval_3 PREFIXED_NAME(findpts_el_eval_3 )

Definition at line 17 of file findpts_el_3.c.

#define findpts_el_free_3 PREFIXED_NAME(findpts_el_free_3 )

Definition at line 15 of file findpts_el_3.c.

#define findpts_el_setup_3 PREFIXED_NAME(findpts_el_setup_3)

Definition at line 14 of file findpts_el_3.c.

#define FLAG_MASK 0x7fu

Definition at line 73 of file findpts_el_3.c.

Referenced by findpt_edge(), findpt_face(), findpt_pt(), findpts_el_3(), and reject_prior_step_q().

#define SET_EDGE	(	j,
		rd,
		rn,
		base
	)

Value:

do { \
    for(d=0;d<3;++d) { SET_EDGE1(j,plus_2_mod_3(j),d,rd,rn,base); \
                       SET_EDGE2(j,d,rd,rn,base); } \
  } while(0)

Referenced by findpts_el_setup_3().

#define SET_EDGE1	(	j,
		k,
		d,
		rd,
		rn,
		base
	)

Value:

for(i=0;i<2;++i) \
      fd->edge[4*j+2*i+0].dxdn2[d] = fd->face[2*k+i].dxdn[d], \
      fd->edge[4*j+2*i+1].dxdn2[d] = fd->face[2*k+i].dxdn[d]+n##rd##rn-n##rd;

#define SET_EDGE2	(	j,
		d,
		rd,
		rn,
		base
	)

Value:

for(i=0;i<4;++i) \
      fd->edge[4*j+i].dxdn1 [d] = fd->sides + base + (9*i  +d)*n##rd, \
      fd->edge[4*j+i].d2xdn1[d] = fd->sides + base + (9*i+3+d)*n##rd, \
      fd->edge[4*j+i].d2xdn2[d] = fd->sides + base + (9*i+6+d)*n##rd;

#define SET_FACE	(	i,
		base,
		n
	)

Value:

do { for(d=0;d<3;++d) \
    fd->face[2*i  ].x[d]    = fd->sides + base +    d *n, \
    fd->face[2*i  ].dxdn[d] = fd->sides + base + (3+d)*n, \
    fd->face[2*i+1].x[d]    = fd->sides + base + (6+d)*n, \
    fd->face[2*i+1].dxdn[d] = fd->sides + base + (9+d)*n; \
  } while(0)

Referenced by findpts_el_setup_3().

#define SETDR ( d )

Value:

do { \
      unsigned f = flags>>(2*d) & 3u; \
      if(f) dr[d] = bnd[2*d+(f-1)] - r0[d]; \
    } while(0)

Referenced by newton_vol().

#define SETR ( d )

Value:

do { \
    unsigned f = flags>>(2*d) & 3u; \
    out->r[d] = f==0 ? r0[d]+dr[d] : ( f==1 ? -1 : 1 ); \
  } while(0)

Referenced by newton_face(), and newton_vol().

#define SETR	(	d,
		d3
	)

Value:

do { \
    const unsigned f = i>>(2*d) & 3u; \
    if(f==0) r[d]=r0[d]+dr[d]; \
    else { \
      if((f&(mask>>(2*d)))==0) r[d]=r0[d]+(f==1?-tr:tr); \
      else                     r[d]=(f==1?-1:1), new_flags |= f<<(2*d3); \
    } \
  } while(0)

Typedef Documentation

typedef void compute_face_data_fun(struct findpts_el_data_3 *fd)

Definition at line 289 of file findpts_el_3.c.

typedef void findpt_fun(struct findpts_el_pt_3 *const out, struct findpts_el_data_3 *const fd, const struct findpts_el_pt_3 *const p, const unsigned pn, const double tol)

Definition at line 915 of file findpts_el_3.c.

Function Documentation

static void compute_edge_data	(	struct findpts_el_data_3 *	fd,
		unsigned	d
	)

static

Definition at line 364 of file findpts_el_3.c.

References D2XDN1, D2XDN2, DXDN1, findpts_el_data_3::edge, get_face(), i, n, findpts_el_data_3::n, nr, ns, plus_1_mod_3(), plus_2_mod_3(), tensor_mtxm, tensor_mxm(), work, findpts_el_data_3::work, findpts_el_data_3::wtend, and x.

Referenced by get_edge().

Here is the call graph for this function:

Here is the caller graph for this function:

static void compute_face_data_rs ( struct findpts_el_data_3 * fd )

static

Definition at line 292 of file findpts_el_3.c.

References findpts_el_data_3::face, findpts_el_data_3::n, nr, ns, out, findpts_el_data_3::sides, tensor_mxm(), work, findpts_el_data_3::work, findpts_el_data_3::wtend, findpts_el_gface_3::x, and findpts_el_data_3::x.

Referenced by get_face().

Here is the call graph for this function:

Here is the caller graph for this function:

static void compute_face_data_st ( struct findpts_el_data_3 * fd )

static

Definition at line 310 of file findpts_el_3.c.

References i, findpts_el_data_3::n, nr, ns, out, findpts_el_data_3::sides, tensor_mtxm, work, findpts_el_data_3::work, findpts_el_data_3::wtend, and findpts_el_data_3::x.

Referenced by get_face().

Here is the caller graph for this function:

static void compute_face_data_tr ( struct findpts_el_data_3 * fd )

static

Definition at line 327 of file findpts_el_3.c.

References i, findpts_el_data_3::n, nr, ns, out, findpts_el_data_3::sides, tensor_mxm(), work, findpts_el_data_3::work, findpts_el_data_3::wtend, and findpts_el_data_3::x.

Referenced by get_face().

Here is the call graph for this function:

Here is the caller graph for this function:

static void compute_pt_data ( struct findpts_el_data_3 * fd )

static

Definition at line 414 of file findpts_el_3.c.

References findpts_el_gedge_3::d2xdn1, findpts_el_gedge_3::d2xdn2, findpts_el_gedge_3::dxdn1, findpts_el_gedge_3::dxdn2, get_edge(), findpts_el_gpt_3::hes, i, findpts_el_gpt_3::jac, findpts_el_data_3::n, nr, findpts_el_data_3::pt, tensor_mtxv(), work, findpts_el_data_3::work, findpts_el_data_3::wtend, x, findpts_el_gedge_3::x, and findpts_el_gpt_3::x.

Referenced by get_pt().

Here is the call graph for this function:

Here is the caller graph for this function:

static unsigned edge_index ( const unsigned x )

static

Definition at line 111 of file findpts_el_3.c.

References se, and y.

Referenced by findpt_edge(), and newton_vol().

Here is the caller graph for this function:

static unsigned face_index ( const unsigned x )

static

Definition at line 109 of file findpts_el_3.c.

References which_bit().

Referenced by findpt_face(), and newton_vol().

Here is the call graph for this function:

Here is the caller graph for this function:

static void findpt_edge	(	struct findpts_el_pt_3 *const	out,
		struct findpts_el_data_3 *const	fd,
		const struct findpts_el_pt_3 *const	p,
		const unsigned	pn,
		const double	tol
	)

static

Definition at line 1034 of file findpts_el_3.c.

References findpts_el_gedge_3::d2xdn1, findpts_el_gedge_3::d2xdn2, findpts_el_gedge_3::dxdn1, findpts_el_gedge_3::dxdn2, edge_index(), FLAG_MASK, findpts_el_pt_3::flags, get_edge(), i, findpts_el_data_3::lag, findpts_el_data_3::lag_data, n, findpts_el_data_3::n, newton_edge(), newton_face(), newton_vol(), plus_1_mod_3(), plus_2_mod_3(), findpts_el_pt_3::r, reject_prior_step_q(), tensor_dot, tensor_mtxv(), findpts_el_data_3::work, findpts_el_pt_3::x, and findpts_el_gedge_3::x.

Referenced by findpts_el_3().

Here is the call graph for this function:

Here is the caller graph for this function:

static void findpt_face	(	struct findpts_el_pt_3 *const	out,
		struct findpts_el_data_3 *const	fd,
		const struct findpts_el_pt_3 *const	p,
		const unsigned	pn,
		const double	tol
	)

static

Definition at line 958 of file findpts_el_3.c.

References findpts_el_gface_3::dxdn, face_index(), FLAG_MASK, findpts_el_pt_3::flags, get_face(), findpts_el_gpt_3::hes, i, findpts_el_data_3::lag, findpts_el_data_3::lag_data, findpts_el_data_3::n, newton_face(), newton_vol(), plus_1_mod_3(), plus_2_mod_3(), findpts_el_pt_3::r, reject_prior_step_q(), tensor_dot, tensor_mtxm, tensor_mxm(), findpts_el_data_3::work, findpts_el_pt_3::x, and findpts_el_gface_3::x.

Referenced by findpts_el_3().

Here is the call graph for this function:

Here is the caller graph for this function:

static void findpt_pt	(	struct findpts_el_pt_3 *const	out,
		struct findpts_el_data_3 *const	fd,
		const struct findpts_el_pt_3 *const	p,
		const unsigned	pn,
		const double	tol
	)

static

Definition at line 1115 of file findpts_el_3.c.

References CONVERGED_FLAG, findpts_el_pt_3::dist2p, FLAG_MASK, findpts_el_pt_3::flags, get_pt(), findpts_el_gpt_3::hes, i, findpts_el_gpt_3::jac, newton_edge(), newton_face(), newton_vol(), point_index(), findpts_el_pt_3::r, reject_prior_step_q(), x, findpts_el_pt_3::x, and findpts_el_gpt_3::x.

Referenced by findpts_el_3().

Here is the call graph for this function:

Here is the caller graph for this function:

static void findpt_vol	(	struct findpts_el_pt_3 *const	out,
		struct findpts_el_data_3 *const	fd,
		const struct findpts_el_pt_3 *const	p,
		const unsigned	pn,
		const double	tol
	)

static

Definition at line 921 of file findpts_el_3.c.

References i, findpts_el_data_3::lag, findpts_el_data_3::lag_data, findpts_el_data_3::n, newton_vol(), nr, ns, findpts_el_pt_3::r, reject_prior_step_q(), tensor_i2(), tensor_ig2(), tensor_mxm(), findpts_el_data_3::work, findpts_el_pt_3::x, and findpts_el_data_3::x.

Referenced by findpts_el_3().

Here is the call graph for this function:

Here is the caller graph for this function:

void findpts_el_3	(	struct findpts_el_data_3 *const	fd,
		const unsigned	npt,
		const double	tol
	)

Definition at line 1211 of file findpts_el_3.c.

References findpt_edge(), findpt_face(), findpt_pt(), findpt_vol(), FLAG_MASK, findpts_el_pt_3::flags, i, findpts_el_pt_3::index, n, findpts_el_data_3::n, num_constrained(), p, findpts_el_data_3::p, pt_flags_to_bin(), pt_flags_to_bin_noC(), findpts_el_pt_3::r, seed(), sum(), findpts_el_pt_3::x, and findpts_el_data_3::x.

Here is the call graph for this function:

void findpts_el_eval_3	(	double *const	out_base,
		const unsigned	out_stride,
		const double *const	r_base,
		const unsigned	r_stride,
		const unsigned	pn,
		const double *const	in,
		struct findpts_el_data_3 *const	fd
	)

Definition at line 1294 of file findpts_el_3.c.

References i, findpts_el_data_3::lag, findpts_el_data_3::lag_data, findpts_el_data_3::n, nr, ns, out, findpts_el_pt_3::r, tensor_i2(), tensor_mxm(), and findpts_el_data_3::work.

Here is the call graph for this function:

void findpts_el_free_3 ( struct findpts_el_data_3 *const fd )

Definition at line 283 of file findpts_el_3.c.

References findpts_el_data_3::p, and findpts_el_data_3::z.

void findpts_el_setup_3	(	struct findpts_el_data_3 *const	fd,
		const unsigned	n[3],
		const unsigned	npt_max
	)

Definition at line 202 of file findpts_el_3.c.

References findpts_el_gface_3::dxdn, findpts_el_data_3::face, gll_lag_setup, gll_lag_size, i, findpts_el_data_3::lag, findpts_el_data_3::lag_data, lobatto_nodes, n, findpts_el_data_3::n, findpts_el_data_3::npt_max, nr, ns, findpts_el_data_3::p, SET_EDGE, SET_FACE, findpts_el_data_3::side_init, findpts_el_data_3::sides, tmalloc, findpts_el_data_3::work, work_size(), findpts_el_data_3::wtend, findpts_el_gface_3::x, and findpts_el_data_3::z.

Here is the call graph for this function:

static const struct findpts_el_gedge_3* get_edge	(	struct findpts_el_data_3 *	fd,
		unsigned	ei
	)

static

Definition at line 404 of file findpts_el_3.c.

References compute_edge_data(), findpts_el_data_3::edge, and findpts_el_data_3::side_init.

Referenced by compute_pt_data(), and findpt_edge().

Here is the call graph for this function:

Here is the caller graph for this function:

static const struct findpts_el_gface_3* get_face	(	struct findpts_el_data_3 *	fd,
		unsigned	fi
	)

static

Definition at line 347 of file findpts_el_3.c.

References compute_face_data_rs(), compute_face_data_st(), compute_face_data_tr(), findpts_el_data_3::face, and findpts_el_data_3::side_init.

Referenced by compute_edge_data(), and findpt_face().

Here is the call graph for this function:

Here is the caller graph for this function:

static const struct findpts_el_gpt_3* get_pt	(	struct findpts_el_data_3 *	fd,
		unsigned	pi
	)

static

Definition at line 459 of file findpts_el_3.c.

References compute_pt_data(), findpts_el_data_3::pt, and findpts_el_data_3::side_init.

Referenced by findpt_pt().

Here is the call graph for this function:

Here is the caller graph for this function:

static void lin_solve_3	(	double	x[3],
		const double	A[9],
		const double	y[3]
	)

static

Definition at line 30 of file findpts_el_3.c.

Referenced by newton_vol().

Here is the caller graph for this function:

static void lin_solve_sym_2	(	double	x[2],
		const double	A[3],
		const double	y[2]
	)

static

Definition at line 51 of file findpts_el_3.c.

Referenced by newton_face(), and newton_vol().

Here is the caller graph for this function:

static void newton_edge	(	struct findpts_el_pt_3 *const	out,
		const double	jac[9],
		const double	rhes,
		const double	resid[3],
		const unsigned	de,
		const unsigned	dn1,
		const unsigned	dn2,
		unsigned	flags,
		const struct findpts_el_pt_3 *const	p,
		const double	tol
	)

static

Definition at line 857 of file findpts_el_3.c.

References A, CONVERGED_FLAG, findpts_el_pt_3::dist2p, EVAL, findpts_el_pt_3::flags, nr, findpts_el_pt_3::r, findpts_el_pt_3::tr, and y.

Referenced by findpt_edge(), and findpt_pt().

Here is the caller graph for this function:

static void newton_face	(	struct findpts_el_pt_3 *const	out,
		const double	jac[9],
		const double	rhes[3],
		const double	resid[3],
		const unsigned	d1,
		const unsigned	d2,
		const unsigned	dn,
		const unsigned	flags,
		const struct findpts_el_pt_3 *const	p,
		const double	tol
	)

static

Definition at line 745 of file findpts_el_3.c.

References A, CONVERGED_FLAG, findpts_el_pt_3::dist2p, EVAL, findpts_el_pt_3::flags, i, lin_solve_sym_2(), findpts_el_pt_3::r, SETR, findpts_el_pt_3::tr, and y.

Referenced by findpt_edge(), findpt_face(), and findpt_pt().

Here is the call graph for this function:

Here is the caller graph for this function:

static void newton_vol	(	struct findpts_el_pt_3 *const	out,
		const double	jac[9],
		const double	resid[3],
		const struct findpts_el_pt_3 *const	p,
		const double	tol
	)