OpenCL產生不正確的計算

我一直在嘗試使用openCL做一些計算，但結果不正確。OpenCL產生不正確的計算

我路輸入三個FLOAT3的看起來像這樣：

[300000,0,0] 
[300000,300000,0] 
[300000,300000,300000]

這個內核：

__kernel void gravitate(__global const float3 *position,__global const float3 *momentum,__global const float3 *mass,__global float3 *newPosition,__global float3 *newMomentum,unsigned int numBodies,unsigned int seconds) 
{ 
    int gid=get_global_id(0); 

    newPosition[gid]=position[gid]*2; 
    newMomentum[gid]=momentum[gid]*2; 
}

我本來試圖模擬一組機構的引力相互作用這就是爲什麼內核看起來像那樣。

對於輸出中，我得到：

[600000.0,0.0,0.0] 
[x,600000.0,0.0] 
[600000.0,x,600000.0]

其中x是在NAN範圍爲10^-44介於0和每內核運行時是不同的值。 newPosition和newMomentum都有相同的錯誤輸出模式。

，我使用看起來像這樣的Python代碼：

def __init__(self): 
     self.context=cl.create_some_context() 
     self.queue=cl.CommandQueue(self.context) 

     f=open("gravitate.cl") 
     self.program=cl.Program(self.context,f.read()).build() 

def simulate(self,seconds): 
     bodyPosition=[] 
     bodyMomentum=[] 
     bodyMass=[] 

     for body in self.objects: 
      bodyPosition+=[body.position.x,body.position.y,body.position.z] 
      bodyMomentum+=[body.momentum.x,body.momentum.y,body.momentum.z] 
      bodyMass+=[body.mass] 

     bodyPosition=numpy.array(bodyPosition).astype(numpy.float32) 
     bodyMomentum=numpy.array(bodyMomentum).astype(numpy.float32) 
     bodyMass=numpy.array(bodyMass).astype(numpy.float32) 

     bodyPositionCl=cl.Buffer(self.context,cl.mem_flags.READ_ONLY|cl.mem_flags.COPY_HOST_PTR,hostbuf=bodyPosition) 
     bodyMomentumCl=cl.Buffer(self.context,cl.mem_flags.READ_ONLY|cl.mem_flags.COPY_HOST_PTR,hostbuf=bodyMomentum) 
     bodyMassCl=cl.Buffer(self.context,cl.mem_flags.READ_ONLY|cl.mem_flags.COPY_HOST_PTR,hostbuf=bodyMass) 

     newBodyPosition=numpy.zeros(len(self.objects)*3).astype(numpy.float32) 
     newBodyMomentum=numpy.zeros(len(self.objects)*3).astype(numpy.float32) 

     newBodyPositionCl=cl.Buffer(self.context,cl.mem_flags.WRITE_ONLY,newBodyPosition.nbytes) 
     newBodyMomentumCl=cl.Buffer(self.context,cl.mem_flags.WRITE_ONLY,newBodyMomentum.nbytes) 

     self.program.gravitate(self.queue,(3,),None,bodyPositionCl,bodyMomentumCl,bodyMassCl,newBodyPositionCl,newBodyMomentumCl,numpy.uint32(len(self.objects)),numpy.uint32(seconds)) 
     cl.enqueue_read_buffer(self.queue,newBodyPositionCl,newBodyPosition).wait() 
     cl.enqueue_read_buffer(self.queue,newBodyMomentumCl,newBodyMomentum).wait()

來源

2011-07-06 user467526

FLOAT3是16字節對齊。請參閱OpenCL 1.1規範6.1.5。

來源

2011-07-06 07:37:23

OpenCL產生不正確的計算

回答

相關問題