Add reconstructArt(unfinished)

src/transmissionReconstruction/reconstruction/reconstruction.cpp

@@ -4,6 +4,11 @@
 #include "Function2D.h"
 #include "Function3D.h"
 #include "Matrix.h"
+#include "config/config.h"
+
+#include "CudaEnvInit.h"
+
+#include <algorithm>
 #include <cmath>
 #include <iostream>
 using namespace Aurora;
@@ -103,4 +108,26 @@ namespace Recon {
         result.senderCoordList = aVSenderCoordList;
         return result;
     }
+    
+    ArtResult reconstructArt(Aurora::Matrix &data, Aurora::Matrix &dataAtt,
+                             Aurora::Matrix &dims, 
+                             Aurora::Matrix &senderList,
+                             Aurora::Matrix &receiverList, Aurora::Matrix &res,
+                             double SOS_IN_WATER)
+    {
+        auto nTotalRays = size(senderList, 2);
+        int numIter = 1;
+        if (transParams::bentReconstruction)
+        {
+            numIter =transParams::bentIter+1; 
+        }
+        for (size_t i = 0; i < transParams::gpuSelectionList.getDataSize(); i++)
+        {
+            std::string msg;
+            if (!resetGPUDevice((int)transParams::gpuSelectionList[i],msg))
+            {
+                std::cerr<<msg<<std::endl;
+            }
+        }
+    }
 } // namespace Recon

src/transmissionReconstruction/reconstruction/reconstruction.h

@@ -9,6 +9,13 @@ struct DiscretizePositionValues{
     Aurora::Matrix dims;
     Aurora::Matrix res;
 };
+struct ArtResult{
+    Aurora::Matrix outSOS;
+    Aurora::Matrix outATT;
+    Aurora::Matrix outAll;
+    bool bentReco;
+    Aurora::Matrix allHitMaps;
+};
 Aurora::Matrix calculateMinimalMaximalTransducerPositions(
    const Aurora::Matrix &aMSenderList, const Aurora::Matrix &aMReceiverList);
 
@@ -20,6 +27,12 @@ void slownessToSOS(Aurora::Matrix & aVF1, double aSOS_IN_WATER);
 
 DiscretizePositionValues discretizePositions(Aurora::Matrix &aVSenderCoordList, Aurora::Matrix &aVReceiverCoordList, double aNumPixelXY);
 
+ArtResult reconstructArt(Aurora::Matrix &data, Aurora::Matrix &dataAtt,
+                         Aurora::Matrix &dims, 
+                         Aurora::Matrix &senderList,
+                         Aurora::Matrix &receiverList, Aurora::Matrix &res,
+                         double SOS_IN_WATER);
+
 //TODO: exitBent, reconstructArt, saveTransmissionReconstruction.m
 
 }