#include "moab/FindPtFuncs.h"

Include dependency graph for tensor.c:

Functions
static void	mxm_cc (const realType A, unsigned na, const realType B, unsigned nb, realType *C, unsigned nc)

static void	mxm_rc (const realType A, unsigned na, const realType B, unsigned nb, realType *C, unsigned nc)

static void	mxm_cr (const realType A, unsigned na, const realType B, unsigned nb, realType *C, unsigned nc)

static void	mxv_c (realType y, unsigned ny, const realType A, const realType *x, unsigned nx)

static void	mxv_r (realType y, unsigned ny, const realType A, const realType *x, unsigned nx)

static realType	inner (const realType u, const realType v, unsigned n)

void	tensor_c1 (const realType R, unsigned mr, unsigned nr, const realType u, realType *v)

void	tensor_r1 (const realType R, unsigned mr, unsigned nr, const realType u, realType *v)

void	tensor_c2 (const realType R, unsigned mr, unsigned nr, const realType S, unsigned ms, unsigned ns, const realType u, realType v, realType *W)

void	tensor_r2 (const realType R, unsigned mr, unsigned nr, const realType S, unsigned ms, unsigned ns, const realType u, realType v, realType *W)

void	tensor_c3 (const realType R, unsigned mr, unsigned nr, const realType S, unsigned ms, unsigned ns, const realType T, unsigned mt, unsigned nt, const realType u, realType v, realType W, realType *Z)

void	tensor_r3 (const realType R, unsigned mr, unsigned nr, const realType S, unsigned ms, unsigned ns, const realType T, unsigned mt, unsigned nt, const realType u, realType v, realType W, realType *Z)

realType	tensor_i1 (const realType Jr, unsigned nr, const realType u)

realType	tensor_i2 (const realType Jr, unsigned nr, const realType Js, unsigned ns, const realType u, realType work)

realType	tensor_i3 (const realType Jr, unsigned nr, const realType Js, unsigned ns, const realType Jt, unsigned nt, const realType u, realType *work)

realType	tensor_ig1 (const realType Jr, const realType Dr, unsigned nr, const realType u, realType g)

realType	tensor_ig2 (const realType Jr, const realType Dr, unsigned nr, const realType Js, const realType Ds, unsigned ns, const realType u, realType g, realType *work)

realType	tensor_ig3 (const realType Jr, const realType Dr, unsigned nr, const realType Js, const realType Ds, unsigned ns, const realType Jt, const realType Dt, unsigned nt, const realType u, realType g, realType *work)

Function Documentation

◆ inner()

static realType inner	(	const realType *	u,
		const realType *	v,
		unsigned	n
	)

static

Definition at line 132 of file tensor.c.

 {
     const realType* u_end = u + n;
     realType sum          = *u++ * *v++;
     while( u != u_end )
     {
         sum += *u++ * *v++;
     }
     return sum;
 }

References moab::sum().

Referenced by run_quality_optimizer(), tensor_i1(), tensor_i2(), tensor_i3(), tensor_ig1(), tensor_ig2(), and tensor_ig3().

◆ mxm_cc()

static void mxm_cc	(	const realType *	A,
		unsigned	na,
		const realType *	B,
		unsigned	nb,
		realType *	C,
		unsigned	nc
	)

static

Definition at line 13 of file tensor.c.

 {
     unsigned i, j, k;
     realType* Ccol       = C;
     const realType* Bcol = B;
     for( j = 0; j < nc; ++j, Ccol += na, Bcol += nb )
     {
         const realType* Acol = A;
         for( i = 0; i < na; ++i )
             Ccol[i] = 0;
         for( k = 0; k < nb; ++k, Acol += na )
             for( i = 0; i < na; ++i )
                 Ccol[i] += Acol[i] * Bcol[k];
     }
 }

Referenced by tensor_c2(), tensor_c3(), tensor_r2(), and tensor_r3().

◆ mxm_cr()

static void mxm_cr	(	const realType *	A,
		unsigned	na,
		const realType *	B,
		unsigned	nb,
		realType *	C,
		unsigned	nc
	)

static

Definition at line 46 of file tensor.c.

 {
     unsigned i, j, k;
     const realType *Acol = A, *Brow = B;
     for( i = 0; i < na * nc; ++i )
         C[i] = 0;
     for( k = 0; k < nb; ++k, Acol += na, Brow += nc )
     {
         realType* Ccol = C;
         for( j = 0; j < nc; ++j, Ccol += na )
             for( i = 0; i < na; ++i )
                 Ccol[i] += Acol[i] * Brow[j];
     }
 }

Referenced by tensor_c2(), and tensor_c3().

◆ mxm_rc()

static void mxm_rc	(	const realType *	A,
		unsigned	na,
		const realType *	B,
		unsigned	nb,
		realType *	C,
		unsigned	nc
	)

static

