Function File: [jx,jy] = Ufvsgcurrent3 (mesh, u, alpha, gamma, eta, beta);

Builds the Scharfetter-Gummel approximation of the vector field

where:

• alpha is an element-wise constant scalar function
• eta, u, gamma are piecewise linear conforming scalar functions
• beta is an element-wise constant vector function

J(u) is an element-wise constant vector function

Instead of passing the vector field beta directly one can pass a piecewise linear conforming scalar function phi as the last input. In such case beta = grad phi is assumed. If phi is a single scalar value beta is assumed to be 0 in the whole domain.

See also: Uscharfettergummel3