petsc-3.9.4 2018-09-11

VecFischer

Evaluates the Fischer-Burmeister function for complementarity problems.

#include "petsctao.h" 
PetscErrorCode VecFischer(Vec X, Vec F, Vec L, Vec U, Vec FB)

Logically Collective on vectors

FB -The Fischer-Burmeister function vector

The Fischer-Burmeister function is defined as

       phi(a,b) := sqrt(a*a + b*b) - a - b

and is used reformulate a complementarity problem as a semismooth system of equations.

	l[i] ==	- infinity, u[i] == infinity -- fb[i] = -f[i]
	l[i] ==	- infinity, u[i] finite -- fb[i] = phi(u[i]-x[i], -f[i])
	l[i] finite, u[i] == infinity	- - fb[i] = phi(x[i]-l[i], f[i])
	l[i] finite < u[i] finite	- - fb[i] = phi(x[i]-l[i], phi(u[i]-x[i], -f[u]))
	otherwise l[i] == u[i]	- - fb[i] = l[i] - x[i]

developer