Sorry - I don't reemeber why it was done this way. It could be MPI
stuff [mpi datatpes needing initialization from fortran?] or fortran
command line stuff?
Wrt fortran command line stuff - perhaps that part can be cleaned up
similar to the f90 interface [by using fortran using fortran code -
and not attempting to call getarg/narg and variants directly from C]
Satish
On Tue, 29 Oct 2013, Barry Smith wrote:
>
> Is your question, why is there a
>
> PETSC_EXTERN void PETSC_STDCALL petscinitialize_(CHAR filename
> PETSC_MIXED_LEN(len),PetscErrorCode *ierr PETSC_END_LEN(len))
> {
> #if defined (PETSC_USE_NARGS)
> short flg,i;
> #else
> int i;
> #if !defined(PETSC_HAVE_PXFGETARG_NEW) && !defined (PETSC_HAVE_PXFGETARG_NEW)
> int j;
> #endif
> #endif
>
> function in zstart.c that has a great deal of code duplication (yuck) with
> the C PetscInitialize() function?
>
> The reason is because a long long long time ago we had great difficulty using
> the exact same code for most of the initialization in C and Fortran for a
> bunch of subtle and machine specific reasons I forgot (Satish may remember).
> We introduced the code duplication and have had no problems since (except
> sometimes forgetting to put new things added to PetscInitialize() into
> petscinitialize_() !) Yes ideally we would eliminate the code duplication
> with a refactorization but the gain would mostly be a good feeling and since
> Satish and I sweated a lot of tears long ago with crazy stuff that went wrong
> until we duplicated code I would vote against trying to put them together.
>
> If your question is something else then please restate your question.
>
> Barry
>
>
>
> On Oct 29, 2013, at 4:03 PM, Matthew Knepley <[email protected]> wrote:
>
> > Why do we have a copy of PetscInitialize() in zstart.c? We just call
> > PetscFinalize()
> > from there.
> >
> > Matt
> >
> > --
> > What most experimenters take for granted before they begin their
> > experiments is infinitely more interesting than any results to which their
> > experiments lead.
> > -- Norbert Wiener
>
>
