1D trapezoidal integration.
Parameters: | x_gpu : pycuda.gpuarray.GPUArray
dx : scalar
handle : int
|
---|---|
Returns: | result : float
|
Examples
>>> import pycuda.autoinit
>>> import pycuda.gpuarray
>>> import numpy as np
>>> import integrate
>>> integrate.init()
>>> x = np.asarray(np.random.rand(10), np.float32)
>>> x_gpu = gpuarray.to_gpu(x)
>>> z = integrate.trapz(x_gpu)
>>> np.allclose(np.trapz(x), z)
True