if (ts < tns):
                                        # exploiting the symmetry
                                        # *CC
                                        pdf[i, 1:Nx1, tns] += fxind[IJ:J] * dt
                                        err[i, 1:Nx1, tns] += err0[IJ:J] * dt2