Definition at line 29 of file tensor.c.

 {
     unsigned i, j, k;
     realType* Ccol       = C;
     const realType* Bcol = B;
     for( j = 0; j < nc; ++j, Ccol += na, Bcol += nb )
     {
         const realType* Arow = A;
         for( i = 0; i < na; ++i, Arow += nb )
         {
             Ccol[i] = 0;
             for( k = 0; k < nb; ++k )
                 Ccol[i] += Arow[k] * Bcol[k];
         }
     }
 }

Referenced by tensor_r2(), and tensor_r3().

◆ mxv_c()

static void mxv_c	(	realType *	y,
		unsigned	ny,
		const realType *	A,
		const realType *	x,
		unsigned	nx
	)

static

Definition at line 89 of file tensor.c.

 {
     realType *yp = y, *y_end = y + ny;
     const realType* x_end = x + nx;
     realType xk           = *x;
     do
     {
         *yp++ = *A++ * xk;
     } while( yp != y_end );
     for( ++x; x != x_end; ++x )
     {
         xk = *x;
         yp = y;
         do
         {
             *yp++ += *A++ * xk;
         } while( yp != y_end );
     }
 }

Referenced by tensor_c1().

◆ mxv_r()

static void mxv_r	(	realType *	y,
		unsigned	ny,
		const realType *	A,
		const realType *	x,
		unsigned	nx
	)

static

Definition at line 109 of file tensor.c.

 {
     realType* y_end       = y + ny;
     const realType* x_end = x + nx;
     do
     {
         const realType* xp = x;
         realType sum       = *A++ * *xp++;
         while( xp != x_end )
         {
             sum += *A++ * *xp++;
         }
         *y++ = sum;
     } while( y != y_end );
 }

References moab::sum().

Referenced by tensor_i2(), tensor_i3(), tensor_ig2(), tensor_ig3(), and tensor_r1().

◆ tensor_c1()

void tensor_c1	(	const realType *	R,
		unsigned	mr,
		unsigned	nr,
		const realType *	u,
		realType *	v
	)

Definition at line 155 of file tensor.c.

 {
     mxv_c( v, mr, R, u, nr );
 }

References mxv_c(), nr, and moab::R.

◆ tensor_c2()

void tensor_c2	(	const realType *	R,
		unsigned	mr,
		unsigned	nr,
		const realType *	S,
		unsigned	ms,
		unsigned	ns,
		const realType *	u,
		realType *	v,
		realType *	W
	)

Definition at line 166 of file tensor.c.

 {
     mxm_cc( R, mr, u, nr, W, ns );
     mxm_cr( W, mr, S, ns, v, ms );
 }

References mxm_cc(), mxm_cr(), nr, and moab::R.

◆ tensor_c3()

void tensor_c3	(	const realType *	R,
		unsigned	mr,
		unsigned	nr,
		const realType *	S,
		unsigned	ms,
		unsigned	ns,
		const realType *	T,
		unsigned	mt,
		unsigned	nt,
		const realType *	u,
		realType *	v,
		realType *	W,
		realType *	Z
	)

Definition at line 197 of file tensor.c.

 {
     unsigned n, mrns = mr * ns, mrms = mr * ms;
     realType* Zp = Z;
     mxm_cc( R, mr, u, nr, W, ns * nt );
     for( n = 0; n < nt; ++n, W += mrns, Zp += mrms )
         mxm_cr( W, mr, S, ns, Zp, ms );
     mxm_cr( Z, mrms, T, nt, v, mt );
 }

References mxm_cc(), mxm_cr(), nr, and moab::R.

◆ tensor_i1()

realType tensor_i1	(	const realType *	Jr,
		unsigned	nr,
		const realType *	u
	)

Definition at line 255 of file tensor.c.

 {
     return inner( Jr, u, nr );
 }

References inner(), and nr.

Referenced by opt_edge_hess_2(), opt_edge_hess_3(), opt_edge_intp_2(), and opt_edge_intp_3().

◆ tensor_i2()

realType tensor_i2	(	const realType *	Jr,
		unsigned	nr,
		const realType *	Js,
		unsigned	ns,
		const realType *	u,
		realType *	work
	)

Definition at line 261 of file tensor.c.

 {
     mxv_r( work, ns, u, Jr, nr );
     return inner( Js, work, ns );
 }

References inner(), mxv_r(), and nr.

Referenced by moab::SpectralQuad::evalFcn(), moab::Element::SpectralQuad::evaluate(), moab::Element::SpectralQuad::evaluate_scalar_field(), findpt_eval_2(), opt_face_hess_3(), and opt_face_intp_3().

◆ tensor_i3()

realType tensor_i3	(	const realType *	Jr,
		unsigned	nr,
		const realType *	Js,
		unsigned	ns,
		const realType *	Jt,
		unsigned	nt,
		const realType *	u,
		realType *	work
	)

Definition at line 273 of file tensor.c.

 {
     realType* work2 = work + nt;
     mxv_r( work2, ns * nt, u, Jr, nr );
     mxv_r( work, nt, work2, Js, ns );
     return inner( Jt, work, nt );
 }

References inner(), mxv_r(), and nr.

