Actual source code: vcreatea.c

petsc-3.7.3 2016-08-01
Report Typos and Errors
  2: #include <../src/sys/classes/viewer/impls/ascii/asciiimpl.h>  /*I     "petscsys.h"   I*/

  4: /* ---------------------------------------------------------------------*/

  6: /*
  7:     The variable Petsc_Viewer_Stdout_keyval is used to indicate an MPI attribute that
  8:   is attached to a communicator, in this case the attribute is a PetscViewer.
  9: */
 10: static PetscMPIInt Petsc_Viewer_Stdout_keyval = MPI_KEYVAL_INVALID;

 14: /*@C
 15:    PetscViewerASCIIGetStdout - Creates a ASCII PetscViewer shared by all processors
 16:                     in a communicator. Error returning version of PETSC_VIEWER_STDOUT_()

 18:    Collective on MPI_Comm

 20:    Input Parameter:
 21: .  comm - the MPI communicator to share the PetscViewer

 23:    Level: beginner

 25:    Notes:
 26:      This should be used in all PETSc source code instead of PETSC_VIEWER_STDOUT_()

 28: .seealso: PETSC_VIEWER_DRAW_(), PetscViewerASCIIOpen(), PETSC_VIEWER_STDERR_, PETSC_VIEWER_STDOUT_WORLD,
 29:           PETSC_VIEWER_STDOUT_SELF

 31: @*/
 32: PetscErrorCode  PetscViewerASCIIGetStdout(MPI_Comm comm,PetscViewer *viewer)
 33: {
 35:   PetscBool      flg;
 36:   MPI_Comm       ncomm;

 39:   PetscSpinlockLock(&PetscViewerASCIISpinLockStdout);
 40:   PetscCommDuplicate(comm,&ncomm,NULL);
 41:   if (Petsc_Viewer_Stdout_keyval == MPI_KEYVAL_INVALID) {
 42:     MPI_Keyval_create(MPI_NULL_COPY_FN,MPI_NULL_DELETE_FN,&Petsc_Viewer_Stdout_keyval,0);
 43:   }
 44:   MPI_Attr_get(ncomm,Petsc_Viewer_Stdout_keyval,(void**)viewer,(PetscMPIInt*)&flg);
 45:   if (!flg) { /* PetscViewer not yet created */
 46:     PetscViewerASCIIOpen(ncomm,"stdout",viewer);
 47:     PetscObjectRegisterDestroy((PetscObject)*viewer);
 48:     MPI_Attr_put(ncomm,Petsc_Viewer_Stdout_keyval,(void*)*viewer);
 49:   }
 50:   PetscCommDestroy(&ncomm);
 51:   PetscSpinlockUnlock(&PetscViewerASCIISpinLockStdout);
 52:   return(0);
 53: }

 57: /*@C
 58:    PETSC_VIEWER_STDOUT_ - Creates a ASCII PetscViewer shared by all processors
 59:                     in a communicator.

 61:    Collective on MPI_Comm

 63:    Input Parameter:
 64: .  comm - the MPI communicator to share the PetscViewer

 66:    Level: beginner

 68:    Notes:
 69:    Unlike almost all other PETSc routines, this does not return
 70:    an error code. Usually used in the form
 71: $      XXXView(XXX object,PETSC_VIEWER_STDOUT_(comm));

 73: .seealso: PETSC_VIEWER_DRAW_(), PetscViewerASCIIOpen(), PETSC_VIEWER_STDERR_, PETSC_VIEWER_STDOUT_WORLD,
 74:           PETSC_VIEWER_STDOUT_SELF

 76: @*/
 77: PetscViewer  PETSC_VIEWER_STDOUT_(MPI_Comm comm)
 78: {
 80:   PetscViewer    viewer;

 83:   PetscViewerASCIIGetStdout(comm,&viewer);
 84:   if (ierr) {PetscError(PETSC_COMM_SELF,__LINE__,"PETSC_VIEWER_STDOUT_",__FILE__,PETSC_ERR_PLIB,PETSC_ERROR_INITIAL," "); return(0);}
 85:   PetscFunctionReturn(viewer);
 86: }

 88: /* ---------------------------------------------------------------------*/

 90: /*
 91:     The variable Petsc_Viewer_Stderr_keyval is used to indicate an MPI attribute that
 92:   is attached to a communicator, in this case the attribute is a PetscViewer.
 93: */
 94: static PetscMPIInt Petsc_Viewer_Stderr_keyval = MPI_KEYVAL_INVALID;

 98: /*@C
 99:    PetscViewerASCIIGetStderr - Creates a ASCII PetscViewer shared by all processors
100:                     in a communicator. Error returning version of PETSC_VIEWER_STDERR_()

102:    Collective on MPI_Comm

104:    Input Parameter:
105: .  comm - the MPI communicator to share the PetscViewer

107:    Level: beginner

109:    Notes:
110:      This should be used in all PETSc source code instead of PETSC_VIEWER_STDERR_()

112: .seealso: PETSC_VIEWER_DRAW_(), PetscViewerASCIIOpen(), PETSC_VIEWER_STDERR_, PETSC_VIEWER_STDERR_WORLD,
113:           PETSC_VIEWER_STDERR_SELF

115: @*/
116: PetscErrorCode  PetscViewerASCIIGetStderr(MPI_Comm comm,PetscViewer *viewer)
117: {
119:   PetscBool      flg;
120:   MPI_Comm       ncomm;

123:   PetscSpinlockLock(&PetscViewerASCIISpinLockStderr);
124:   PetscCommDuplicate(comm,&ncomm,NULL);
125:   if (Petsc_Viewer_Stderr_keyval == MPI_KEYVAL_INVALID) {
126:     MPI_Keyval_create(MPI_NULL_COPY_FN,MPI_NULL_DELETE_FN,&Petsc_Viewer_Stderr_keyval,0);
127:   }
128:   MPI_Attr_get(ncomm,Petsc_Viewer_Stderr_keyval,(void**)viewer,(PetscMPIInt*)&flg);
129:   if (!flg) { /* PetscViewer not yet created */
130:     PetscViewerASCIIOpen(ncomm,"stderr",viewer);
131:     PetscObjectRegisterDestroy((PetscObject)*viewer);
132:     MPI_Attr_put(ncomm,Petsc_Viewer_Stderr_keyval,(void*)*viewer);
133:   }
134:   PetscCommDestroy(&ncomm);
135:   PetscSpinlockUnlock(&PetscViewerASCIISpinLockStderr);
136:   return(0);
137: }

