PetscErrorCode PCGAMGCoarsen_Classical(PC pc,Mat *G,PetscCoarsenData **agg_lists)
{
PetscErrorCode ierr;
MatCoarsen crs;
MPI_Comm fcomm = ((PetscObject)pc)->comm;
PetscFunctionBegin;
/* construct the graph if necessary */
if (!G) {
SETERRQ(fcomm,PETSC_ERR_ARG_WRONGSTATE,"Must set Graph in PC in PCGAMG before coarsening");
}
ierr = MatCoarsenCreate(fcomm,&crs);CHKERRQ(ierr);
ierr = MatCoarsenSetFromOptions(crs);CHKERRQ(ierr);
ierr = MatCoarsenSetAdjacency(crs,*G);CHKERRQ(ierr);
ierr = MatCoarsenSetStrictAggs(crs,PETSC_TRUE);CHKERRQ(ierr);
ierr = MatCoarsenApply(crs);CHKERRQ(ierr);
ierr = MatCoarsenGetData(crs,agg_lists);CHKERRQ(ierr);
ierr = MatCoarsenDestroy(&crs);CHKERRQ(ierr);
PetscFunctionReturn(0);
}
void PETSC_STDCALL matcoarsensetfromoptions_(MatCoarsen *coarser, int *__ierr ){
*__ierr = MatCoarsenSetFromOptions(*coarser);
}
PETSC_EXTERN void PETSC_STDCALL matcoarsensetfromoptions_(MatCoarsen coarser, int *__ierr ){
*__ierr = MatCoarsenSetFromOptions(
(MatCoarsen)PetscToPointer((coarser) ));
}