Referenced by moab::Element::SpectralHex::evaluate(), moab::element_utility::Spectral_hex_map< _Matrix >::evaluate(), moab::Element::SpectralHex::evaluate_scalar_field(), moab::element_utility::Spectral_hex_map< _Matrix >::evaluate_scalar_field(), findpt_eval_3(), and moab::ElemUtil::hex_eval().

◆ tensor_ig1()

realType tensor_ig1	(	const realType *	Jr,
		const realType *	Dr,
		unsigned	nr,
		const realType *	u,
		realType *	g
	)

Definition at line 304 of file tensor.c.

 {
     *g = inner( Dr, u, nr );
     return inner( Jr, u, nr );
 }

References inner(), and nr.

Referenced by opt_edge_intp_2(), and opt_edge_intp_3().

◆ tensor_ig2()

realType tensor_ig2	(	const realType *	Jr,
		const realType *	Dr,
		unsigned	nr,
		const realType *	Js,
		const realType *	Ds,
		unsigned	ns,
		const realType *	u,
		realType *	g,
		realType *	work
	)

Definition at line 311 of file tensor.c.

 {
     realType *a = work, *ar = a + ns;
     mxv_r( a, ns, u, Jr, nr );
     mxv_r( ar, ns, u, Dr, nr );
     g[0] = inner( Js, ar, ns );
     g[1] = inner( Ds, a, ns );
     return inner( Js, a, ns );
 }

References inner(), mxv_r(), and nr.

Referenced by obbox_calc_2(), opt_area_intp_2(), opt_face_hess_3(), and opt_face_intp_3().

◆ tensor_ig3()

realType tensor_ig3	(	const realType *	Jr,
		const realType *	Dr,
		unsigned	nr,
		const realType *	Js,
		const realType *	Ds,
		unsigned	ns,
		const realType *	Jt,
		const realType *	Dt,
		unsigned	nt,
		const realType *	u,
		realType *	g,
		realType *	work
	)

Definition at line 330 of file tensor.c.

 {
     unsigned nsnt = ns * nt;
     realType *a = work, *ar = a + nsnt, *b = ar + nsnt, *br = b + ns, *bs = br + ns;
     mxv_r( a, nsnt, u, Jr, nr );
     mxv_r( ar, nsnt, u, Dr, nr );
     mxv_r( b, nt, a, Js, ns );
     mxv_r( br, nt, ar, Js, ns );
     mxv_r( bs, nt, a, Ds, ns );
     g[0] = inner( Jt, br, nt );
     g[1] = inner( Jt, bs, nt );
     g[2] = inner( Dt, b, nt );
     return inner( Jt, b, nt );
 }

References inner(), mxv_r(), and nr.

Referenced by obbox_calc_3(), and opt_vol_intp_3().

◆ tensor_r1()

void tensor_r1	(	const realType *	R,
		unsigned	mr,
		unsigned	nr,
		const realType *	u,
		realType *	v
	)

Definition at line 160 of file tensor.c.

 {
     mxv_r( v, mr, R, u, nr );
 }

References mxv_r(), nr, and moab::R.

◆ tensor_r2()

void tensor_r2	(	const realType *	R,
		unsigned	mr,
		unsigned	nr,
		const realType *	S,
		unsigned	ms,
		unsigned	ns,
		const realType *	u,
		realType *	v,
		realType *	W
	)

Definition at line 181 of file tensor.c.

 {
     mxm_rc( R, mr, u, nr, W, ns );
     mxm_cc( W, mr, S, ns, v, ms );
 }

References mxm_cc(), mxm_rc(), nr, and moab::R.

◆ tensor_r3()

void tensor_r3	(	const realType *	R,
		unsigned	mr,
		unsigned	nr,
		const realType *	S,
		unsigned	ms,
		unsigned	ns,
		const realType *	T,
		unsigned	mt,
		unsigned	nt,
		const realType *	u,
		realType *	v,
		realType *	W,
		realType *	Z
	)

Definition at line 221 of file tensor.c.

 {
     unsigned n, mrns = mr * ns, mrms = mr * ms;
     realType* Zp = Z;
     mxm_rc( R, mr, u, nr, W, ns * nt );
     for( n = 0; n < nt; ++n, W += mrns, Zp += mrms )
         mxm_cc( W, mr, S, ns, Zp, ms );
     mxm_cc( Z, mrms, T, nt, v, mt );
 }

References mxm_cc(), mxm_rc(), nr, and moab::R.

Functions

Function Documentation

◆ inner()

◆ mxm_cc()

◆ mxm_cr()

◆ mxm_rc()

◆ mxv_c()

◆ mxv_r()

◆ tensor_c1()

◆ tensor_c2()

◆ tensor_c3()

◆ tensor_i1()

◆ tensor_i2()

◆ tensor_i3()

◆ tensor_ig1()

◆ tensor_ig2()

◆ tensor_ig3()

◆ tensor_r1()

◆ tensor_r2()

◆ tensor_r3()