void DoInitialBlockOfFlowsAllBeadFit(WorkerInfoQueueItem &item)
{
//printf("=====> All bead fit job on CPU\n");
BkgModelWorkInfo *info = (BkgModelWorkInfo *) (item.private_data);
FlowBlockSequence::const_iterator flowBlock =
info->inception_state->bkg_control.signal_chunks.flow_block_sequence.BlockAtFlow( info->flow );
info->bkgObj->FitAllBeadsForInitialFlowBlock( KEY_LEN, flowBlock->size(), info->table, flowBlock->begin() );
info->type = INITIAL_FLOW_BLOCK_REMAIN_REGIONAL_FIT;
info->pq->GetCpuQueue()->PutItem(item);
}
void BkgFitterTracker::DumpBkgModelBeadInfo (char *results_folder, int flow, bool last_flow,
bool debug_bead_only, FlowBlockSequence::const_iterator flow_block) const
{
// get some regional data for the entire chip as debug
// should be triggered by bkgmodel
if ( last_flow || flow + 1 == flow_block->end() )
{
DumpBkgModelBeadParams (results_folder, flow, debug_bead_only, flow_block->size() );
DumpBkgModelBeadOffset (results_folder, flow, debug_bead_only);
}
}
void DoPostFitProcessing(WorkerInfoQueueItem &item) {
//printf("=====> Post Processing job on CPU\n");
BkgModelWorkInfo *info = (BkgModelWorkInfo *) (item.private_data);
FlowBlockSequence::const_iterator flowBlock =
info->inception_state->bkg_control.signal_chunks.flow_block_sequence.BlockAtFlow( info->flow );
bool ewscale_correct = info->img->isEmptyWellAmplitudeAvailable();
int flow_block_id = info->inception_state->bkg_control.signal_chunks.flow_block_sequence.FlowBlockIndex( info->flow );
info->bkgObj->PreWellCorrectionFactors( ewscale_correct, flowBlock->size(), flowBlock->begin() ); // correct data for anything needed
// export to wells,debug, etc, reset for new set of traces
info->bkgObj->ExportAllAndReset(info->flow, info->last, flowBlock->size(), info->polyclonal_filter_opts, flow_block_id, flowBlock->begin() );
}
void DoInitialBlockOfFlowsRemainingRegionalFit(WorkerInfoQueueItem &item)
{
//printf("=====> Remaining fit steps job on CPU\n");
BkgModelWorkInfo *info = (BkgModelWorkInfo *) (item.private_data);
FlowBlockSequence::const_iterator flowBlock =
info->inception_state->bkg_control.signal_chunks.flow_block_sequence.BlockAtFlow( info->flow );
info->bkgObj->RemainingFitStepsForInitialFlowBlock( KEY_LEN, flowBlock->size(), info->table, flowBlock->begin() );
info->type = SINGLE_FLOW_FIT;
if (info->pq->GetGpuQueue() && info->pq->performGpuSingleFlowFitting())
info->pq->GetGpuQueue()->PutItem(item);
else
info->pq->GetCpuQueue()->PutItem(item);
}
void BkgFitterTracker::DumpBkgModelRegionInfo (char *results_folder, int flow, bool last_flow,
FlowBlockSequence::const_iterator flow_block) const
{
// get some regional data for the entire chip as debug
// should be triggered by bkgmodel
if ( last_flow || flow + 1 == flow_block->end() )
{
DumpBkgModelRegionParameters (results_folder, flow, flow_block->size() );
DumpBkgModelDarkMatter (results_folder, flow);
DumpBkgModelEmphasisTiming (results_folder, flow);
DumpBkgModelEmptyTrace (results_folder, flow, flow_block->size() );
DumpBkgModelInitVals (results_folder, flow);
}
}
void DoMultiFlowRegionalFit (WorkerInfoQueueItem &item) {
BkgModelWorkInfo *info = (BkgModelWorkInfo *) (item.private_data);
FlowBlockSequence::const_iterator flowBlock =
info->inception_state->bkg_control.signal_chunks.flow_block_sequence.BlockAtFlow( info->flow );
info->bkgObj->InitializeFlowBlock( flowBlock->size() );
info->bkgObj->ProcessImage (info->img, info->flow, info->flow - flowBlock->begin(),
flowBlock->size() );
// execute block if necessary
if (info->bkgObj->TestAndTriggerComputation (info->last))
{
info->bkgObj->MultiFlowRegionalFitting (info->flow, info->last,
info->flow_key, flowBlock->size(), info->table, flowBlock->begin() ); // CPU based regional fit
if ( info->inception_state->bkg_control.signal_chunks.flow_block_sequence.
HasFlowInFirstFlowBlock( info->flow ))
{
info->type = INITIAL_FLOW_BLOCK_ALLBEAD_FIT;
//if (info->pq->GetGpuQueue() && info->pq->performGpuMultiFlowFitting())
// info->pq->GetGpuQueue()->PutItem(item);
//else
// info->pq->GetCpuQueue()->PutItem(item);
info->QueueControl->AssignMultiFLowFitItemToQueue(item);
}
else
{
info->type = SINGLE_FLOW_FIT;
//if (info->pq->GetGpuQueue() && info->pq->performGpuSingleFlowFitting())
// info->pq->GetGpuQueue()->PutItem(item);
//else
// info->pq->GetCpuQueue()->PutItem(item);
info->QueueControl->AssignSingleFLowFitItemToQueue(item);
}
}
}
void BkgParamH5::IncrementalWrite ( int flow, bool last_flow, FlowBlockSequence::const_iterator flow_block,
int flow_block_id )
{
// single check: are we done with a compute block?
// try to isolate logic as much as possible - we should only see one check for writing, ever
if ( last_flow || flow == flow_block->end() - 1 )
{
IncrementalWriteBeads ( flow, flow_block_id );
IncrementalWriteRegions ( flow, flow_block_id );
IncrementalWriteBestRegion ( flow, last_flow );
IncrementalWriteRegionSamples ( flow, last_flow );
IncrementalWrite_xyflow ( last_flow );
}
}
