Vectorizing the forward mode of ADOL-C on a GPU ... - Autodiff.org

. K. Kulshreshtha, A. Koniaeva 9 / 13 <strong>Vectorizing</strong> <strong>ADOL</strong>-C using CUDA Euro AD 10.06.2013Implementation.Implementation__global__ void kernel_func(double* inx, size_t n,double* outy, size_t m, double* dery);__device__ void evalf(adtlc::adouble* x, size_t n,adtlc::adouble* y, size_t m);• Say derivatives are evaluated at N = (N x × N y × N z ) points inp = (p x × p y × p z ) directions simultaneously• Use Nb = {N x , N y , N z } blocks <strong>of</strong> Nt = {p x , p y , p z } threads each• inx has n × N entries• outy has m × N entries• dery has m × p × N entries