openmp/src/wave.cpp
std::complex<double> 
gradIncident(const double nvec[], const double kvec[], 
             const double point[]) {
    return J1*(nvec[0]*kvec[0] + nvec[1]*kvec[1])*incident(kvec, point);
}
/**
 * Normal gradient of the incident wave, assumes incident wave is exp(1j * kvec.x)
 *
 * @param nvec normal vector pointing inwards
 * @param kvec incident wave vector
 * @param point (source) point
 * @return amplitude
 */
Normal gradient of the incident wave, assumes incident wave is exp(1j * kvec.x)
@param nvec normal vector pointing inwards
@param kvec incident wave vector
@param point (source) point
@return amplitude
Normal gradient of the incident wave assumes incident wave is exp ( 1j * kvec . x ) @param nvec normal vector pointing inwards @param kvec incident wave vector @param point ( source ) point @return amplitude
std::complex<double> 
gradIncident(const double nvec[], const double kvec[], 
             const double point[]) {
    return J1*(nvec[0]*kvec[0] + nvec[1]*kvec[1])*incident(kvec, point);
}
std :: complex < double > gradIncident ( const double nvec [ ] , const double kvec [ ] , const double point [ ] ) { return J1 * ( nvec [ 0 ] * kvec [ 0 ] + nvec [ 1 ] * kvec [ 1 ] ) * incident ( kvec , point ) ; }
Normal gradient of the incident wave, assumes incident wave is exp(1j * kvec.x)
 @param nvec normal vector pointing inwards
 @param kvec incident wave vector
 @param point (source) point
 @return amplitude
Normal gradient of the incident wave assumes incident wave is exp ( 1j * kvec . x ) @param nvec normal vector pointing inwards @param kvec incident wave vector @param point ( source ) point @return amplitude