/*        Generated by TAPENADE     (INRIA, Tropics team)
    Tapenade 3.7 (r4786) - 21 Feb 2013 15:53
*/
#include "cstd.h"
extern void pushcontrol1b(int cc);
extern void popcontrol1b(int *cc);

/*
  Differentiation of acd_2d_8 in reverse (adjoint) mode:
   gradient     of useful results: **uc **up
   with respect to varying inputs: **csq **uc **up
   RW status of diff variables: **csq:out **uc:incr **up:in-out
   Plus diff mem management of: csq:in *csq:in uc:in *uc:in up:in
                *up:in
*/
void acd_2d_8_b(float **uc, float **ucb, float **up, float **upb, float **csq,
        float **csqb, int *s, int *e, float c0, float *c1, float *c2, float *
        c3, float *c4, int *lbc, int *rbc) {
    int i0, i1;
    int branch;
    float tmp4b;
    float tmpb;
    float tmp0b;
    float tmp3b;
    float tempb;
    float tmp2b;
    float tmp1b;
    /* boundary conditions - note that uc[-1][i]=0 etc. */
    if (lbc[1])
        pushcontrol1b(0);
    else
        pushcontrol1b(1);
    if (rbc[1])
        pushcontrol1b(0);
    else
        pushcontrol1b(1);
    if (lbc[0])
        pushcontrol1b(0);
    else
        pushcontrol1b(1);
    if (rbc[0])
        for (i1 = e[1]; i1 > s[1]-1; --i1) {
            upb[i1][e[0] - 2] = upb[i1][e[0] - 2] - upb[i1][e[0] + 4];
            upb[i1][e[0] + 4] = 0.0;
            upb[i1][e[0] - 1] = upb[i1][e[0] - 1] - upb[i1][e[0] + 3];
            upb[i1][e[0] + 3] = 0.0;
            upb[i1][e[0] - 0] = upb[i1][e[0] - 0] - upb[i1][e[0] + 2];
            upb[i1][e[0] + 2] = 0.0;
        }
    popcontrol1b(&branch);
    if (branch == 0)
        for (i1 = e[1]; i1 > s[1]-1; --i1) {
            upb[i1][s[0] + 2] = upb[i1][s[0] + 2] - upb[i1][s[0] - 4];
            upb[i1][s[0] - 4] = 0.0;
            upb[i1][s[0] + 1] = upb[i1][s[0] + 1] - upb[i1][s[0] - 3];
            upb[i1][s[0] - 3] = 0.0;
            upb[i1][s[0] + 0] = upb[i1][s[0] + 0] - upb[i1][s[0] - 2];
            upb[i1][s[0] - 2] = 0.0;
        }
    popcontrol1b(&branch);
    if (branch == 0)
        for (i0 = e[0]; i0 > s[0]-1; --i0) {
            tmp4b = upb[e[1] + 4][i0];
            upb[e[1] + 4][i0] = 0.0;
            upb[e[1] - 2][i0] = upb[e[1] - 2][i0] - tmp4b;
            tmp3b = upb[e[1] + 3][i0];
            upb[e[1] + 3][i0] = 0.0;
            upb[e[1] - 1][i0] = upb[e[1] - 1][i0] - tmp3b;
            tmp2b = upb[e[1] + 2][i0];
            upb[e[1] + 2][i0] = 0.0;
            upb[e[1] - 0][i0] = upb[e[1] - 0][i0] - tmp2b;
        }
    popcontrol1b(&branch);
    if (branch == 0)
        for (i0 = e[0]; i0 > s[0]-1; --i0) {
            tmp1b = upb[s[1] - 4][i0];
            upb[s[1] - 4][i0] = 0.0;
            upb[s[1] + 2][i0] = upb[s[1] + 2][i0] - tmp1b;
            tmp0b = upb[s[1] - 3][i0];
            upb[s[1] - 3][i0] = 0.0;
            upb[s[1] + 1][i0] = upb[s[1] + 1][i0] - tmp0b;
            tmpb = upb[s[1] - 2][i0];
            upb[s[1] - 2][i0] = 0.0;
            upb[s[1] + 0][i0] = upb[s[1] + 0][i0] - tmpb;
        }
/*    **csqb = 0.0; */
    for (i1 = e[1]; i1 > s[1]-1; --i1)
        for (i0 = e[0]; i0 > s[0]-1; --i0) {
            float lap = c0*uc[i1][i0]	+
            c1[0]*(uc[i1][i0+1] + uc[i1][i0-1]) +
            c1[1]*(uc[i1+1][i0] + uc[i1-1][i0]) +
            c2[0]*(uc[i1][i0+2] + uc[i1][i0-2]) +
            c2[1]*(uc[i1+2][i0] + uc[i1-2][i0]) +
            c3[0]*(uc[i1][i0+3] + uc[i1][i0-3]) +
            c3[1]*(uc[i1+3][i0] + uc[i1-3][i0]) +
            c4[0]*(uc[i1][i0+4] + uc[i1][i0-4]) +
            c4[1]*(uc[i1+4][i0] + uc[i1-4][i0]);
            
            tempb = csq[i1][i0]*upb[i1][i0];
            ucb[i1][i0] = ucb[i1][i0] + c0*tempb + 2.0*upb[i1][i0];
            csqb[i1][i0] = csqb[i1][i0] + lap*upb[i1][i0];
            ucb[i1][i0 + 1] = ucb[i1][i0 + 1] + c1[0]*tempb;
            ucb[i1][i0 - 1] = ucb[i1][i0 - 1] + c1[0]*tempb;
            ucb[i1 + 1][i0] = ucb[i1 + 1][i0] + c1[1]*tempb;
            ucb[i1 - 1][i0] = ucb[i1 - 1][i0] + c1[1]*tempb;
            ucb[i1][i0 + 2] = ucb[i1][i0 + 2] + c2[0]*tempb;
            ucb[i1][i0 - 2] = ucb[i1][i0 - 2] + c2[0]*tempb;
            ucb[i1 + 2][i0] = ucb[i1 + 2][i0] + c2[1]*tempb;
            ucb[i1 - 2][i0] = ucb[i1 - 2][i0] + c2[1]*tempb;
            ucb[i1][i0 + 3] = ucb[i1][i0 + 3] + c3[0]*tempb;
            ucb[i1][i0 - 3] = ucb[i1][i0 - 3] + c3[0]*tempb;
            ucb[i1 + 3][i0] = ucb[i1 + 3][i0] + c3[1]*tempb;
            ucb[i1 - 3][i0] = ucb[i1 - 3][i0] + c3[1]*tempb;
            ucb[i1][i0 + 4] = ucb[i1][i0 + 4] + c4[0]*tempb;
            ucb[i1][i0 - 4] = ucb[i1][i0 - 4] + c4[0]*tempb;
            ucb[i1 + 4][i0] = ucb[i1 + 4][i0] + c4[1]*tempb;
            ucb[i1 - 4][i0] = ucb[i1 - 4][i0] + c4[1]*tempb;
            upb[i1][i0] = -upb[i1][i0];
        }
}