THE TWO-GRID AND MULTIGRID DISCRETIZATIONS OF THE C 0 IPG METHOD FOR BIHARMONIC EIGENVALUE PROBLEM