Skip to content
Snippets Groups Projects
Commit a5af5dc6 authored by Tim Behrens's avatar Tim Behrens
Browse files

*** empty log message ***

parent f5d0ef77
No related branches found
No related tags found
No related merge requests found
...@@ -58,13 +58,8 @@ namespace FIBRE{ ...@@ -58,13 +58,8 @@ namespace FIBRE{
const ColumnVector& m_alpha; const ColumnVector& m_alpha;
const ColumnVector& m_beta; const ColumnVector& m_beta;
const Matrix& m_bvals; const Matrix& m_bvals;
<<<<<<< fibre.h
xfibresOptions& opts;
public:
=======
xfibresOptions& opts; xfibresOptions& opts;
public: public:
>>>>>>> 1.10
//constructors:: //constructors::
Fibre( const ColumnVector& alpha, const ColumnVector& beta, Fibre( const ColumnVector& alpha, const ColumnVector& beta,
const Matrix& bvals,const float& d): const Matrix& bvals,const float& d):
...@@ -228,7 +223,7 @@ namespace FIBRE{ ...@@ -228,7 +223,7 @@ namespace FIBRE{
else{ else{
//m_f_prior=-(log(m_lam) + (m_lam-1)*log(1-m_f)); //m_f_prior=-(log(m_lam) + (m_lam-1)*log(1-m_f));
if(m_lam_jump) if(m_lam_jump)
m_f_prior=log(1-f)+2*log(fabs(log(1-f))); //marginalised with uniform prior on lambda m_f_prior=log(1-m_f)+2*log(fabs(log(1-m_f))); //marginalised with uniform prior on lambda
else else
m_f_prior=0; m_f_prior=0;
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment