if (_ldk.getStencil()!=LocalDiffusionKernel.Stencil.D21) { for (int i2=n2-1; i2>0; --i2) { for (int i1=n1-1; i1>0; --i1) { if (_ldk.getStencil()!=LocalDiffusionKernel.Stencil.D21) lsf.applySmoothS(r,r); lsf.apply(_tensors,_c,s,r,q);
if (_ldk.getStencil()!=LocalDiffusionKernel.Stencil.D21) { for (int i3=n3-1; i3>0; --i3) { for (int i2=n2-1; i2>0; --i2) { if (_ldk.getStencil()!=LocalDiffusionKernel.Stencil.D21) lsf.applySmoothS(r,r); lsf.apply(_tensors,_c,s,r,q);