Ben Hirsch edited CUDA_Implementation.tex  about 11 years ago

Commit id: 64f7d4de4c207d1c2a04a61bd4342b679ec4838c

deletions | additions      

       

\begin{quote}  int idx;  Vector Vector\  currenttDataLocus (X(0, i), numObs); Vector Vector\  currentTestLocus (test_X(0, i), testObs); float oldSamplei = solSample\[i\];  float rhsModeli = (currentDataLocus.dot(RInverseY)) + diagLhs(i) * oldSamplei;  rhsModeli *= invVarResidual;  float Lhs = diagLhs\[i\] diagLhs[i]  * invVarResidual + invVarEffects; float probDeltaOne = 1.0f / (1.0f + exp(logDeltaZero + 0.5f * (log(Lhs) + logVarEffects - rhsModeli*rhsModeli/Lhs) - logPiComp));  if (matvec::ranf() < probDeltaOne) \{