velocity_potential.cc

#include <vtestbed/config/opencl.hh>
#include <vtestbed/config/real_type.hh>
#include <vtestbed/core/derivative.hh>
#include <vtestbed/core/math.hh>
#include <vtestbed/core/velocity_potential.hh>
#include <vtestbed/core/wave_numbers.hh>

using vtb::core::Policy;

template <class T>
void
vtb::core::Linear_velocity_potential_solver_base<T>::second_function(
    const Grid<T,4>& grid_tzxy,
    const Array3<T>& zeta,
    Array3<C>& result
) {
    if (this->assumptions().small_amplitude()) {
        result = -derivative<T,3>(zeta, grid_tzxy.delta(0), 0);
    } else {
        using blitz::Range;
        using blitz::sqrt;
        using blitz::pow;
        Array3<T> zeta_t(derivative<T,3>(zeta, grid_tzxy.delta(0), 0));
        Array3<T> zeta_x(derivative<T,3>(zeta, grid_tzxy.delta(2), 1));
        Array3<T> zeta_y(derivative<T,3>(zeta, grid_tzxy.delta(3), 2));
        Array3<T> sqrt_zeta(sqrt(T{1} + pow2(zeta_x) + pow2(zeta_y)));
        const C i(0,1);
        result = zeta_t /
            (i*((zeta_x + zeta_y)/sqrt_zeta - zeta_x - zeta_y) -
             T{1}/sqrt_zeta);
    }
}


template class
vtb::core::Velocity_potential_solver<VTB_REAL_TYPE,3>;

template class
vtb::core::Linear_velocity_potential_solver_base<VTB_REAL_TYPE>;