private final static void modFloat(int sa, final PyMultiarray a, int sb, final PyMultiarray b, int sr, final PyMultiarray r, final int d) {
        final int dsa = a.strides[d], dsb = b.strides[d], dsr = r.strides[d], maxSr = sr
                + r.dimensions[d] * r.strides[d];
        if (d == r.dimensions.length - 1) {
            final float[] aData = (float[]) a.data, bData = (float[]) b.data;
            final float[] rData = (float[]) r.data;
            for (; sr != maxSr; sa += dsa, sb += dsb, sr += dsr) {
                rData[sr] = (aData[sa] % bData[sb]);
            }
        } else {
            for (; sr != maxSr; sa += dsa, sb += dsb, sr += dsr) {
                PyMultiarray.modFloat(sa, a, sb, b, sr, r, d + 1);
            }
        }
    }