141: /*@C
142:    PETSC_VIEWER_STDERR_ - Creates a ASCII PetscViewer shared by all processors
143:                     in a communicator.

145:    Collective on MPI_Comm

147:    Input Parameter:
148: .  comm - the MPI communicator to share the PetscViewer

150:    Level: beginner

152:    Note:
153:    Unlike almost all other PETSc routines, this does not return
154:    an error code. Usually used in the form
155: $      XXXView(XXX object,PETSC_VIEWER_STDERR_(comm));

157: .seealso: PETSC_VIEWER_DRAW_, PetscViewerASCIIOpen(), PETSC_VIEWER_STDOUT_, PETSC_VIEWER_STDOUT_WORLD,
158:           PETSC_VIEWER_STDOUT_SELF, PETSC_VIEWER_STDERR_WORLD, PETSC_VIEWER_STDERR_SELF
159: @*/
160: PetscViewer  PETSC_VIEWER_STDERR_(MPI_Comm comm)
161: {
163:   PetscViewer    viewer;

166:   PetscViewerASCIIGetStderr(comm,&viewer);
167:   if (ierr) {PetscError(PETSC_COMM_SELF,__LINE__,"PETSC_VIEWER_STDERR_",__FILE__,PETSC_ERR_PLIB,PETSC_ERROR_INITIAL," "); return(0);}
168:   PetscFunctionReturn(viewer);
169: }


