simple generator

src/generator-rf.h

@@ -12,6 +12,7 @@ namespace mri {
     void SetT2sDecaysNorm(std::vector<double> T2s_decays_norm);
     void SetTimer(long total_periods, double period_ms, int samples_per_period);
     void SetVoxels(double voxel_num, double phase_init_, double phase_range_, double T2s_scale);
+    const std::vector<double>& GetTimesteps(){return timesteps_};
   private:
     std::vector<double> amplitudes_;
     std::vector<double> phases_;
@@ -33,7 +34,6 @@ namespace mri {
     long total_samples_ = total_periods_ * samples_per_period_;
 
     std::vector<double> timesteps_;
-
-  };
+  };  // end of class GeneratorRF
 }  // end of namespace mri
 #endif  // SRC_GENERATOR_RF_H_

src/lock-in.h

@@ -1,4 +1,9 @@
 #ifndef SRC_LOCK_IN_H_
 #define SRC_LOCK_IN_H_
-
+namespace mri {
+  class LockIn {
+  public:
+    void Analyse();
+  };  // end of class LockIn
+}  // end of namespace mri
 #endif  // SRC_LOCK_IN_H_

src/main.cc

@@ -45,8 +45,10 @@ int main(int argc, char *argv[]) {
     generator_rf.SetVoxels(voxel_num, phase_init, phase_range, T2s_scale);
     generator_rf.SetTimer(total_periods, period_ms, rf_samples_per_period);
     generator_rf.RandomizeVoxels();
-    std::vector<double> signal_rf;
+    std::vector<double> signal_rf, magnitude_rf, phase_rf;
     generator_rf.EvaluateSignal(signal_rf);
+    auto timesteps = generator_rf.GetTimesteps();
+    std::cout<< timesteps.size() <<std::endl;
     // for (auto value:signal_rf){
     //   std::cout << value <<std::endl;
     // }