campreilly/UnderSeaModelingLibrary — Issue #2: Investigate use of ViennaCL to implement GPU implementation
Author: campreilly
Created: 2012-03-06 19:17:14+00:00
State: open
Labels: enhancement

See http://viennacl.sourceforge.net/ for code and docs.

* Need to balance slow operations for individual accesses with fast matrix operations.
* Several of the intermediate vars in wave_element could be GPU based. Because they are only used in matrix form. Could be a serious benefit to wave_element::update(). Examples: _c2_r, _dc_c, _sin_theta, _cot_theta, pos_gradient?, ndir_gradient?, sound_speed?, sound_gradient?.
* The operations in ode_integ::ab3_pos() and ode_integ::ab3_ndir() could also run very quickly iff the constants were also GPU scalars.
* Would it make sense for GPU version of wvectors to be split into separate rho, theta, and phi components?