inv_beamformer_compute.h

Go to the documentation of this file.

//=============================================================================================================
 
#ifndef INV_BEAMFORMER_COMPUTE_H
#define INV_BEAMFORMER_COMPUTE_H
 
//=============================================================================================================
// INCLUDES
//=============================================================================================================
 
#include "../inv_global.h"
#include "inv_beamformer_settings.h"
 
//=============================================================================================================
// EIGEN INCLUDES
//=============================================================================================================
 
#include <Eigen/Core>
 
//=============================================================================================================
// DEFINE NAMESPACE INVLIB
//=============================================================================================================
 
namespace INVLIB
{
 
//=============================================================================================================
class INVSHARED_EXPORT InvBeamformerCompute
{
public:
    //=========================================================================================================
    static bool computeBeamformer(const Eigen::MatrixXd& G,
                                  const Eigen::MatrixXd& Cm,
                                  double reg,
                                  int nOrient,
                                  BeamformerWeightNorm weightNorm,
                                  BeamformerPickOri pickOri,
                                  bool reduceRank,
                                  BeamformerInversion invMethod,
                                  const Eigen::MatrixX3d& nn,
                                  Eigen::MatrixXd& W,
                                  Eigen::MatrixX3d& maxPowerOri);
 
    //=========================================================================================================
    static Eigen::VectorXd computePower(const Eigen::MatrixXd& Cm,
                                        const Eigen::MatrixXd& W,
                                        int nOrient);
 
    //=========================================================================================================
    static Eigen::MatrixXd symMatPow(const Eigen::MatrixXd& X, double p, bool reduceRank = false);
 
private:
    //=========================================================================================================
    static void regPinv(const Eigen::MatrixXd& C,
                        double reg,
                        Eigen::MatrixXd& CInv,
                        double& loadingFactor,
                        int& rankOut);
 
    //=========================================================================================================
    static void reduceLeadfieldRank(Eigen::MatrixXd& Gk);
};
 
} // NAMESPACE INVLIB
 
#endif // INV_BEAMFORMER_COMPUTE_H