function p=DDGOXhole_driftdiffusion(mesh,Dsides,nin,pin,V,up,tn,tp,n0,p0)
DESCRIPTION
%
p=DDGhole_driftdiffusion(mesh,Dsides,nin,pin,V,up,tn,tp,n0,p0)
IN:
v = electric potential
mesh = integration domain
nin = initial guess and BCs for electron density
pin = hole density (to compute SRH recombination)
OUT:
p = updated hole density
%