Function norm4

examples/pde/bhrt.cpp:561–563

}

af::array norm4(const af::array& pos, const af::array& vel) {
    return dot_product(pos, vel, vel);
}

af::array partials(const af::array& pos4, uint32_t index, double rel_diff,
                   double abs_diff) {