1-dimension wave equation을 이용한 모델링을 구현중입니다.

포트란코드를 C++코드로 변환하고 있는데... 어려움이 있어 질문드립니다...
★로 표시되어 있는 부분에서 * 를 저는 곱셈연산자로 사용하였는데, 컴파일하게되면

        error C2112 : '-' : 포인터 빼기에는 정수 계열 또는 포인터 피연산자가 있어야 합니다.

라고 오류가 뜨네요.
포인터배열로 선언한것도 아닌데 a1에서는 오류가 없고 a3, a4 계산에서만 오류가 뜨는 이유를 모르겠네요.

===========================================================================================

#include <stdio.h>
#include <fstream>
#include <iostream>
#include <cmath>
using namespace std;

int main(){
    #define nx 401
    #define nz 201
    #define nt 3000   

    int ix, iz, it, srcx, srcz;
    double dt, dx, dz, a1, a2, a3, a4, a5, temp;

    double u1[nx+2][nz+2]={0}, u2[nx+2][nz+2]={0}, u3[nx+2][nz+2]={0};
    double v[nx+1][nz+1]={0}, f[nx+1][nz+1]={0}, w[nt+1];

    (중략)

        for (iz=1;iz<=nz;iz++){
            for (ix=1;ix<=nx;ix++){
                a1=2*u2[ix][iz] - u1[ix][iz];
                a2=(v[ix][iz]*v[ix][iz])*(dt * dt);
                a3=(u2[ix+1,iz]-2*u2[ix][iz]+u2[ix-1][iz])/(dx * dx); ★
                a4=(u2[ix,iz+1]-2*u2[ix][iz]+u2[ix][iz-1])/(dz * dz); ★
                a5=f[ix][iz]*w[it]/dx/dz;
                u3[ix][iz]=a1+a2*(a3+a4+a5);
            }
        }

    (이하생략)

}


추가적으로 포트란90의 경우에는, 단순하게
        u1=u2; u2=u3; u3=0.
의 코드로 배열간의 값 이동과 초기화가 가능한데, C++의 경우 아래와 같이 반복문을 통한 방법 외에
배열간의 값 이동 및 초기화를 빠른 속도로 해낼 수 있는 방법이 있는지 궁금합니다.

        for (ix=1;ix<=nx;ix++){
            for (iz=1;iz<=nz;iz++){
                u1[ix][iz]=u2[ix][iz];
                u2[ix][iz]=u3[ix][iz];
                u3[ix][iz]=0;
            }
        }





ps. 혹시 CUDA C에 대해 공부하시는 분들 계신지 궁금합니다...!!