Made up contact algorithms for Q3D and S2D solvers and added denum = 0 logic to ShellProduct + fixed ShellProduct bug