CompactSchemeInitialize.c File Reference

Initializes the compact schemes. More...

#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <interpolation.h>
#include <mpivars.h>
#include <hypar.h>

Go to the source code of this file.

Functions
int	CompactSchemeInitialize (void *s, void *m, char *type)

Detailed Description

Initializes the compact schemes.

Author

Debojyoti Ghosh

Definition in file CompactSchemeInitialize.c.

Function Documentation

CompactSchemeInitialize()

int CompactSchemeInitialize	(	void *	s,
		void *	m,
		char *	type
	)

This function initializes the compact finite-difference methods: allocates the arrays to store the tridiagonal system.

Parameters

s	Solver object of type HyPar
m	MPI object of type MPIVariables
type	Type of interpolation

Definition at line 17 of file CompactSchemeInitialize.c.

{
  HyPar         *solver = (HyPar*) s;
  MPIVariables  *mpi    = (MPIVariables*) m;
  CompactScheme *compact   = (CompactScheme*) solver->compact;

  int nvars = solver->nvars;
  int ndims = solver->ndims;

  int size = 1, d;
  for (d=0; d<solver->ndims; d++) size *= (solver->dim_local[d]+1);
  size *= solver->nvars;
  if (!strcmp(solver->interp_type,_CHARACTERISTIC_)) size *= solver->nvars;

  compact->A = (double*) calloc (size, sizeof(double));
  compact->B = (double*) calloc (size, sizeof(double));
  compact->C = (double*) calloc (size, sizeof(double));
  compact->R = (double*) calloc (size, sizeof(double));

  compact->sendbuf = (double*) calloc (size, sizeof(double));
  compact->recvbuf = (double*) calloc (size, sizeof(double));

  return(0);
}