nugrid.h File Reference

#include "local_definitions.h"

Go to the source code of this file.

Classes
struct	center_grid
struct	log_grid
struct	NUgrid

Enumerations
enum	grid_type { LINEAR , GENERAL , CENTER , LOG }

Functions
NUgrid *	create_center_grid (double start, double end, double ratio, int num_points)
NUgrid *	create_general_grid (double *points, int num_points)
NUgrid *	create_log_grid (double start, double end, int num_points)
void	destroy_grid (NUgrid *grid)

Enumeration Type Documentation

grid_type

enum grid_type

Enumerator
LINEAR
GENERAL
CENTER
LOG

Definition at line 26 of file nugrid.h.

{ LINEAR, GENERAL, CENTER, LOG } grid_type;

Function Documentation

create_center_grid()

NUgrid * create_center_grid	(	double	start,
		double	end,
		double	ratio,
		int	num_points )

Definition at line 83 of file nugrid.cpp.

{
  center_grid *grid = new center_grid;
  if (grid != NULL) {
    assert (ratio > 1.0);
    grid->start       = start;
    grid->end         = end;
    grid->center      = 0.5*(start + end);
    grid->num_points  = num_points;
    grid->half_points = num_points/2;
    grid->odd = ((num_points % 2) == 1);
    grid->b = log(ratio) / (double)(grid->half_points-1);
    grid->bInv = 1.0/grid->b;
    grid->points = new double[num_points];
    if (grid->odd) {
      grid->even_half = 0.0;  
      grid->odd_one   = 1;
      grid->a = 0.5*(end-start)/expm1(grid->b*grid->half_points);
      grid->aInv = 1.0/grid->a;
      for (int i=-grid->half_points; i<=grid->half_points; i++) {
    double sign;
    if (i<0) 
      sign = -1.0;
    else
      sign =  1.0;
    grid->points[i+grid->half_points] = 
      sign * grid->a*expm1(grid->b*abs(i))+grid->center;
      }
    }
    else {
      grid->even_half = 0.5;  
      grid->odd_one   = 0;
      grid->a = 
    0.5*(end-start)/expm1(grid->b*(-0.5+grid->half_points));
      grid->aInv = 1.0/grid->a;
      for (int i=-grid->half_points; i<grid->half_points; i++) {
    double sign;
    if (i<0) sign = -1.0; 
    else     sign =  1.0;
    grid->points[i+grid->half_points] = 
      sign * grid->a*expm1(grid->b*fabs(0.5+i)) + grid->center;
      }
    }
    grid->reverse_map = center_grid_reverse_map;
    grid->code = CENTER;
  }
  return (NUgrid*) grid;
}

create_general_grid()

NUgrid * create_general_grid	(	double *	points,
		int	num_points )

Definition at line 155 of file nugrid.cpp.

{
  NUgrid* grid = new NUgrid;
  if (grid != NULL) {
    grid->reverse_map = general_grid_reverse_map;
    grid->code = GENERAL;
    grid->points = new double[num_points];
    grid->start = points[0];
    grid->end   = points[num_points-1];
    grid->num_points = num_points;
    for (int i=0; i<num_points; i++) 
      grid->points[i] = points[i];
    grid->code = GENERAL;
  }
  return grid;
}

create_log_grid()

NUgrid * create_log_grid	(	double	start,
		double	end,
		int	num_points )

Definition at line 135 of file nugrid.cpp.

{
  log_grid *grid = new log_grid;
  grid->code = LOG;
  grid->start = start;
  grid->end = end;
  grid->num_points = num_points;
  grid->points = new double[num_points];
  grid->a = 1.0/(double)(num_points-1)*log(end/start);
  grid->ainv = 1.0/grid->a;
  grid->startinv = 1.0/start;
  for (int i=0; i<num_points; i++)
    grid->points[i] = start*exp(grid->a*(double)i);
  grid->reverse_map = log_grid_reverse_map;
  return (NUgrid*) grid;
}

References log_grid::a, log_grid::ainv, log_grid::code, log_grid::end, LOG, log_grid_reverse_map(), log_grid::num_points, log_grid::points, log_grid::reverse_map, log_grid::start, and log_grid::startinv.

destroy_grid()

void destroy_grid

(

NUgrid *

grid

)

Definition at line 173 of file nugrid.cpp.

{
  delete[] grid->points;
  delete grid;
}

References NUgrid::points.