/*        Generated by TAPENADE     (INRIA, Ecuador team)
    Tapenade 3.10 (r5717) - 30 Jul 2015 16:03
*/
#include "cstd.h"

/*
  Differentiation of acd_2d in forward (tangent) mode:
   variations   of useful results: **up
   with respect to varying inputs: **csq **uc **up
   RW status of diff variables: **csq:in **uc:in **up:in-out
   Plus diff mem management of: lap:in csq:in *csq:in uc:in *uc:in
                up:in *up:in
*/
void acd_2d_d(float **uc, float **ucd, float **up, float **upd, float **csq, 
        float **csqd, int *s, int *e, float **c, int k, int *lbc, int *rbc, 
        float *lap, float *lapd) 
{
    int i0, i1, ioff;
    int s0 = s[0];
    int e0 = e[0];
    float c0 = c[0][0] + c[1][0];

    for (i1 = s[1]; i1 < e[1]+1; ++i1) {
        // offset independent part of stencil
        //#pragma ivdep
        for (i0 = s0; i0 < e0+1; ++i0) {
            upd[i1][i0] = 2.0*ucd[i1][i0] - upd[i1][i0];
            up[i1][i0] = 2.0*uc[i1][i0] - up[i1][i0];
            lapd[i0] = c0*ucd[i1][i0];
            lap[i0] = c0*uc[i1][i0];
        }
        // offset dependent part of stencil
        for (ioff = 1; ioff < k+1; ++ioff)
            //#pragma ivdep
            for (i0 = s0; i0 < e0+1; ++i0) {
                lapd[i0] = lapd[i0] + c[0][ioff]*(ucd[i1][i0+ioff]+ucd[i1][i0-
                    ioff]) + c[1][ioff]*(ucd[i1+ioff][i0]+ucd[i1-ioff][i0]);
                lap[i0] += c[0][ioff]*(uc[i1][i0+ioff]+uc[i1][i0-ioff]) + c[1]
                [ioff]*(uc[i1+ioff][i0]+uc[i1-ioff][i0]);
            }
        for (i0 = s0; i0 < e0+1; ++i0) {
            upd[i1][i0] = upd[i1][i0] + csqd[i1][i0]*lap[i0] + csq[i1][i0]*
                lapd[i0];
            up[i1][i0] += csq[i1][i0]*lap[i0];
        }
    }

    /* boundary conditions - note that uc[-1][i]=0 etc. */
    if (lbc[1])
        for (ioff = 0; ioff < k-1; ++ioff)
            //#pragma ivdep
            for (i0 = s0; i0 < e0+1; ++i0) {
                upd[s[1] - ioff - 2][i0] = -upd[s[1]+ioff][i0];
                up[s[1] - ioff - 2][i0] = -up[s[1]+ioff][i0];
            }
    if (rbc[1])
        for (ioff = 0; ioff < k-1; ++ioff)
            //#pragma ivdep
            for (i0 = s0; i0 < e0+1; ++i0) {
                upd[e[1] + 2 + ioff][i0] = -upd[e[1]-ioff][i0];
                up[e[1] + 2 + ioff][i0] = -up[e[1]-ioff][i0];
            }
    if (lbc[0])
        for (i1 = s[1]; i1 < e[1]+1; ++i1)
            for (ioff = 0; ioff < k-1; ++ioff) {
                upd[i1][s[0] - ioff - 2] = -upd[i1][s[0]+ioff];
                up[i1][s[0] - ioff - 2] = -up[i1][s[0]+ioff];
            }
    if (rbc[0])
        for (i1 = s[1]; i1 < e[1]+1; ++i1)
            for (ioff = 0; ioff < k-1; ++ioff) {
                upd[i1][e[0] + ioff + 2] = -upd[i1][e[0]-ioff];
                up[i1][e[0] + ioff + 2] = -up[i1][e[0]-ioff];
            }
}