CCTK_REAL_VEC kzero = ToReal(0.0);
    CCTK_REAL_VEC kone = ToReal(1.0);
    CCTK_REAL_VEC ktwo = ToReal(2.0);
    CCTK_REAL_VEC xx = vec_load(x[index]);
    CCTK_REAL_VEC x2 = kmul(xx, xx);

    CCTK_REAL_VEC xinv  = kdiv(kone, xx);
    CCTK_REAL_VEC x2inv = kdiv(kone, x2);

    CCTK_REAL_VEC absx = kfabs(xx);
    CCTK_REAL_VEC eps  = ToReal(1e-8);
    CCTK_BOOLEAN_VEC origin = kcmplt(absx, eps);

    PDupwindNthAnti2At11 = kzero;
    PDupwindNthSymm2At11 = kzero;
    PDupwindNthAnti2At12 = kifthen(origin, ksub(PDupwindNthAnti1At11, PDupwindNthAntiAt22), kmul(ksub(At11L, At22L), xinv));
    PDupwindNthSymm2At12 = kifthen(origin, ksub(PDupwindNthSymm1At11, PDupwindNthSymmAt22), kmul(ksub(At11L, At22L), xinv));
    PDupwindNthAnti2At13 = kzero;
    PDupwindNthSymm2At13 = kzero;
    PDupwindNthAnti2At22 = kzero;
    PDupwindNthSymm2At22 = kzero;
    PDupwindNthAnti2At23 = kifthen(origin, PDupwindNth1AntiAt13, kmul(At13L, xinv));
    PDupwindNthSymm2At23 = kifthen(origin, PDupwindNth1SymmAt13, kmul(At13L, xinv));
    PDupwindNthAnti2At33 = kzero;
    PDupwindNthSymm2At33 = kzero;

    PDupwindNthAnti2B1 = kzero;
    PDupwindNthSymm2B1 = kzero;
    PDupwindNthAnti2B2 = kifthen(origin, PDstandardNthAnti1B1, kmul(B1L, xinv));
    PDupwindNthSymm2B2 = kifthen(origin, PDstandardNthSymm1B1, kmul(B1L, xinv));
    PDupwindNthAnti2B3 = kzero;
    PDupwindNthSymm2B3 = kzero;

    PDupwindNthAnti2beta1 = kzero;
    PDupwindNthSymm2beta1 = kzero;
    PDupwindNthAnti2beta2 = kifthen(origin, PDstandardNthAnti1beta1, kmul(beta1L, xinv));
    PDupwindNthSymm2beta2 = kifthen(origin, PDstandardNthSymm1beta1, kmul(beta1L, xinv));
    PDupwindNthAnti2beta3 = kzero;
    PDupwindNthSymm2beta3 = kzero;

    PDupwindNthAnti2gt11 = kzero;
    PDupwindNthSymm2gt11 = kzero;
    PDupwindNthAnti2gt12 = kifthen(origin, ksub(PDstandardNthAnti1gt11, PDstandardNthAnti1gt22), kmul(ksub(gt11L, gt22L), xinv));
    PDupwindNthSymm2gt12 = kifthen(origin, ksub(PDstandardNthSymm1gt11, PDstandardNthSymm1gt22), kmul(ksub(gt11L, gt22L), xinv));
    PDupwindNthAnti2gt13 = kzero;
    PDupwindNthSymm2gt13 = kzero;
    PDupwindNthAnti2gt22 = kzero;
    PDupwindNthSymm2gt22 = kzero;
    PDupwindNthAnti2gt23 = kifthen(origin, PDstandardNthAnti1gt13, kmul(gt13L, xinv));
    PDupwindNthSymm2gt23 = kifthen(origin, PDstandardNthSymm1gt13, kmul(gt13L, xinv));
    PDupwindNthAnti2gt33 = kzero;
    PDupwindNthSymm2gt33 = kzero;

    PDupwindNthAnti2phi = kzero;
    PDupwindNthSymm2phi = kzero;
    PDupwindNthAnti2trK = kzero;
    PDupwindNthSymm2trK = kzero;
    PDupwindNthAnti2Xt2 = kifthen(origin, PDstandardNthAnti1Xt1, kmul(Xt1L, xinv));
    PDupwindNthSymm2Xt2 = kifthen(origin, PDstandardNthSymm1Xt1, kmul(Xt1L, xinv));