Numerical study of the square Ising frustration model