Function sph_to_cart_position

examples/pde/bhrt.cpp:405–417

* @brief Transform the position vectors from spherical to cartesian coordinates * * @param pos * @return af::array */

 * @return af::array
 */
af::array sph_to_cart_position(const af::array& pos) {
    af::array r = pos(0, af::span);
    af::array o = pos(1, af::span);
    af::array p = pos(2, af::span);

    af::array x = r * af::sin(o) * af::cos(p);
    af::array y = r * af::sin(o) * af::sin(p);
    af::array z = r * af::cos(o);

    af::array transformed_pos = af::join(0, x, y, z);

    return transformed_pos;
}

/**
 * @brief Computes the inverse of a 4x4 matrix with the layout