triangle_felippa_rule.c File Reference

#include <stdlib.h>
#include <stdio.h>
#include <time.h>
#include <std::string.h>
#include <math.h>
#include "triangle_felippa_rule.h"

Go to the source code of this file.

Functions
void	comp_next (int n, int k, int a[], int *more, int *h, int *t)
double *	monomial_value (int m, int n, int e[], double x[])
void	r8vec_copy (int n, double a1[], double a2[])
double	r8vec_dot_product (int n, double a1[], double a2[])
void	subcomp_next (int n, int k, int a[], int *more, int *h, int *t)
double	triangle_unit_monomial (int expon[2])
void	triangle_unit_o01 (double w[], double xy[])
void	triangle_unit_o03 (double w[], double xy[])
void	triangle_unit_o03b (double w[], double xy[])
void	triangle_unit_o06 (double w[], double xy[])
void	triangle_unit_o06b (double w[], double xy[])
void	triangle_unit_o07 (double w[], double xy[])
void	triangle_unit_o12 (double w[], double xy[])
double	triangle_unit_volume ()

Function Documentation

comp_next()

void comp_next	(	int	n,
		int	k,
		int	a[],
		int *	more,
		int *	h,
		int *	t )

Definition at line 11 of file triangle_felippa_rule.c.

{
  int i;
 
  if ( !( *more ) )
  {
    *t = n;
    *h = 0;
    a[0] = n;
    for ( i = 1; i < k; i++ )
    {
       a[i] = 0;
    }
  }
  else
  {
    if ( 1 < *t )
    {
      *h = 0;
    }
    *h = *h + 1;
    *t = a[*h-1];
    a[*h-1] = 0;
    a[0] = *t - 1;
    a[*h] = a[*h] + 1;
  }
 
  *more = ( a[k-1] != n );
 
  return;
}

monomial_value()

double * monomial_value	(	int	m,
		int	n,
		int	e[],
		double	x[] )

Definition at line 153 of file triangle_felippa_rule.c.

{
  int i;
  int j;
  double *v;
 
  v = ( double * ) malloc ( n * sizeof ( double ) );
 
  for ( j = 0; j < n; j++ )
  {
    v[j] = 1.0;
  }
 
  for ( i = 0; i < m; i++ )
  {
    if ( 0 != e[i] )
    {
      for ( j = 0; j < n; j++ )
      {
        v[j] = v[j] * pow ( x[i+j*m], e[i] );
      }
    }
  }
 
  return v;
}

r8vec_copy()

void r8vec_copy	(	int	n,
		double	a1[],
		double	a2[] )

Definition at line 220 of file triangle_felippa_rule.c.

{
  int i;
 
  for ( i = 0; i < n; i++ )
  {
    a2[i] = a1[i];
  }
  return;
}

r8vec_dot_product()

double r8vec_dot_product	(	int	n,
		double	a1[],
		double	a2[] )

Definition at line 263 of file triangle_felippa_rule.c.

{
  int i;
  double value;
 
  value = 0.0;
  for ( i = 0; i < n; i++ )
  {
    value = value + a1[i] * a2[i];
  }
  return value;
}

subcomp_next()

void subcomp_next	(	int	n,
		int	k,
		int	a[],
		int *	more,
		int *	h,
		int *	t )

Definition at line 304 of file triangle_felippa_rule.c.

{
  int i;
  static int more2 = 0;
  static int n2 = 0;
/*
  The first computation.
*/
  if ( !( *more ) )
  {
    n2 = 0;
 
    for ( i = 0; i < k; i++ )
    {
      a[i] = 0;
    }
    more2 = 0;
    *h = 0;
    *t = 0;
 
    *more = 1;
  }
/*
  Do the next element at the current value of N.
*/
  else if ( more2 )
  {
    comp_next ( n2, k, a, &more2, h, t );
  }
  else
  {
    more2 = 0;
    n2 = n2 + 1;
 
    comp_next ( n2, k, a, &more2, h, t );
  }
/*
  Termination occurs if MORE2 = FALSE and N2 = N.
*/
  if ( !more2 && n2 == n )
  {
    *more = 0;
  }
 
  return;
}

triangle_unit_monomial()

double triangle_unit_monomial

(

int

expon[2]

)

Definition at line 453 of file triangle_felippa_rule.c.

{
  int i;
  int k;
  double value;
/*
  The first computation ends with VALUE = 1.0;
*/
  value = 1.0;
 
/* k = 0;
 
  The first loop simply computes 1 so we short circuit it!
 
 for ( i = 1; i <= expon[0]; i++ )
 {
   k = k + 1;
   value = value * ( double ) ( i ) / ( double ) ( k );
 }*/
 
  k = expon[0];
 
  for ( i = 1; i <= expon[1]; i++ )
  {
    k = k + 1;
    value = value * ( double ) ( i ) / ( double ) ( k );
  }
 
  k = k + 1;
  value = value / ( double ) ( k );
 
  k = k + 1;
  value = value / ( double ) ( k );
 
  return value;
}

triangle_unit_o01()

void triangle_unit_o01	(	double	w[],
		double	xy[] )

Definition at line 534 of file triangle_felippa_rule.c.

{
  int order = 1;
 
  double w_save[1] = {
    1.0 };
  double xy_save[2*1] = {
    0.33333333333333333333,  0.33333333333333333333 };
 
  r8vec_copy ( order, w_save, w );
  r8vec_copy ( 2*order, xy_save, xy );
 
  return;
}

triangle_unit_o03()

void triangle_unit_o03	(	double	w[],
		double	xy[] )

