///////////////////// SHIFT MATRICES ////////////////////// matrix Px;{ int[int] I=(0:Nx*Ny-1); int[int] J=(-1:Nx*Ny-2); real[int] C(Nx*Ny);C=1.; for(int i=0;i