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);

    PDstandardNth2gt11  = kzero;
    PDstandardNth22gt11 = kifthen(origin, PDstandardNth11gt22,
                                          ksub(kmul(PDstandardNth1gt11, xinv), kmul(ktwo, kmul(ksub(gt11L, gt22L), x2inv))));
    PDstandardNth12gt11 = kzero;
    PDstandardNth23gt11 = kzero;
    PDstandardNth2gt12  = kifthen(origin, ksub(PDstandardNth1gt11, PDstandardNth1gt22),
                                          kmul(ksub(gt11L, gt22L), xinv));
    PDstandardNth22gt12 = kzero;
    PDstandardNth12gt12 = kifthen(origin, kmul(ksub(PDstandardNth11gt11, PDstandardNth11gt22), khalf),
                                          ksub(kmul(ksub(PDstandardNth1gt11, PDstandardNth1gt22), xinv), kmul(ksub(gt11L, gt22L), x2inv)));
    PDstandardNth23gt12 = kifthen(origin, ksub(PDstandardNth13gt11, PDstandardNth13gt22),
                                          kmul(ksub(PDstandardNth3gt11, PDstandardNth3gt22), xinv));
    PDstandardNth2gt13  = kzero;
    PDstandardNth22gt13 = kifthen(origin, kmul(PDstandardNth11gt13, khalf),
                                          ksub(kmul(PDstandardNth1gt13, xinv), kmul(gt13L, x2inv)));
    PDstandardNth12gt13 = kzero;
    PDstandardNth23gt13 = kzero;
    PDstandardNth2gt22  = kzero;
    PDstandardNth22gt22 = kifthen(origin, PDstandardNth11gt11,
                                          kadd(kmul(PDstandardNth1gt22, xinv), kmul(ktwo, kmul(ksub(gt11L, gt22L), x2inv))));
    PDstandardNth12gt22 = kzero;
    PDstandardNth23gt22 = kzero;
    PDstandardNth2gt23  = kifthen(origin, PDstandardNth1gt13, kmul(gt13L, xinv));
    PDstandardNth22gt23 = kzero;
    PDstandardNth12gt23 = kifthen(origin, kmul(PDstandardNth11gt13, khalf),
                                          ksub(kmul(PDstandardNth1gt13, xinv), kmul(gt13L, x2inv)));
    PDstandardNth23gt23 = kifthen(origin, PDstandardNth13gt13, kmul(PDstandardNth3gt13, xinv));
    PDstandardNth2gt33  = kzero;
    PDstandardNth22gt33 = kifthen(origin, PDstandardNth11gt33, kmul(PDstandardNth1gt33, xinv));
    PDstandardNth12gt33 = kzero;
    PDstandardNth23gt33 = kzero;

    PDstandardNth2phi  = kzero;
    PDstandardNth22phi = kifthen(origin, PDstandardNth11phi, kmul(PDstandardNth1phi, xinv));
    PDstandardNth12phi = kzero;
    PDstandardNth23phi = kzero;
    PDstandardNth2Xt1  = kzero;
    PDstandardNth2Xt2  = kifthen(origin, PDstandardNth1Xt1, kmul(Xt1L, xinv));
    PDstandardNth2Xt3  = kzero;
    PDstandardNth2trK  = kzero;

    //PDstandardNth2At11 = origin ? -2 * PDstandardNth1At12 : -2 * At12[index] / xx;
    //PDstandardNth2At13 = origin ? -PDstandardNth1At23 : - At23[index] / xx;
    //PDstandardNth2At22 = origin ? 2 * PDstandardNth1At12 : 2 * At12[index] / xx;
    PDstandardNth2At11 = kzero;
    PDstandardNth2At12 = kifthen(origin, ksub(PDstandardNth1At11, PDstandardNth1At22),
                                         kmul(ksub(At11L, At22L), xinv));
    PDstandardNth2At13 = kzero;
    PDstandardNth2At22 = kzero;
    PDstandardNth2At23 = kifthen(origin, PDstandardNth1At13, kmul(At13L, xinv));
    PDstandardNth2At33 = kzero;