#ifndef PostEddyAlignShellsFunctions_h
#define PostEddyAlignShellsFunctions_h

#include "EddyHelperClasses.h"
#include "ECScanClasses.h"
#include "DiffusionGP.h"
#include "b0Predictor.h"
#include "EddyUtils.h"
#include "EddyCommandLineOptions.h"

namespace EDDY {

/****************************************************************//**
*
* \brief This class does Post Eddy Alignment of Shells (PEAS)
*
* 
*
********************************************************************/ 

class PEASUtils
{
public:
  /// Attempts to align the shells after eddy has done all within shell alignments
  static void PostEddyAlignShells(// Input
				  const EddyCommandLineOptions&   clo,
				  bool                            upe, // Update parameter estimates or not
				  // Input/Output
				  ECScanManager&                  sm);
  /// Aligns shells through PE-direction translation _only_
  static void PostEddyAlignShellsAlongPE(// Input
					 const EddyCommandLineOptions&   clo,
					 bool                            upe, // Update parameter estimates or not
					 // Input/Output
					 ECScanManager&                  sm);
  /// Generate text file with MI values for different movement parameters. For debugging purposes only.
  static void WritePostEddyBetweenShellMIValues(// Input
						const EddyCommandLineOptions&     clo,
						const ECScanManager&              sm,
						const std::vector<unsigned int>&  n,
						const std::vector<double>&        first,
						const std::vector<double>&        last,
						const std::string&                bfname);
private:
  static NEWIMAGE::volume<float> get_mean_scan(// Input
					       const EddyCommandLineOptions&     clo,      
					       const ECScanManager&              sm,
					       const std::vector<unsigned int>&  indx,
					       // Output
					       NEWIMAGE::volume<float>&          mask);

  static void get_mean_scan_wrapper(unsigned int                      first_index,
				    unsigned int                      last_index,
				    const ECScanManager&              sm,
				    const std::vector<unsigned int>&  indx,
				    bool                              very_verbose,
				    // Output
				    NEWIMAGE::volume<float>&          mean,
				    NEWIMAGE::volume<float>&          mask);

  static NEWMAT::ColumnVector register_volumes(// Input
					       const NEWIMAGE::volume<float>& ref,
					       const NEWIMAGE::volume<float>& ima,
					       const NEWIMAGE::volume<float>& mask,
					       // Output
					       NEWIMAGE::volume<float>&       rima);

  static NEWMAT::ColumnVector register_volumes_along_PE(// Input
							const NEWIMAGE::volume<float>& ref,
							const NEWIMAGE::volume<float>& ima,
							const NEWIMAGE::volume<float>& mask,
							unsigned int                   pe_dir,
							// Output
							NEWIMAGE::volume<float>&       rima);

  static void register_volumes_wrapper(// Input
				       unsigned int                                     first_reg,
				       unsigned int                                     last_reg,
				       const std::vector<std::pair<int,unsigned int> >& regs,
				       bool                                             pe_only,
				       unsigned int                                     pe_dir,
				       const NEWIMAGE::volume<float>&                   b0,
				       const std::vector<NEWIMAGE::volume<float> >&     dwis,
				       const NEWIMAGE::volume<float>&                   mask,
				       const std::vector<double>&                       grpb,
				       bool                                             vv,
				       // Output
				       std::vector<NEWMAT::ColumnVector>&               mov_par,
				       std::vector<std::vector<NEWMAT::ColumnVector> >& cmov_par);

  static std::string format_mp(const NEWMAT::ColumnVector& mp);

  static std::vector<NEWMAT::ColumnVector> collate_mov_par_estimates_for_use(const std::vector<NEWMAT::ColumnVector>&                mp,
									     const std::vector<std::vector<NEWMAT::ColumnVector> >&  cmp,
									     const NEWIMAGE::volume<float>&                          ima);

  static void write_post_eddy_align_shells_report(const std::vector<NEWMAT::ColumnVector>&                mi_dmp,
						  const std::vector<NEWMAT::ColumnVector>&                mi_ump,
						  const std::vector<std::vector<NEWMAT::ColumnVector> >&  mi_cmp,
						  const std::vector<NEWMAT::ColumnVector>&                b0_dmp,
						  const std::vector<NEWMAT::ColumnVector>&                b0_ump,
						  const std::vector<double>&                              grpb,
						  bool                                                    upe,
						  const EddyCommandLineOptions&                           clo);

  static void write_post_eddy_align_shells_along_PE_report(const std::vector<NEWMAT::ColumnVector>&                mi_dmp,
							   const std::vector<NEWMAT::ColumnVector>&                mi_ump,
							   const std::vector<std::vector<NEWMAT::ColumnVector> >&  mi_cmp,
							   const std::vector<double>&                              grpb,
							   bool                                                    upe,
							   const EddyCommandLineOptions&                           clo);

  static void update_mov_par_estimates(// Input
				       const NEWMAT::ColumnVector&       mp,
				       const std::vector<unsigned int>&  indx,
				       // Input/output
				       ECScanManager&                    sm);

  static void align_shells_using_MI(// Input
				    const EddyCommandLineOptions&                      clo,
				    bool                                               pe_only,
				    // Input/Output
				    ECScanManager&                                     sm,
				    // Output
				    std::vector<double>&                               grpb,
				    std::vector<NEWMAT::ColumnVector>&                 mov_par,
				    std::vector<std::vector<NEWMAT::ColumnVector> >&   cmov_par,
				    std::vector<NEWMAT::ColumnVector>&                 mp_for_updates);

  static void write_between_shell_MI_values(const NEWIMAGE::volume<float>&    b0_mean,
					    const NEWIMAGE::volume<float>&    dwi_means,
					    const NEWIMAGE::volume<float>&    mask,
					    const std::string&                fname,
					    const std::vector<unsigned int>&  n,
					    const std::vector<double>&        first,
					    const std::vector<double>&        last);

  static void align_shells_using_interspersed_B0_scans(// Input
						       const EddyCommandLineOptions&                      clo,
						       // Input/Output
						       ECScanManager&                                     sm,
						       // Output
						       std::vector<double>&                               grpb,
						       std::vector<NEWMAT::ColumnVector>&                 mov_par,
						       std::vector<NEWMAT::ColumnVector>&                 mp_for_updates);
}; // End of class PEASUtils

} // End namespace EDDY

#endif // end #ifndef PostEddyAlignShellsFunctions_h