// Begin Code to Initialize and Use GPU
{
// Set the limits of the calculations to the interior points
    iter float2 it<insizey,insizex> = iter ( float2(0.0f, 0.0f), 
                float2((float)insizex, (float)insizey));

    // Initialize the streams
    float Obs<1,1>;
    float Hy<insizey,insizex>, Hx<insizey,insizex>;
    float o_Ez<insizey,insizex>, o_Hy<insizey,insizex>;
    float o_Hx<insizey,insizex>;
    float Ceze<insizey,insizex>, Chxh<insizey,insizex>;
    float Chyh<insizey,insizex>;
    float Cezh<insizey,insizex>, Chxe<insizey,insizex>;
    float Chye<insizey,insizex>, Cs<insizey,insizex>;
    float Ez<insizey,insizex>;
    float Kex<insizex>, Khx<insizex>;
    float Bex<insizex>, Bhx<insizex>;
    float Cex<insizex>, Chx<insizex>;
    float Key<insizey>, Khy<insizey>;
    float Bey<insizey>, Bhy<insizey>;
    float Cey<insizey>, Chy<insizey>;
    float psiEzx<insizey,insizex>, o_psiEzx<insizey,insizex>;
    float psiEzy<insizey,insizex>, o_psiEzy<insizey,insizex>;
    float psiHxy<insizey,insizex>, o_psiHxy<insizey,insizex>;
    float psiHyx<insizey,insizex>, o_psiHyx<insizey,insizex>;


    // Copy our C arrays into the Streams so the GPU 
    // can process them
    streamRead(Ez, aEz);
    streamRead(Hy, aHy);
    streamRead(Hx, aHx);
    streamRead(o_Ez, aEz);
    streamRead(o_Hy, aHy);
    streamRead(o_Hx, aHx);
    streamRead(Ceze, aCeze);
    streamRead(Chxh, aChxh);
    streamRead(Chyh, aChyh);
    streamRead(Cezh, aCezh);
    streamRead(Chxe, aChxe);
    streamRead(Chye, aChye);
    streamRead(Cs, aCs);

    streamRead(Kex, aKex);
    streamRead(Khx, aKhx);
    streamRead(Bex, abex);
    streamRead(Bhx, abhx);
    streamRead(Cex, acex);
    streamRead(Chx, achx);

    streamRead(Key, aKey);
    streamRead(Khy, aKhy);
    streamRead(Bey, abey);
    streamRead(Bhy, abhy);
    streamRead(Cey, acey);
    streamRead(Chy, achy);

    streamRead(psiEzx, apsi);
    streamRead(o_psiEzx, apsi);
    streamRead(psiEzy, apsi);
    streamRead(o_psiEzy, apsi);
    streamRead(psiHxy, apsi);
    streamRead(o_psiHxy, apsi);
    streamRead(psiHyx, apsi);
    streamRead(o_psiHyx, apsi);