SubChannel1PhaseProblem.C File Reference

Go to the source code of this file.

Classes
struct	Ctx

Functions
PetscErrorCode	formFunction (SNES, Vec x, Vec f, void *ctx)

Function Documentation

formFunction()

PetscErrorCode formFunction	(	SNES	snes,
		Vec	x,
		Vec	f,
		void *	ctx
	)

Definition at line 33 of file SubChannel1PhaseProblem.C.

{
  const PetscScalar * xx;
  PetscScalar * ff;
  PetscInt size;
 
  PetscFunctionBegin;
  Ctx * cc = static_cast<Ctx *>(ctx);
  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);
}