15 PetscErrorCode PipeCreate(MPI_Comm comm, Pipe *pipe)  in PipeCreate()  argument
18   PetscCall(PetscNew(pipe));  in PipeCreate()
28 PetscErrorCode PipeDestroy(Pipe *pipe)  in PipeDestroy()  argument
31   if (!*pipe) PetscFunctionReturn(PETSC_SUCCESS);  in PipeDestroy()
33   PetscCall(PipeDestroyJacobian(*pipe));  in PipeDestroy()
34   PetscCall(VecDestroy(&(*pipe)->x));  in PipeDestroy()
35   PetscCall(DMDestroy(&(*pipe)->da));  in PipeDestroy()
50 PetscErrorCode PipeSetParameters(Pipe pipe, PetscReal length, PetscReal D, PetscReal a, PetscReal f…  in PipeSetParameters()  argument
53   pipe->length = length;  in PipeSetParameters()
54   pipe->D      = D;  in PipeSetParameters()
55   pipe->a      = a;  in PipeSetParameters()
56   pipe->fric   = fric;  in PipeSetParameters()
63 PetscErrorCode PipeSetUp(Pipe pipe)  in PipeSetUp()  argument
68 …PetscCall(DMDACreate1d(PETSC_COMM_SELF, DM_BOUNDARY_GHOSTED, pipe->nnodes, 2, 1, NULL, &pipe->da));  in PipeSetUp()
69   PetscCall(DMSetFromOptions(pipe->da));  in PipeSetUp()
70   PetscCall(DMSetUp(pipe->da));  in PipeSetUp()
71   PetscCall(DMDASetFieldName(pipe->da, 0, "Q"));  in PipeSetUp()
72   PetscCall(DMDASetFieldName(pipe->da, 1, "H"));  in PipeSetUp()
73   PetscCall(DMDASetUniformCoordinates(pipe->da, 0, pipe->length, 0, 0, 0, 0));  in PipeSetUp()
74   PetscCall(DMCreateGlobalVector(pipe->da, &pipe->x));  in PipeSetUp()
76   PetscCall(DMDAGetLocalInfo(pipe->da, &info));  in PipeSetUp()
78   pipe->rad = pipe->D / 2;  in PipeSetUp()
79   pipe->A   = PETSC_PI * pipe->rad * pipe->rad;  in PipeSetUp()
80   pipe->R   = pipe->fric / (2 * pipe->D * pipe->A);  in PipeSetUp()
98 PetscErrorCode PipeCreateJacobian(Pipe pipe, Mat *Jin, Mat *J[])  in PipeCreateJacobian()  argument
107     pipe->jacobian = Jin;  in PipeCreateJacobian()
114   PetscCall(DMSetMatrixStructureOnly(pipe->da, PETSC_TRUE));  in PipeCreateJacobian()
115   PetscCall(DMCreateMatrix(pipe->da, &Jpipe[0]));  in PipeCreateJacobian()
116   PetscCall(DMSetMatrixStructureOnly(pipe->da, PETSC_FALSE));  in PipeCreateJacobian()
156   pipe->jacobian = Jpipe;  in PipeCreateJacobian()
160 PetscErrorCode PipeDestroyJacobian(Pipe pipe)  in PipeDestroyJacobian()  argument
162   Mat     *Jpipe = pipe->jacobian;  in PipeDestroyJacobian()