Definition at line 592 of file triangle_felippa_rule.c.

{
  int order = 3;
 
  double w_save[3] = {
    0.33333333333333333333,
    0.33333333333333333333,
    0.33333333333333333333 };
  double xy_save[2*3] = {
    0.66666666666666666667,  0.16666666666666666667,
    0.16666666666666666667,  0.66666666666666666667,
    0.16666666666666666667,  0.16666666666666666667 };
 
  r8vec_copy ( order, w_save, w );
  r8vec_copy ( 2*order, xy_save, xy );
 
  return;
}

triangle_unit_o03b()

void triangle_unit_o03b	(	double	w[],
		double	xy[] )

Definition at line 654 of file triangle_felippa_rule.c.

{
  int order = 3;
 
  double w_save[3] = {
    0.33333333333333333333,
    0.33333333333333333333,
    0.33333333333333333333 };
  double xy_save[2*3] = {
    0.0,  0.5,
    0.5,  0.0,
    0.5,  0.5 };
 
  r8vec_copy ( order, w_save, w );
  r8vec_copy ( 2*order, xy_save, xy );
 
  return;
}

triangle_unit_o06()

void triangle_unit_o06	(	double	w[],
		double	xy[] )

Definition at line 716 of file triangle_felippa_rule.c.

{
  int order = 6;
 
  double w_save[6] = {
    0.22338158967801146570,
    0.22338158967801146570,
    0.22338158967801146570,
    0.10995174365532186764,
    0.10995174365532186764,
    0.10995174365532186764 };
  double xy_save[2*6] = {
    0.10810301816807022736,  0.44594849091596488632,
    0.44594849091596488632,  0.10810301816807022736,
    0.44594849091596488632,  0.44594849091596488632,
    0.81684757298045851308,  0.091576213509770743460,
    0.091576213509770743460,  0.81684757298045851308,
    0.091576213509770743460,  0.091576213509770743460 };
 
  r8vec_copy ( order, w_save, w );
  r8vec_copy ( 2*order, xy_save, xy );
 
  return;
}

triangle_unit_o06b()

void triangle_unit_o06b	(	double	w[],
		double	xy[] )

Definition at line 784 of file triangle_felippa_rule.c.

{
  int order = 6;
 
  double w_save[6] = {
    0.30000000000000000000,
    0.30000000000000000000,
    0.30000000000000000000,
    0.033333333333333333333,
    0.033333333333333333333,
    0.033333333333333333333 };
  double xy_save[2*6] = {
    0.66666666666666666667,  0.16666666666666666667,
    0.16666666666666666667,  0.66666666666666666667,
    0.16666666666666666667,  0.16666666666666666667,
    0.0,  0.5,
    0.5,  0.0,
    0.5,  0.5 };
 
  r8vec_copy ( order, w_save, w );
  r8vec_copy ( 2*order, xy_save, xy );
 
  return;
}

triangle_unit_o07()

void triangle_unit_o07	(	double	w[],
		double	xy[] )

Definition at line 852 of file triangle_felippa_rule.c.

{
  int order = 7;
 
  double w_save[7] = {
    0.12593918054482715260,
    0.12593918054482715260,
    0.12593918054482715260,
    0.13239415278850618074,
    0.13239415278850618074,
    0.13239415278850618074,
    0.22500000000000000000 };
  double xy_save[2*7] = {
    0.79742698535308732240,  0.10128650732345633880,
    0.10128650732345633880,  0.79742698535308732240,
    0.10128650732345633880,  0.10128650732345633880,
    0.059715871789769820459,  0.47014206410511508977,
    0.47014206410511508977,  0.059715871789769820459,
    0.47014206410511508977,  0.47014206410511508977,
    0.33333333333333333333,  0.33333333333333333333 };
 
  r8vec_copy ( order, w_save, w );
  r8vec_copy ( 2*order, xy_save, xy );
 
  return;
}

triangle_unit_o12()

void triangle_unit_o12	(	double	w[],
		double	xy[] )

Definition at line 922 of file triangle_felippa_rule.c.

{
  int order = 12;
 
  double w_save[12] = {
     0.050844906370206816921,
     0.050844906370206816921,
     0.050844906370206816921,
     0.11678627572637936603,
     0.11678627572637936603,
     0.11678627572637936603,
     0.082851075618373575194,
     0.082851075618373575194,
     0.082851075618373575194,
     0.082851075618373575194,
     0.082851075618373575194,
     0.082851075618373575194 };
  double xy_save[2*12] = {
    0.87382197101699554332,  0.063089014491502228340,
    0.063089014491502228340,  0.87382197101699554332,
    0.063089014491502228340,  0.063089014491502228340,
    0.50142650965817915742,  0.24928674517091042129,
    0.24928674517091042129,  0.50142650965817915742,
    0.24928674517091042129,  0.24928674517091042129,
    0.053145049844816947353,  0.31035245103378440542,
    0.31035245103378440542,  0.053145049844816947353,
    0.053145049844816947353,  0.63650249912139864723,
    0.31035245103378440542,  0.63650249912139864723,
    0.63650249912139864723,  0.053145049844816947353,
    0.63650249912139864723,  0.31035245103378440542 };
 
  r8vec_copy ( order, w_save, w );
  r8vec_copy ( 2*order, xy_save, xy );
 
  return;
}

triangle_unit_volume()

double triangle_unit_volume

(

)

Definition at line 1002 of file triangle_felippa_rule.c.

{
  double volume;
 
  volume = 1.0 / 2.0;
 
  return volume;
}