p=DDGhole_driftdiffusion(mesh,Dsides,nin,pin,V,un,tn,tp,n0,p0,weight) IN: V = electric potential mesh = integration domain nin = electron density in the past pin = hole density in the past + initial guess n0,p0 = equilibrium densities tn,tp = carrier lifetimes weight = BDF weights up = mobility OUT: p = updated hole density