15c6c1daeSBarry Smith 25c6c1daeSBarry Smith #include <petsc-private/viewerimpl.h> /*I "petscsys.h" I*/ 35c6c1daeSBarry Smith 45c6c1daeSBarry Smith PetscFList PetscViewerList = 0; 55c6c1daeSBarry Smith 62bf49c77SBarry Smith PetscErrorCode PetscOptionsFindPair_Private(const char[],const char[],char *[],PetscBool*); 72bf49c77SBarry Smith #undef __FUNCT__ 82bf49c77SBarry Smith #define __FUNCT__ "PetscOptionsGetViewer" 92bf49c77SBarry Smith /*@C 102bf49c77SBarry Smith PetscOptionsGetViewer - Gets a viewer appropriate for the type indicated by the user 112bf49c77SBarry Smith 122bf49c77SBarry Smith Collective on MPI_Comm 132bf49c77SBarry Smith 142bf49c77SBarry Smith Input Parameters: 152bf49c77SBarry Smith + comm - the communicator to own the viewer 162bf49c77SBarry Smith . pre - the string to prepend to the name or PETSC_NULL 172bf49c77SBarry Smith - name - the option one is seeking 182bf49c77SBarry Smith 192bf49c77SBarry Smith Output Parameter: 202bf49c77SBarry Smith + viewer - the viewer 212bf49c77SBarry Smith - set - PETSC_TRUE if found, else PETSC_FALSE 222bf49c77SBarry Smith 232bf49c77SBarry Smith Level: intermediate 242bf49c77SBarry Smith 252bf49c77SBarry Smith Notes: If no value is provided ascii:stdout is used 267f677774SBarry Smith $ ascii[:[filename][:format]] defaults to stdout - format can be one of info, info_detailed, or matlab, for example ascii::info prints just the info 277f677774SBarry Smith $ about the object to standard out 282bf49c77SBarry Smith $ binary[:filename] defaults to binaryoutput 292bf49c77SBarry Smith $ draw 302bf49c77SBarry Smith $ socket[:port] defaults to the standard output port 312bf49c77SBarry Smith 322bf49c77SBarry Smith Use PetscOptionsRestoreViewer() after using the viewer, otherwise a memory leak may occur 332bf49c77SBarry Smith 342bf49c77SBarry Smith .seealso: PetscOptionsGetReal(), PetscOptionsHasName(), PetscOptionsGetString(), 352bf49c77SBarry Smith PetscOptionsGetIntArray(), PetscOptionsGetRealArray(), PetscOptionsBool() 362bf49c77SBarry Smith PetscOptionsInt(), PetscOptionsString(), PetscOptionsReal(), PetscOptionsBool(), 372bf49c77SBarry Smith PetscOptionsName(), PetscOptionsBegin(), PetscOptionsEnd(), PetscOptionsHead(), 382bf49c77SBarry Smith PetscOptionsStringArray(),PetscOptionsRealArray(), PetscOptionsScalar(), 392bf49c77SBarry Smith PetscOptionsBoolGroupBegin(), PetscOptionsBoolGroup(), PetscOptionsBoolGroupEnd(), 402bf49c77SBarry Smith PetscOptionsList(), PetscOptionsEList(), PetscOptionsRestoreViewer() 412bf49c77SBarry Smith @*/ 422bf49c77SBarry Smith PetscErrorCode PetscOptionsGetViewer(MPI_Comm comm,const char pre[],const char name[],PetscViewer *viewer,PetscBool *set) 432bf49c77SBarry Smith { 442bf49c77SBarry Smith char *value; 452bf49c77SBarry Smith PetscErrorCode ierr; 462bf49c77SBarry Smith PetscBool flag; 472bf49c77SBarry Smith 482bf49c77SBarry Smith PetscFunctionBegin; 492bf49c77SBarry Smith PetscValidCharPointer(name,3); 502bf49c77SBarry Smith 512bf49c77SBarry Smith if (set) *set = PETSC_FALSE; 522bf49c77SBarry Smith ierr = PetscOptionsFindPair_Private(pre,name,&value,&flag);CHKERRQ(ierr); 532bf49c77SBarry Smith if (flag) { 542bf49c77SBarry Smith if (set) *set = PETSC_TRUE; 552bf49c77SBarry Smith if (!value) { 562bf49c77SBarry Smith ierr = PetscViewerASCIIGetStdout(comm,viewer);CHKERRQ(ierr); 572bf49c77SBarry Smith } else { 587f677774SBarry Smith char *cvalue,*loc,*loc2 = PETSC_NULL; 592bf49c77SBarry Smith PetscInt cnt; 607f677774SBarry Smith const char *viewers[] = {PETSCVIEWERASCII,PETSCVIEWERBINARY,PETSCVIEWERDRAW,PETSCVIEWERSOCKET,0}; 612bf49c77SBarry Smith 622bf49c77SBarry Smith ierr = PetscStrallocpy(value,&cvalue);CHKERRQ(ierr); 632bf49c77SBarry Smith ierr = PetscStrchr(cvalue,':',&loc);CHKERRQ(ierr); 642bf49c77SBarry Smith if (loc) {*loc = 0; loc++;} 657f677774SBarry Smith ierr = PetscStrendswithwhich(*cvalue ? cvalue : "ascii",viewers,&cnt);CHKERRQ(ierr); 662bf49c77SBarry Smith if (cnt == 4) SETERRQ1(comm,PETSC_ERR_ARG_OUTOFRANGE,"Unknown viewer type: %s",cvalue); 672bf49c77SBarry Smith if (!loc) { 687f677774SBarry Smith if (cnt == 0) { 692bf49c77SBarry Smith ierr = PetscViewerASCIIGetStdout(comm,viewer);CHKERRQ(ierr); 702bf49c77SBarry Smith } 717f677774SBarry Smith if (cnt == 1) { 727f677774SBarry Smith *viewer = PETSC_VIEWER_BINARY_(comm);CHKERRQ(ierr); 732bf49c77SBarry Smith } 747f677774SBarry Smith if (cnt == 2) { 757f677774SBarry Smith *viewer = PETSC_VIEWER_DRAW_(comm);CHKERRQ(ierr); 767f677774SBarry Smith } 777f677774SBarry Smith if (cnt == 3) { 787f677774SBarry Smith *viewer = PETSC_VIEWER_SOCKET_(comm);CHKERRQ(ierr); 797f677774SBarry Smith } 807f677774SBarry Smith } else { 817f677774SBarry Smith ierr = PetscStrchr(loc,':',&loc2);CHKERRQ(ierr); 827f677774SBarry Smith if (loc2) {*loc2 = 0; loc2++;} 837f677774SBarry Smith if (loc2 && !*loc) { /* ASCII format without file name */ 847f677774SBarry Smith ierr = PetscViewerASCIIGetStdout(comm,viewer);CHKERRQ(ierr); 857f677774SBarry Smith } else { 862bf49c77SBarry Smith ierr = PetscViewerCreate(comm,viewer);CHKERRQ(ierr); 877f677774SBarry Smith ierr = PetscViewerSetType(*viewer,*cvalue ? cvalue : "ascii");CHKERRQ(ierr); 882bf49c77SBarry Smith ierr = PetscViewerFileSetMode(*viewer,FILE_MODE_WRITE);CHKERRQ(ierr); 892bf49c77SBarry Smith ierr = PetscViewerFileSetName(*viewer,loc);CHKERRQ(ierr); 902bf49c77SBarry Smith ierr = PetscObjectComposeFunction((PetscObject)*viewer,"PetscOptionsDestroyViewer","PetscViewerDestroy",(void (*)(void))PetscViewerDestroy);CHKERRQ(ierr); 917f677774SBarry Smith } 927f677774SBarry Smith } 932bf49c77SBarry Smith ierr = PetscViewerSetUp(*viewer);CHKERRQ(ierr); 94*2f234a98SBarry Smith if (loc2 && *loc2) { 95*2f234a98SBarry Smith ierr = PetscStrtoupper(loc2);CHKERRQ(ierr); 96*2f234a98SBarry Smith ierr = PetscStrendswithwhich(loc2,PetscViewerFormats,&cnt);CHKERRQ(ierr); 97*2f234a98SBarry Smith if (!PetscViewerFormats[cnt]) SETERRQ1(PETSC_COMM_SELF,PETSC_ERR_SUP,"Unknown viewer format %s",loc2);CHKERRQ(ierr); 98*2f234a98SBarry Smith ierr = PetscViewerPushFormat(*viewer,(PetscViewerFormat)cnt);CHKERRQ(ierr); 99*2f234a98SBarry Smith ierr = PetscObjectComposeFunction((PetscObject)*viewer,"PetscOptionsPopViewer","PetscViewerPopFormat",(void (*)(void))PetscViewerPopFormat);CHKERRQ(ierr); 100*2f234a98SBarry Smith } 101*2f234a98SBarry Smith ierr = PetscFree(cvalue);CHKERRQ(ierr); 1022bf49c77SBarry Smith } 1032bf49c77SBarry Smith } 1042bf49c77SBarry Smith PetscFunctionReturn(0); 1052bf49c77SBarry Smith } 1062bf49c77SBarry Smith 1072bf49c77SBarry Smith #undef __FUNCT__ 1082bf49c77SBarry Smith #define __FUNCT__ "PetscOptionsRestoreViewer" 1092bf49c77SBarry Smith /*@C 1102bf49c77SBarry Smith PetscOptionsRestoresViewer - Restores a viewer obtained with PetscOptionsGetViewer() 1112bf49c77SBarry Smith 1122bf49c77SBarry Smith Collective on PetscViewer 1132bf49c77SBarry Smith 1142bf49c77SBarry Smith Input Parameters: 1152bf49c77SBarry Smith . viewer - the viewer 1162bf49c77SBarry Smith 1172bf49c77SBarry Smith Level: intermediate 1182bf49c77SBarry Smith 1192bf49c77SBarry Smith .seealso: PetscOptionsGetReal(), PetscOptionsHasName(), PetscOptionsGetString(), 1202bf49c77SBarry Smith PetscOptionsGetIntArray(), PetscOptionsGetRealArray(), PetscOptionsBool() 1212bf49c77SBarry Smith PetscOptionsInt(), PetscOptionsString(), PetscOptionsReal(), PetscOptionsBool(), 1222bf49c77SBarry Smith PetscOptionsName(), PetscOptionsBegin(), PetscOptionsEnd(), PetscOptionsHead(), 1232bf49c77SBarry Smith PetscOptionsStringArray(),PetscOptionsRealArray(), PetscOptionsScalar(), 1242bf49c77SBarry Smith PetscOptionsBoolGroupBegin(), PetscOptionsBoolGroup(), PetscOptionsBoolGroupEnd(), 1252bf49c77SBarry Smith PetscOptionsList(), PetscOptionsEList(), PetscOptionsRestoreViewer() 1262bf49c77SBarry Smith @*/ 1272bf49c77SBarry Smith PetscErrorCode PetscOptionsRestoreViewer(PetscViewer viewer) 1282bf49c77SBarry Smith { 129*2f234a98SBarry Smith PetscErrorCode ierr,(*f)(PetscViewer*),(*g)(PetscViewer); 1302bf49c77SBarry Smith 1312bf49c77SBarry Smith PetscFunctionBegin; 132*2f234a98SBarry Smith ierr = PetscObjectQueryFunction((PetscObject)viewer,"PetscOptionsPopViewer",(void (**)(void))&g);CHKERRQ(ierr); 133*2f234a98SBarry Smith if (g) { 134*2f234a98SBarry Smith ierr = (*g)(viewer);CHKERRQ(ierr); 135*2f234a98SBarry Smith } 1362bf49c77SBarry Smith ierr = PetscObjectQueryFunction((PetscObject)viewer,"PetscOptionsDestroyViewer",(void (**)(void))&f);CHKERRQ(ierr); 1372bf49c77SBarry Smith if (f) { 1382bf49c77SBarry Smith ierr = (*f)(&viewer);CHKERRQ(ierr); 1392bf49c77SBarry Smith } 1402bf49c77SBarry Smith PetscFunctionReturn(0); 1412bf49c77SBarry Smith } 1422bf49c77SBarry Smith 1435c6c1daeSBarry Smith #undef __FUNCT__ 1445c6c1daeSBarry Smith #define __FUNCT__ "PetscViewerCreate" 1455c6c1daeSBarry Smith /*@ 1465c6c1daeSBarry Smith PetscViewerCreate - Creates a viewing context 1475c6c1daeSBarry Smith 1485c6c1daeSBarry Smith Collective on MPI_Comm 1495c6c1daeSBarry Smith 1505c6c1daeSBarry Smith Input Parameter: 1515c6c1daeSBarry Smith . comm - MPI communicator 1525c6c1daeSBarry Smith 1535c6c1daeSBarry Smith Output Parameter: 1545c6c1daeSBarry Smith . inviewer - location to put the PetscViewer context 1555c6c1daeSBarry Smith 1565c6c1daeSBarry Smith Level: advanced 1575c6c1daeSBarry Smith 1585c6c1daeSBarry Smith Concepts: graphics^creating PetscViewer 1595c6c1daeSBarry Smith Concepts: file input/output^creating PetscViewer 1605c6c1daeSBarry Smith Concepts: sockets^creating PetscViewer 1615c6c1daeSBarry Smith 1625c6c1daeSBarry Smith .seealso: PetscViewerDestroy(), PetscViewerSetType(), PetscViewerType 1635c6c1daeSBarry Smith 1645c6c1daeSBarry Smith @*/ 1655c6c1daeSBarry Smith PetscErrorCode PetscViewerCreate(MPI_Comm comm,PetscViewer *inviewer) 1665c6c1daeSBarry Smith { 1675c6c1daeSBarry Smith PetscViewer viewer; 1685c6c1daeSBarry Smith PetscErrorCode ierr; 1695c6c1daeSBarry Smith 1705c6c1daeSBarry Smith PetscFunctionBegin; 1715c6c1daeSBarry Smith *inviewer = 0; 1725c6c1daeSBarry Smith #if !defined(PETSC_USE_DYNAMIC_LIBRARIES) 1735c6c1daeSBarry Smith ierr = PetscViewerInitializePackage(PETSC_NULL);CHKERRQ(ierr); 1745c6c1daeSBarry Smith #endif 1755c6c1daeSBarry Smith ierr = PetscHeaderCreate(viewer,_p_PetscViewer,struct _PetscViewerOps,PETSC_VIEWER_CLASSID,-1,"PetscViewer","PetscViewer","Viewer",comm,PetscViewerDestroy,0);CHKERRQ(ierr); 1765c6c1daeSBarry Smith *inviewer = viewer; 1775c6c1daeSBarry Smith viewer->data = 0; 1785c6c1daeSBarry Smith PetscFunctionReturn(0); 1795c6c1daeSBarry Smith } 1805c6c1daeSBarry Smith 1815c6c1daeSBarry Smith #undef __FUNCT__ 1825c6c1daeSBarry Smith #define __FUNCT__ "PetscViewerSetType" 1835c6c1daeSBarry Smith /*@C 1845c6c1daeSBarry Smith PetscViewerSetType - Builds PetscViewer for a particular implementation. 1855c6c1daeSBarry Smith 1865c6c1daeSBarry Smith Collective on PetscViewer 1875c6c1daeSBarry Smith 1885c6c1daeSBarry Smith Input Parameter: 1895c6c1daeSBarry Smith + viewer - the PetscViewer context 1905c6c1daeSBarry Smith - type - for example, "ASCII" 1915c6c1daeSBarry Smith 1925c6c1daeSBarry Smith Options Database Command: 1935c6c1daeSBarry Smith . -draw_type <type> - Sets the type; use -help for a list 1945c6c1daeSBarry Smith of available methods (for instance, ascii) 1955c6c1daeSBarry Smith 1965c6c1daeSBarry Smith Level: advanced 1975c6c1daeSBarry Smith 1985c6c1daeSBarry Smith Notes: 1995c6c1daeSBarry Smith See "include/petscviewer.h" for available methods (for instance, 2005c6c1daeSBarry Smith PETSC_VIEWER_SOCKET) 2015c6c1daeSBarry Smith 2025c6c1daeSBarry Smith .seealso: PetscViewerCreate(), PetscViewerGetType(), PetscViewerType 2035c6c1daeSBarry Smith @*/ 2045c6c1daeSBarry Smith PetscErrorCode PetscViewerSetType(PetscViewer viewer,PetscViewerType type) 2055c6c1daeSBarry Smith { 2065c6c1daeSBarry Smith PetscErrorCode ierr,(*r)(PetscViewer); 2075c6c1daeSBarry Smith PetscBool match; 2085c6c1daeSBarry Smith 2095c6c1daeSBarry Smith PetscFunctionBegin; 2105c6c1daeSBarry Smith PetscValidHeaderSpecific(viewer,PETSC_VIEWER_CLASSID,1); 2115c6c1daeSBarry Smith PetscValidCharPointer(type,2); 2125c6c1daeSBarry Smith CHKMEMQ; 2135c6c1daeSBarry Smith ierr = PetscObjectTypeCompare((PetscObject)viewer,type,&match);CHKERRQ(ierr); 2145c6c1daeSBarry Smith if (match) PetscFunctionReturn(0); 2155c6c1daeSBarry Smith 2165c6c1daeSBarry Smith /* cleanup any old type that may be there */ 2175c6c1daeSBarry Smith if (viewer->data) { 2185c6c1daeSBarry Smith ierr = (*viewer->ops->destroy)(viewer);CHKERRQ(ierr); 2195c6c1daeSBarry Smith viewer->ops->destroy = PETSC_NULL; 2205c6c1daeSBarry Smith viewer->data = 0; 2215c6c1daeSBarry Smith } 2225c6c1daeSBarry Smith ierr = PetscMemzero(viewer->ops,sizeof(struct _PetscViewerOps));CHKERRQ(ierr); 2235c6c1daeSBarry Smith 2245c6c1daeSBarry Smith ierr = PetscFListFind(PetscViewerList,((PetscObject)viewer)->comm,type,PETSC_TRUE,(void (**)(void)) &r);CHKERRQ(ierr); 2255c6c1daeSBarry Smith if (!r) SETERRQ1(PETSC_COMM_SELF,PETSC_ERR_ARG_UNKNOWN_TYPE,"Unknown PetscViewer type given: %s",type); 2265c6c1daeSBarry Smith 2275c6c1daeSBarry Smith ierr = PetscObjectChangeTypeName((PetscObject)viewer,type);CHKERRQ(ierr); 2285c6c1daeSBarry Smith ierr = (*r)(viewer);CHKERRQ(ierr); 2295c6c1daeSBarry Smith PetscFunctionReturn(0); 2305c6c1daeSBarry Smith } 2315c6c1daeSBarry Smith 2325c6c1daeSBarry Smith #undef __FUNCT__ 2335c6c1daeSBarry Smith #define __FUNCT__ "PetscViewerRegisterDestroy" 2345c6c1daeSBarry Smith /*@C 2355c6c1daeSBarry Smith PetscViewerRegisterDestroy - Frees the list of PetscViewer methods that were 2365c6c1daeSBarry Smith registered by PetscViewerRegisterDynamic(). 2375c6c1daeSBarry Smith 2385c6c1daeSBarry Smith Not Collective 2395c6c1daeSBarry Smith 2405c6c1daeSBarry Smith Level: developer 2415c6c1daeSBarry Smith 2425c6c1daeSBarry Smith .seealso: PetscViewerRegisterDynamic(), PetscViewerRegisterAll() 2435c6c1daeSBarry Smith @*/ 2445c6c1daeSBarry Smith PetscErrorCode PetscViewerRegisterDestroy(void) 2455c6c1daeSBarry Smith { 2465c6c1daeSBarry Smith PetscErrorCode ierr; 2475c6c1daeSBarry Smith 2485c6c1daeSBarry Smith PetscFunctionBegin; 2495c6c1daeSBarry Smith ierr = PetscFListDestroy(&PetscViewerList);CHKERRQ(ierr); 2505c6c1daeSBarry Smith PetscFunctionReturn(0); 2515c6c1daeSBarry Smith } 2525c6c1daeSBarry Smith 2535c6c1daeSBarry Smith #undef __FUNCT__ 2545c6c1daeSBarry Smith #define __FUNCT__ "PetscViewerRegister" 2555c6c1daeSBarry Smith PetscErrorCode PetscViewerRegister(const char *sname,const char *path,const char *name,PetscErrorCode (*function)(PetscViewer)) 2565c6c1daeSBarry Smith { 2575c6c1daeSBarry Smith PetscErrorCode ierr; 2585c6c1daeSBarry Smith char fullname[PETSC_MAX_PATH_LEN]; 2595c6c1daeSBarry Smith 2605c6c1daeSBarry Smith PetscFunctionBegin; 2615c6c1daeSBarry Smith ierr = PetscFListConcat(path,name,fullname);CHKERRQ(ierr); 2625c6c1daeSBarry Smith ierr = PetscFListAdd(&PetscViewerList,sname,fullname,(void (*)(void))function);CHKERRQ(ierr); 2635c6c1daeSBarry Smith PetscFunctionReturn(0); 2645c6c1daeSBarry Smith } 2655c6c1daeSBarry Smith 2665c6c1daeSBarry Smith #undef __FUNCT__ 2675c6c1daeSBarry Smith #define __FUNCT__ "PetscViewerSetFromOptions" 2685c6c1daeSBarry Smith /*@C 2695c6c1daeSBarry Smith PetscViewerSetFromOptions - Sets the graphics type from the options database. 2705c6c1daeSBarry Smith Defaults to a PETSc X windows graphics. 2715c6c1daeSBarry Smith 2725c6c1daeSBarry Smith Collective on PetscViewer 2735c6c1daeSBarry Smith 2745c6c1daeSBarry Smith Input Parameter: 2755c6c1daeSBarry Smith . PetscViewer - the graphics context 2765c6c1daeSBarry Smith 2775c6c1daeSBarry Smith Level: intermediate 2785c6c1daeSBarry Smith 2795c6c1daeSBarry Smith Notes: 2805c6c1daeSBarry Smith Must be called after PetscViewerCreate() before the PetscViewer is used. 2815c6c1daeSBarry Smith 2825c6c1daeSBarry Smith Concepts: PetscViewer^setting options 2835c6c1daeSBarry Smith 2845c6c1daeSBarry Smith .seealso: PetscViewerCreate(), PetscViewerSetType(), PetscViewerType 2855c6c1daeSBarry Smith 2865c6c1daeSBarry Smith @*/ 2875c6c1daeSBarry Smith PetscErrorCode PetscViewerSetFromOptions(PetscViewer viewer) 2885c6c1daeSBarry Smith { 2895c6c1daeSBarry Smith PetscErrorCode ierr; 2905c6c1daeSBarry Smith char vtype[256]; 2915c6c1daeSBarry Smith PetscBool flg; 2925c6c1daeSBarry Smith 2935c6c1daeSBarry Smith PetscFunctionBegin; 2945c6c1daeSBarry Smith PetscValidHeaderSpecific(viewer,PETSC_VIEWER_CLASSID,1); 2955c6c1daeSBarry Smith 2965c6c1daeSBarry Smith if (!PetscViewerList) { 2975c6c1daeSBarry Smith ierr = PetscViewerRegisterAll(PETSC_NULL);CHKERRQ(ierr); 2985c6c1daeSBarry Smith } 2995c6c1daeSBarry Smith ierr = PetscObjectOptionsBegin((PetscObject)viewer);CHKERRQ(ierr); 3005c6c1daeSBarry Smith ierr = PetscOptionsList("-viewer_type","Type of PetscViewer","None",PetscViewerList,(char *)(((PetscObject)viewer)->type_name?((PetscObject)viewer)->type_name:PETSCVIEWERASCII),vtype,256,&flg);CHKERRQ(ierr); 3015c6c1daeSBarry Smith if (flg) { 3025c6c1daeSBarry Smith ierr = PetscViewerSetType(viewer,vtype);CHKERRQ(ierr); 3035c6c1daeSBarry Smith } 3045c6c1daeSBarry Smith /* type has not been set? */ 3055c6c1daeSBarry Smith if (!((PetscObject)viewer)->type_name) { 3065c6c1daeSBarry Smith ierr = PetscViewerSetType(viewer,PETSCVIEWERASCII);CHKERRQ(ierr); 3075c6c1daeSBarry Smith } 3085c6c1daeSBarry Smith if (viewer->ops->setfromoptions) { 3095c6c1daeSBarry Smith ierr = (*viewer->ops->setfromoptions)(viewer);CHKERRQ(ierr); 3105c6c1daeSBarry Smith } 3115c6c1daeSBarry Smith 3125c6c1daeSBarry Smith /* process any options handlers added with PetscObjectAddOptionsHandler() */ 3135c6c1daeSBarry Smith ierr = PetscObjectProcessOptionsHandlers((PetscObject)viewer);CHKERRQ(ierr); 3145c6c1daeSBarry Smith ierr = PetscOptionsEnd();CHKERRQ(ierr); 3155c6c1daeSBarry Smith 3165c6c1daeSBarry Smith PetscFunctionReturn(0); 3175c6c1daeSBarry Smith } 318