172: PetscMPIInt Petsc_Viewer_keyval = MPI_KEYVAL_INVALID;
175: /*
176:    Called with MPI_Comm_free() is called on a communicator that has a viewer as an attribute. The viewer is not actually destroyed because that is managed by
177:    PetscObjectDestroyRegisterAll(). PetscViewerASCIIGetStdout() registers the viewer with PetscObjectDestroyRegister() to be destroyed when PetscFinalize() is called.

179:   This is called by MPI, not by users.

181: */
182: PETSC_EXTERN PetscMPIInt MPIAPI Petsc_DelViewer(MPI_Comm comm,PetscMPIInt keyval,void *attr_val,void *extra_state)
183: {

187:   PetscInfo1(0,"Removing viewer data attribute in an MPI_Comm %ld\n",(long)comm);if (ierr) PetscFunctionReturn((PetscMPIInt)ierr);
188:   PetscFunctionReturn(MPI_SUCCESS);
189: }

193: /*@C
194:    PetscViewerASCIIOpen - Opens an ASCII file as a PetscViewer.

196:    Collective on MPI_Comm

198:    Input Parameters:
199: +  comm - the communicator
200: -  name - the file name

202:    Output Parameter:
203: .  lab - the PetscViewer to use with the specified file

205:    Level: beginner

207:    Notes:
208:    This PetscViewer can be destroyed with PetscViewerDestroy().

210:    The MPI communicator used here must match that used by the object one is viewing. For example if the 
211:    Mat was created with a PETSC_COMM_WORLD, then the Viewer must be created with PETSC_COMM_WORLD

213:    As shown below, PetscViewerASCIIOpen() is useful in conjunction with
214:    MatView() and VecView()
215: .vb
216:      PetscViewerASCIIOpen(PETSC_COMM_WORLD,"mat.output",&viewer);
217:      MatView(matrix,viewer);
218: .ve

220:   Concepts: PetscViewerASCII^creating
221:   Concepts: printf
222:   Concepts: printing
223:   Concepts: accessing remote file
224:   Concepts: remote file

226: .seealso: MatView(), VecView(), PetscViewerDestroy(), PetscViewerBinaryOpen(),
227:           PetscViewerASCIIGetPointer(), PetscViewerPushFormat(), PETSC_VIEWER_STDOUT_, PETSC_VIEWER_STDERR_,
228:           PETSC_VIEWER_STDOUT_WORLD, PETSC_VIEWER_STDOUT_SELF,
229: @*/
230: PetscErrorCode  PetscViewerASCIIOpen(MPI_Comm comm,const char name[],PetscViewer *lab)
231: {
232:   PetscErrorCode  ierr;
233:   PetscViewerLink *vlink,*nv;
234:   PetscBool       flg,eq;
235:   size_t          len;

238:   PetscStrlen(name,&len);
239:   if (!len) {
240:     PetscViewerASCIIGetStdout(comm,lab);
241:     PetscObjectReference((PetscObject)*lab);
242:     return(0);
243:   }
244:   PetscSpinlockLock(&PetscViewerASCIISpinLockOpen);
245:   if (Petsc_Viewer_keyval == MPI_KEYVAL_INVALID) {
246:     MPI_Keyval_create(MPI_NULL_COPY_FN,Petsc_DelViewer,&Petsc_Viewer_keyval,(void*)0);
247:   }
248:   /*
249:        It would be better to move this code to PetscFileSetName() but since it must return a preexiting communicator
250:      we cannot do that, since PetscFileSetName() takes a communicator that already exists.

252:       Plus if the original communicator that created the file has since been close this will not detect the old
253:       communictor and hence will overwrite the old data. It may be better to simply remove all this code
254:   */
255:   /* make sure communicator is a PETSc communicator */
256:   PetscCommDuplicate(comm,&comm,NULL);
257:   /* has file already been opened into a viewer */
258:   MPI_Attr_get(comm,Petsc_Viewer_keyval,(void**)&vlink,(PetscMPIInt*)&flg);
259:   if (flg) {
260:     while (vlink) {
261:       PetscStrcmp(name,((PetscViewer_ASCII*)(vlink->viewer->data))->filename,&eq);
262:       if (eq) {
263:         PetscObjectReference((PetscObject)vlink->viewer);
264:         *lab = vlink->viewer;
265:         PetscCommDestroy(&comm);
266:         PetscSpinlockUnlock(&PetscViewerASCIISpinLockOpen);
267:         return(0);
268:       }
269:       vlink = vlink->next;
270:     }
271:   }
272:   PetscViewerCreate(comm,lab);
273:   PetscViewerSetType(*lab,PETSCVIEWERASCII);
274:   if (name) {
275:     PetscViewerFileSetName(*lab,name);
276:   }
277:   /* save viewer into communicator if needed later */
278:   PetscNew(&nv);
279:   nv->viewer = *lab;
280:   if (!flg) {
281:     MPI_Attr_put(comm,Petsc_Viewer_keyval,nv);
282:   } else {
283:     MPI_Attr_get(comm,Petsc_Viewer_keyval,(void**)&vlink,(PetscMPIInt*)&flg);
284:     if (vlink) {
285:       while (vlink->next) vlink = vlink->next;
286:       vlink->next = nv;
287:     } else {
288:       MPI_Attr_put(comm,Petsc_Viewer_keyval,nv);
289:     }
290:   }
291:   PetscCommDestroy(&comm);
292:   PetscSpinlockUnlock(&PetscViewerASCIISpinLockOpen);
293:   return(0);
294: }

