Actual source code: zvecnestf.c

  1: #include <petsc/private/fortranimpl.h>
  2: #include <petscvec.h>

  4: #if defined(PETSC_HAVE_FORTRAN_CAPS)
  5:   #define vecnestgetsubvecs_ VECNESTGETSUBVECS
  6:   #define veccreatenest_     VECCREATENEST
  7: #elif !defined(PETSC_HAVE_FORTRAN_UNDERSCORE)
  8:   #define vecnestgetsubvecs_ vecnestgetsubvecs
  9:   #define veccreatenest_     veccreatenest
 10: #endif

 12: PETSC_EXTERN void vecnestgetsubvecs_(Vec *X, PetscInt *N, Vec *sx, PetscErrorCode *ierr)
 13: {
 14:   Vec     *tsx;
 15:   PetscInt i, n;
 16:   CHKFORTRANNULLINTEGER(N);
 17:   *ierr = VecNestGetSubVecs(*X, &n, &tsx);
 18:   if (*ierr) return;
 19:   if (N) *N = n;
 20:   CHKFORTRANNULLOBJECT(sx);
 21:   if (sx) {
 22:     for (i = 0; i < n; i++) sx[i] = tsx[i];
 23:   }
 24: }

 26: PETSC_EXTERN void veccreatenest_(MPI_Fint *comm, PetscInt *nb, IS is[], Vec x[], Vec *Y, int *ierr)
 27: {
 28:   CHKFORTRANNULLOBJECT(is);
 29:   CHKFORTRANNULLOBJECT(x);
 30:   *ierr = VecCreateNest(MPI_Comm_f2c(*comm), *nb, is, x, Y);
 31: }