#include "interpolate.h" 
#include "CVDButil.h" 

float lerpConstant( constant float * in, int size, float pos);

kernel void kernelName( 
                 int P_length, 
                 global float * P ,
 global void * vdb_dens 
)
{
    int idx = get_global_id(0);
    if (idx >= P_length)
        return;
    const global cnanovdb_griddata *grid = (const global cnanovdb_griddata *) vdb_dens;
    cnanovdb_readaccessor acc;
    //check for vdb type
    if(cnanovdb_griddata_validF(grid)) {
        cnanovdb_getreadaccessor(grid, &acc);
        float3 sampleP;
        sampleP.x = P[(idx*3) +0];
        sampleP.y = P[(idx*3) +1];
        sampleP.z = P[(idx*3) +2];
        float dens = cnanovdb_sampleF(grid, &acc, sampleP);
        P[(idx*3) +1] += dens;
    }
}