298: /*@C
299:    PetscViewerASCIIOpenWithFILE - Given an open file creates an ASCII viewer that prints to it.

301:    Collective on MPI_Comm

303:    Input Parameters:
304: +  comm - the communicator
305: -  fd - the FILE pointer

307:    Output Parameter:
308: .  lab - the PetscViewer to use with the specified file

310:    Level: beginner

312:    Notes:
313:    This PetscViewer can be destroyed with PetscViewerDestroy(), but the fd will NOT be closed.

315:    If a multiprocessor communicator is used (such as PETSC_COMM_WORLD),
316:    then only the first processor in the group uses the file.  All other
317:    processors send their data to the first processor to print.

319:   Concepts: PetscViewerASCII^creating
320:   Concepts: printf
321:   Concepts: printing
322:   Concepts: accessing remote file
323:   Concepts: remote file

325: .seealso: MatView(), VecView(), PetscViewerDestroy(), PetscViewerBinaryOpen(),
326:           PetscViewerASCIIGetPointer(), PetscViewerPushFormat(), PETSC_VIEWER_STDOUT_, PETSC_VIEWER_STDERR_,
327:           PETSC_VIEWER_STDOUT_WORLD, PETSC_VIEWER_STDOUT_SELF, PetscViewerASCIIOpen()
328: @*/
329: PetscErrorCode  PetscViewerASCIIOpenWithFILE(MPI_Comm comm,FILE *fd,PetscViewer *lab)
330: {

334:   PetscViewerCreate(comm,lab);
335:   PetscViewerSetType(*lab,PETSCVIEWERASCII);
336:   PetscViewerASCIISetFILE(*lab,fd);
337:   return(0);
338: }

342: PetscErrorCode  PetscViewerASCIISetFILE(PetscViewer viewer,FILE *fd)
343: {
344:   PetscViewer_ASCII *vascii = (PetscViewer_ASCII*)viewer->data;

347:   vascii->fd        = fd;
348:   vascii->closefile = PETSC_FALSE;
349:   return(0);
350: }