screenpos = (mul(matModelWorldProj, pos).xy*0.5+0.5);

or what should also work is VPOS: http://msdn.microsoft.com/en-us/library/windows/desktop/bb509647(v=vs.85).aspx#VPOS