static PetscErrorCode Tao_ASLS_FunctionGradient(TaoLineSearch ls, Vec X, PetscReal *fcn, Vec G, void *ptr)
{
Tao tao = (Tao)ptr;
TAO_SSLS *asls = (TAO_SSLS *)tao->data;
PetscErrorCode ierr;
PetscFunctionBegin;
ierr = TaoComputeConstraints(tao, X, tao->constraints);CHKERRQ(ierr);
ierr = VecFischer(X,tao->constraints,tao->XL,tao->XU,asls->ff);CHKERRQ(ierr);
ierr = VecNorm(asls->ff,NORM_2,&asls->merit);CHKERRQ(ierr);
*fcn = 0.5*asls->merit*asls->merit;
ierr = TaoComputeJacobian(tao,tao->solution,tao->jacobian,tao->jacobian_pre);CHKERRQ(ierr);
ierr = MatDFischer(tao->jacobian, tao->solution, tao->constraints,tao->XL, tao->XU, asls->t1, asls->t2,asls->da, asls->db);CHKERRQ(ierr);
ierr = VecPointwiseMult(asls->t1, asls->ff, asls->db);CHKERRQ(ierr);
ierr = MatMultTranspose(tao->jacobian,asls->t1,G);CHKERRQ(ierr);
ierr = VecPointwiseMult(asls->t1, asls->ff, asls->da);CHKERRQ(ierr);
ierr = VecAXPY(G,1.0,asls->t1);CHKERRQ(ierr);
PetscFunctionReturn(0);
}
PetscErrorCode Tao_SSLS_FunctionGradient(TaoLineSearch ls, Vec X, PetscReal *fcn, Vec G, void *ptr)
{
Tao tao = (Tao)ptr;
TAO_SSLS *ssls = (TAO_SSLS *)tao->data;
PetscErrorCode ierr;
PetscFunctionBegin;
ierr = TaoComputeConstraints(tao, X, tao->constraints);CHKERRQ(ierr);
ierr = VecFischer(X,tao->constraints,tao->XL,tao->XU,ssls->ff);CHKERRQ(ierr);
ierr = VecNorm(ssls->ff,NORM_2,&ssls->merit);CHKERRQ(ierr);
*fcn = 0.5*ssls->merit*ssls->merit;
ierr = TaoComputeJacobian(tao,tao->solution,tao->jacobian,tao->jacobian_pre);CHKERRQ(ierr);
ierr = MatDFischer(tao->jacobian, tao->solution, tao->constraints,tao->XL, tao->XU, ssls->t1, ssls->t2,ssls->da, ssls->db);CHKERRQ(ierr);
ierr = MatDiagonalScale(tao->jacobian,ssls->db,NULL);CHKERRQ(ierr);
ierr = MatDiagonalSet(tao->jacobian,ssls->da,ADD_VALUES);CHKERRQ(ierr);
ierr = MatMultTranspose(tao->jacobian,ssls->ff,G);CHKERRQ(ierr);
PetscFunctionReturn(0);
}
