/* store.c -- store molecule coefficients in Jacobian */
/* $Header$ */

#include "cctk.h"
#include "../InterpLocalUniform.h"

#include "structs.h"
#include "store.h"

/******************************************************************************/

/*
 * 1-D routines
 */

#undef  COEFF
#define COEFF(mi)	Jacobian_ptr[Jacobian_mi_stride*mi]

void LocalInterp_store_1d_cube2
      (fp factor, const struct coeffs_struct_1d_cube_size2 *coeffs,
       fp Jacobian_ptr[],
       int Jacobian_mi_stride)
{
#include "1d.cube.size2/store-coeffs.c"
}

void LocalInterp_store_1d_cube3
      (fp factor, const struct coeffs_struct_1d_cube_size3 *coeffs,
       fp Jacobian_ptr[],
       int Jacobian_mi_stride)
{
#include "1d.cube.size3/store-coeffs.c"
}

void LocalInterp_store_1d_cube4
      (fp factor, const struct coeffs_struct_1d_cube_size4 *coeffs,
       fp Jacobian_ptr[],
       int Jacobian_mi_stride)
{
#include "1d.cube.size4/store-coeffs.c"
}

void LocalInterp_store_1d_cube5
      (fp factor, const struct coeffs_struct_1d_cube_size5 *coeffs,
       fp Jacobian_ptr[],
       int Jacobian_mi_stride)
{
#include "1d.cube.size5/store-coeffs.c"
}

void LocalInterp_store_1d_cube6
      (fp factor, const struct coeffs_struct_1d_cube_size6 *coeffs,
       fp Jacobian_ptr[],
       int Jacobian_mi_stride)
{
#include "1d.cube.size6/store-coeffs.c"
}

void LocalInterp_store_1d_cube7
      (fp factor, const struct coeffs_struct_1d_cube_size7 *coeffs,
       fp Jacobian_ptr[],
       int Jacobian_mi_stride)
{
#include "1d.cube.size7/store-coeffs.c"
}

/******************************************************************************/

/*
 * 2-D routines
 */

#undef  COEFF
#define COEFF(mi,mj)	Jacobian_ptr[   Jacobian_mi_stride*mi \
				      + Jacobian_mj_stride*mj ]

void LocalInterp_store_2d_cube2
      (fp factor, const struct coeffs_struct_2d_cube_size2 *coeffs,
       fp Jacobian_ptr[],
       int Jacobian_mi_stride, int Jacobian_mj_stride)
{
#include "2d.cube.size2/store-coeffs.c"
}

void LocalInterp_store_2d_cube3
      (fp factor, const struct coeffs_struct_2d_cube_size3 *coeffs,
       fp Jacobian_ptr[],
       int Jacobian_mi_stride, int Jacobian_mj_stride)
{
#include "2d.cube.size3/store-coeffs.c"
}

void LocalInterp_store_2d_cube4
      (fp factor, const struct coeffs_struct_2d_cube_size4 *coeffs,
       fp Jacobian_ptr[],
       int Jacobian_mi_stride, int Jacobian_mj_stride)
{
#include "2d.cube.size4/store-coeffs.c"
}

void LocalInterp_store_2d_cube5
      (fp factor, const struct coeffs_struct_2d_cube_size5 *coeffs,
       fp Jacobian_ptr[],
       int Jacobian_mi_stride, int Jacobian_mj_stride)
{
#include "2d.cube.size5/store-coeffs.c"
}

void LocalInterp_store_2d_cube6
      (fp factor, const struct coeffs_struct_2d_cube_size6 *coeffs,
       fp Jacobian_ptr[],
       int Jacobian_mi_stride, int Jacobian_mj_stride)
{
#include "2d.cube.size6/store-coeffs.c"
}

/******************************************************************************/

/*
 * 3-D routines
 */

#undef  COEFF
#define COEFF(mi,mj,mk)	Jacobian_ptr[   Jacobian_mi_stride*mi \
				      + Jacobian_mj_stride*mj \
				      + Jacobian_mk_stride*mk ]

void LocalInterp_store_3d_cube2
      (fp factor, const struct coeffs_struct_3d_cube_size2 *coeffs,
       fp Jacobian_ptr[],
       int Jacobian_mi_stride, int Jacobian_mj_stride, int Jacobian_mk_stride)
{
#include "3d.cube.size2/store-coeffs.c"
}

void LocalInterp_store_3d_cube3
      (fp factor, const struct coeffs_struct_3d_cube_size3 *coeffs,
       fp Jacobian_ptr[],
       int Jacobian_mi_stride, int Jacobian_mj_stride, int Jacobian_mk_stride)
{
#include "3d.cube.size3/store-coeffs.c"
}

void LocalInterp_store_3d_cube4
      (fp factor, const struct coeffs_struct_3d_cube_size4 *coeffs,
       fp Jacobian_ptr[],
       int Jacobian_mi_stride, int Jacobian_mj_stride, int Jacobian_mk_stride)
{
#include "3d.cube.size4/store-coeffs.c"
}

void LocalInterp_store_3d_cube5
      (fp factor, const struct coeffs_struct_3d_cube_size5 *coeffs,
       fp Jacobian_ptr[],
       int Jacobian_mi_stride, int Jacobian_mj_stride, int Jacobian_mk_stride)
{
#include "3d.cube.size5/store-coeffs.c"
}

void LocalInterp_store_3d_cube6
      (fp factor, const struct coeffs_struct_3d_cube_size6 *coeffs,
       fp Jacobian_ptr[],
       int Jacobian_mi_stride, int Jacobian_mj_stride, int Jacobian_mk_stride)
{
#include "3d.cube.size6/store-coeffs.c"
}