330e465d87
versions of StructGrid and StructStencil, and follows drift in naming conventions, calling conventions, etc., since original CI_struct code was written.
226 lines
7.4 KiB
C
226 lines
7.4 KiB
C
/*BHEADER**********************************************************************
|
|
* (c) 1997 The Regents of the University of California
|
|
*
|
|
* See the file COPYRIGHT_and_DISCLAIMER for a complete copyright
|
|
* notice, contact person, and disclaimer.
|
|
*
|
|
* $Revision$
|
|
*********************************************************************EHEADER*/
|
|
/******************************************************************************
|
|
*
|
|
* HYPRE_StructInterfaceMatrix interface
|
|
*
|
|
*****************************************************************************/
|
|
|
|
#include "headers.h"
|
|
|
|
/*--------------------------------------------------------------------------
|
|
* HYPRE_StructInterfaceMatrixCreate
|
|
*--------------------------------------------------------------------------*/
|
|
|
|
int
|
|
HYPRE_StructInterfaceMatrixCreate( MPI_Comm context,
|
|
HYPRE_StructGrid grid,
|
|
HYPRE_StructStencil stencil,
|
|
HYPRE_StructInterfaceMatrix *matrix )
|
|
{
|
|
*matrix = (HYPRE_StructInterfaceMatrix)
|
|
hypre_NewStructInterfaceMatrix( context,
|
|
(hypre_StructGrid *) grid,
|
|
(hypre_StructStencil *) stencil );
|
|
return 0;
|
|
}
|
|
|
|
/*--------------------------------------------------------------------------
|
|
* HYPRE_StructInterfaceMatrixDestroy
|
|
*--------------------------------------------------------------------------*/
|
|
|
|
int
|
|
HYPRE_StructInterfaceMatrixDestroy( HYPRE_StructInterfaceMatrix matrix )
|
|
{
|
|
return( hypre_FreeStructInterfaceMatrix( (hypre_StructInterfaceMatrix *) matrix ) );
|
|
}
|
|
|
|
/*--------------------------------------------------------------------------
|
|
* HYPRE_SetStructInterfaceMatrixCoeffs
|
|
*--------------------------------------------------------------------------*/
|
|
|
|
int
|
|
HYPRE_SetStructInterfaceMatrixCoeffs( HYPRE_StructInterfaceMatrix matrix,
|
|
int *grid_index,
|
|
double *coeffs )
|
|
{
|
|
hypre_StructInterfaceMatrix *new_matrix;
|
|
hypre_Index new_grid_index;
|
|
|
|
int d;
|
|
int ierr;
|
|
|
|
new_matrix = (hypre_StructInterfaceMatrix *) matrix;
|
|
for (d = 0;
|
|
d < hypre_StructGridDim(hypre_StructInterfaceMatrixStructGrid(new_matrix));
|
|
d++)
|
|
{
|
|
hypre_IndexD(new_grid_index, d) = grid_index[d];
|
|
}
|
|
|
|
ierr = hypre_SetStructInterfaceMatrixCoeffs( new_matrix, new_grid_index, coeffs );
|
|
|
|
return (ierr);
|
|
}
|
|
|
|
/*--------------------------------------------------------------------------
|
|
* HYPRE_StructInterfaceMatrixSetValues
|
|
*--------------------------------------------------------------------------*/
|
|
|
|
/*
|
|
Not currently supported.
|
|
*/
|
|
|
|
int
|
|
HYPRE_StructInterfaceMatrixSetValues( HYPRE_StructInterfaceMatrix matrix,
|
|
int *grid_index,
|
|
int num_stencil_indices,
|
|
int *stencil_indices,
|
|
double *coeffs )
|
|
{
|
|
return( 0 );
|
|
}
|
|
|
|
/*--------------------------------------------------------------------------
|
|
* HYPRE_StructInterfaceMatrixSetBoxValues
|
|
*--------------------------------------------------------------------------*/
|
|
|
|
/*
|
|
Don't completely understand this function. Why is new_matrix needed?
|
|
*/
|
|
|
|
int
|
|
HYPRE_StructInterfaceMatrixSetBoxValues( HYPRE_StructInterfaceMatrix matrix,
|
|
int *lower_grid_index,
|
|
int *upper_grid_index,
|
|
int num_stencil_indices,
|
|
int *stencil_indices,
|
|
double *coeffs )
|
|
|
|
{
|
|
hypre_StructInterfaceMatrix *new_matrix;
|
|
hypre_Index new_lower_grid_index;
|
|
hypre_Index new_upper_grid_index;
|
|
|
|
int d;
|
|
int ierr;
|
|
|
|
new_matrix = (hypre_StructInterfaceMatrix *) matrix;
|
|
for (d = 0;
|
|
d < hypre_StructGridDim(hypre_StructInterfaceMatrixStructGrid(new_matrix));
|
|
d++)
|
|
{
|
|
hypre_IndexD(new_lower_grid_index, d) = lower_grid_index[d];
|
|
hypre_IndexD(new_upper_grid_index, d) = upper_grid_index[d];
|
|
}
|
|
|
|
ierr = hypre_SetStructInterfaceMatrixBoxValues( (hypre_StructInterfaceMatrix *) matrix ,
|
|
new_lower_grid_index,
|
|
new_upper_grid_index,
|
|
num_stencil_indices,
|
|
stencil_indices,
|
|
coeffs );
|
|
|
|
|
|
return (ierr);
|
|
}
|
|
|
|
/*--------------------------------------------------------------------------
|
|
* HYPRE_StructInterfaceMatrixInitialize
|
|
*--------------------------------------------------------------------------*/
|
|
|
|
int
|
|
HYPRE_StructInterfaceMatrixInitialize( HYPRE_StructInterfaceMatrix matrix )
|
|
/* No-op */
|
|
{
|
|
return( 0 );
|
|
}
|
|
|
|
/*--------------------------------------------------------------------------
|
|
* HYPRE_StructInterfaceMatrixAssemble
|
|
*--------------------------------------------------------------------------*/
|
|
|
|
int
|
|
HYPRE_StructInterfaceMatrixAssemble( HYPRE_StructInterfaceMatrix matrix )
|
|
{
|
|
return( hypre_AssembleStructInterfaceMatrix( (hypre_StructInterfaceMatrix *) matrix ) );
|
|
}
|
|
|
|
/*--------------------------------------------------------------------------
|
|
* HYPRE_StructInterfaceMatrixGetData
|
|
*--------------------------------------------------------------------------*/
|
|
|
|
void *
|
|
HYPRE_StructInterfaceMatrixGetData( HYPRE_StructInterfaceMatrix matrix )
|
|
{
|
|
return( hypre_StructInterfaceMatrixData( (hypre_StructInterfaceMatrix *) matrix ) );
|
|
}
|
|
|
|
/*--------------------------------------------------------------------------
|
|
* HYPRE_StructInterfaceMatrixPrint
|
|
*--------------------------------------------------------------------------*/
|
|
|
|
int
|
|
HYPRE_StructInterfaceMatrixPrint( HYPRE_StructInterfaceMatrix matrix )
|
|
{
|
|
return( hypre_PrintStructInterfaceMatrix( (hypre_StructInterfaceMatrix *) matrix ) );
|
|
}
|
|
|
|
/*--------------------------------------------------------------------------
|
|
* HYPRE_SetStructInterfaceMatrixStorageType
|
|
*--------------------------------------------------------------------------*/
|
|
|
|
int
|
|
HYPRE_SetStructInterfaceMatrixStorageType( HYPRE_StructInterfaceMatrix struct_matrix,
|
|
int type )
|
|
{
|
|
return( hypre_SetStructInterfaceMatrixStorageType(
|
|
(hypre_StructInterfaceMatrix *) struct_matrix, type ) );
|
|
}
|
|
|
|
/*--------------------------------------------------------------------------
|
|
* HYPRE_StructInterfaceMatrixSetSymmetric
|
|
*--------------------------------------------------------------------------*/
|
|
|
|
int
|
|
HYPRE_StructInterfaceMatrixSetSymmetric( HYPRE_StructInterfaceMatrix struct_matrix,
|
|
int type )
|
|
{
|
|
return( hypre_SetStructInterfaceMatrixSymmetric(
|
|
(hypre_StructInterfaceMatrix *) struct_matrix, type ) );
|
|
}
|
|
|
|
/*--------------------------------------------------------------------------
|
|
* HYPRE_StructInterfaceMatrixSetNumGhost
|
|
*--------------------------------------------------------------------------*/
|
|
|
|
int
|
|
HYPRE_StructInterfaceMatrixSetNumGhost( HYPRE_StructInterfaceMatrix struct_matrix,
|
|
int *num_ghost )
|
|
/* This does not mean anything for a PETSc matrix under the hood,
|
|
so this is a no-op for now */
|
|
{
|
|
return( 0 );
|
|
}
|
|
|
|
/*--------------------------------------------------------------------------
|
|
* HYPRE_StructInterfaceMatrixGetGrid
|
|
*--------------------------------------------------------------------------*/
|
|
|
|
/*
|
|
Not currently supported.
|
|
*/
|
|
|
|
int
|
|
HYPRE_StructInterfaceMatrixGetGrid( HYPRE_StructInterfaceMatrix matrix,
|
|
HYPRE_StructGrid *grid)
|
|
{
|
|
return( 0 );
|
|
}
|