InterWrapper1PhaseProblem.C File Reference

Go to the source code of this file.

Classes
struct	problemInfo
	problem information to be used in the PETSc problem More...

Functions
PetscErrorCode	formFunctionIW (SNES, Vec x, Vec f, void *info)
	creates the residual function to be used in the PETCs snes context More...

Function Documentation

formFunctionIW()

PetscErrorCode formFunctionIW	(	SNES	,
		Vec	x,
		Vec	f,
		void *	info
	)

creates the residual function to be used in the PETCs snes context

Definition at line 29 of file InterWrapper1PhaseProblem.C.

{
  const PetscScalar * xx;
  PetscScalar * ff;
  PetscInt size;

  PetscFunctionBegin;
  problemInfo * cc = static_cast<problemInfo *>(info);
  LibmeshPetscCallQ(VecGetSize(x, &size));

  libMesh::DenseVector<Real> solution_seed(size, 0.0);
  LibmeshPetscCallQ(VecGetArrayRead(x, &xx));
  for (PetscInt i = 0; i < size; i++)
    solution_seed(i) = xx[i];

  LibmeshPetscCallQ(VecRestoreArrayRead(x, &xx));

  libMesh::DenseVector<Real> Wij_residual_vector =
      cc->schp->residualFunction(cc->iblock, solution_seed);

  LibmeshPetscCallQ(VecGetArray(f, &ff));
  for (int i = 0; i < size; i++)
    ff[i] = Wij_residual_vector(i);

  LibmeshPetscCallQ(VecRestoreArray(f, &ff));
  PetscFunctionReturn(LIBMESH_PETSC_SUCCESS);
}