#ifndef lint
static char vcid[] = "$Id: isltog.c,v 1.3 1996/11/27 16:56:49 bsmith Exp balay $";
#endif

#include "sys.h"   /*I "sys.h" I*/
#include "is.h"    /*I "is.h"  I*/

#undef __FUNCTION__  
#define __FUNCTION__ "ISLocalToGlobalMappingCreate"
/*@
    ISLocalToGlobalMappingCreate - Creates a mapping between a local (0 to n)
      ordering and a global parallel ordering.

   Input Parameters:
.    n - the number of local elements
.    indices - the global index for each local element

   Output Parameters:
.    mapping - new mapping data structure

.keywords: IS, local to global mapping

.seealso: ISLocalToGlobalMappingDestroy(), 
@*/
int ISLocalToGlobalMappingCreate(int n, int *indices,ISLocalToGlobalMapping *mapping)
{
  PetscValidIntPointer(indices);
  PetscValidPointer(mapping);

  *mapping = PetscNew(struct _ISLocalToGlobalMapping); CHKPTRQ(*mapping);
  (*mapping)->refcnt  = 1;
  (*mapping)->indices = (int *) PetscMalloc((n+1)*sizeof(int));CHKPTRQ((*mapping)->indices);
  PetscMemcpy((*mapping)->indices,indices,n*sizeof(int));
  return 0;
}
  
#undef __FUNCTION__  
#define __FUNCTION__ "ISLocalToGlobalMappingDestroy"
/*@
    ISLocalToGlobalMappingDestroy - Destroys a mapping between a local (0 to n)
      ordering and a global parallel ordering.

   Input Parameters:
.    mapping - mapping data structure

.keywords: IS, local to global mapping

.seealso: ISLocalToGlobalMappingCreate(), 
@*/
int ISLocalToGlobalMappingDestroy(ISLocalToGlobalMapping mapping)
{
  PetscValidPointer(mapping);
  if (--mapping->refcnt) return 0;

  PetscFree(mapping->indices);
  PetscFree(mapping);
  return 0;
}
  
#undef __FUNCTION__  
#define __FUNCTION__ "ISLocalToGlobalMappingApplyIS"
/*@
    ISLocalToGlobalMappingApplyIS - Creates a new IS using the global numbering
      defined in an ISLocalToGlobalMapping from an IS in the local numbering.

   Input Parameters:
.   ISLocalToGlobalMapping - mapping between local and global numbering
.   is - index set in local numbering

   Output Parameters:
.   newis - index set in global numbering

.seealso: ISLocalToGlobalMappingApply(), ISLocalToGlobalMappingCreate(),
          ISLocalToGlobalMappingDestroy()

@*/
int ISLocalToGlobalMappingApplyIS(ISLocalToGlobalMapping mapping, IS is, IS *newis)
{
  int ierr,n,i,*idxin,*idxmap,*idxout;
  PetscValidPointer(mapping);
  PetscValidHeaderSpecific(is,IS_COOKIE);
  PetscValidPointer(newis);

  ierr   = ISGetSize(is,&n); CHKERRQ(ierr);
  ierr   = ISGetIndices(is,&idxin); CHKERRQ(ierr);
  idxmap = mapping->indices;
  
  idxout = (int *) PetscMalloc((n+1)*sizeof(int));CHKPTRQ(idxout);
  for ( i=0; i<n; i++ ) {
    idxout[i] = idxmap[idxin[i]];
  }
  ierr = ISCreateGeneral(MPI_COMM_SELF,n,idxout,newis); CHKERRQ(ierr);
  PetscFree(idxout);
  return 0;
}

#undef __FUNCTION__  
#define __FUNCTION__ "ISLocalToGlobalMappingApply"
/*MC
       ISLocalToGlobalMappingApply - Takes a list of integers in local numbering
              and converts them to global numbering.

   Synopsis:
   void ISLocalToGlobalMappingApply(ISLocalToGlobalMapping mapping,int N,int *in,int *out);

   Input Parameters:
.  mapping - the local to global mapping context
.  N - number of integers
.  in - input indices in local numbering

   Output Parameter:
.  out - indices in global numbering



.seealso: ISLocalToGlobalMappingCreate(),ISLocalToGlobalMappingDestroy(), 
          ISLocalToGlobalMappingApplyIS(),AOCreateDebug(),AOApplicationToPetsc(),
          AOPetscToApplication()

.keywords: local to global, mapping
M*/