Merge pull request #337 from zhujun98/dev

Avoid destructing Reconstructor during reconstructing

docs/docs/installation.md

@@ -32,6 +32,8 @@ make -j12 && make install
 
 - Windows
 
+TBD
+
 #### Installing the GUI client
 
 - Linux / MacOS
@@ -51,6 +53,8 @@ make -j12 && make install
 
 - Windows
 
+TBD
+
 ## Package managers
 
 TBD

docs/docs/performance_consideration.md

@@ -4,6 +4,11 @@ Currently, the bottleneck of the data processing pipeline is the preprocessing p
 is still implemented in CPU. The total number of threads of the current benchmark machine is 
 only 32. We expect a significant performance boost on a GPU node with a decent number of threads.
 
+The 3 GB/s throughput was measured using a projection size of 512 x 512 pixels and 500 projections per scan.
+The throughput will decrease as the total size of data per scan increases since the performance does
+not scale linearly. If you need a high throughput with large projections sizes and/or large numbers of projections per
+scan, it is recommended to down-sample the projections, which can be configured in the GUI.
+
 ## Visualization
 
 The data rate in RECASTX is really high not only in terms of the raw 

recon/src/application.cpp

@@ -528,9 +528,6 @@ void Application::initParams() {
     group_size_ = (scan_mode_ == rpc::ScanMode_Mode_CONTINUOUS || angle_count_ == 0)
             ? scan_update_interval_ : angle_count_;
     proj_mediator_->setFilter(group_size_);
-
-    sino_initialized_ = false;
-    gpu_buffer_index_ = 0;
 }
 
 void Application::maybeInitFlatFieldBuffer(uint32_t row_count, uint32_t col_count) {
@@ -601,6 +598,13 @@ void Application::initReconstructor(uint32_t col_count, uint32_t row_count) {
     };
 
     double_buffering_ = scan_mode_ == rpc::ScanMode_Mode_DYNAMIC;
+
+    std::lock_guard lck(recon_mtx_);
+
+    sino_uploaded_ = false;
+    sino_initialized_ = false;
+    gpu_buffer_index_ = 0;
+
     recon_.reset();
     recon_ = recon_factory_->create(proj_geom, slice_geom, volume_geom, double_buffering_);
 

recon/src/reconstructor.cpp

@@ -135,7 +135,10 @@ void ParallelBeamReconstructor::reconstructSlice(Orientation x, int buffer_idx,
         py = scale.cwiseProduct(rot * py);
         d -= 0.5f * (num_cols * px + num_rows * py);
 
-        vec_buf_[i] = {r[0],  r[1],  r[2],  d[0],  d[1],  d[2], px[0], px[1], px[2], py[0], py[1], py[2]};
+        vec_buf_[i] = {r[0],   r[1],  r[2],
+                       d[0],   d[1],  d[2],
+                       px[0], px[1], px[2],
+                       py[0], py[1], py[2]};
         ++i;
     }
 
@@ -234,7 +237,10 @@ void ConeBeamReconstructor::reconstructSlice(Orientation x, int buffer_idx, Tens
         t1 = scale.cwiseProduct(rot * t1);
         t2 = scale.cwiseProduct(rot * t2);
 
-        vec_buf_[i] = {s[0],  s[1],  s[2],  d[0],  d[1],  d[2], t1[0], t1[1], t1[2], t2[0], t2[1], t2[2]};
+        vec_buf_[i] = {s[0],   s[1],  s[2],
+                       d[0],   d[1],  d[2],
+                       t1[0], t1[1], t1[2],
+                       t2[0], t2[1], t2[2]};
         ++i;
     }