PulseChiSqSNNLS Class Reference

#include <PulseChiSqSNNLS.h>

Public Types
typedef BXVector::Index	Index

Public Member Functions
const BXVector &	BXs () const
double	ChiSq () const
void	disableErrorCalculation ()
bool	DoFit (const SampleVector &samples, const SampleMatrix &samplecov, const BXVector &bxs, const FullSampleVector &fullpulse, const FullSampleMatrix &fullpulsecov, const SampleGainVector &gains=-1 *SampleGainVector::Ones(), const SampleGainVector &badSamples=SampleGainVector::Zero())
const PulseVector &	Errors () const
const SampleMatrix &	invcov () const
	PulseChiSqSNNLS ()
const SamplePulseMatrix &	pulsemat () const
void	setMaxIters (int n)
void	setMaxIterWarnings (bool b)
const PulseVector &	X () const
	~PulseChiSqSNNLS ()

Protected Member Functions
double	ComputeApproxUncertainty (unsigned int ipulse)
double	ComputeChiSq ()
bool	Minimize (const SampleMatrix &samplecov, const FullSampleMatrix &fullpulsecov)
bool	NNLS ()
void	NNLSConstrainParameter (Index minratioidx)
void	NNLSUnconstrainParameter (Index idxp)
bool	OnePulseMinimize ()
bool	updateCov (const SampleMatrix &samplecov, const FullSampleMatrix &fullpulsecov)

Protected Attributes
PulseVector	_ampvec
PulseVector	_ampvecmin
BXVector	_bxs
BXVector	_bxsmin
double	_chisq
bool	_computeErrors
SampleDecompLLT	_covdecomp
SampleMatrix	_covdecompLinv
PulseVector	_errvec
SampleMatrix	_invcov
int	_maxiters
bool	_maxiterwarnings
unsigned int	_nP
unsigned int	_npulsetot
PulseDecompLDLT	_pulsedecomp
SamplePulseMatrix	_pulsemat
SampleVector	_sampvec
PulseMatrix	_topleft_work
PulseVector	ampvecpermtest
PulseMatrix	aTamat
PulseVector	aTbvec
SamplePulseMatrix	invcovp
PulseVector	updatework

Detailed Description

Definition at line 10 of file PulseChiSqSNNLS.h.

Member Typedef Documentation

Index

typedef BXVector::Index PulseChiSqSNNLS::Index

Definition at line 12 of file PulseChiSqSNNLS.h.

Constructor & Destructor Documentation

PulseChiSqSNNLS()

PulseChiSqSNNLS::PulseChiSqSNNLS

(

)

Definition at line 55 of file PulseChiSqSNNLS.cc.

: _chisq(0.), _computeErrors(true), _maxiters(50), _maxiterwarnings(true) {}

~PulseChiSqSNNLS()

PulseChiSqSNNLS::~PulseChiSqSNNLS

(

)

Definition at line 57 of file PulseChiSqSNNLS.cc.

{}

Member Function Documentation

BXs()

const BXVector& PulseChiSqSNNLS::BXs ( void  ) const

inline

Definition at line 30 of file PulseChiSqSNNLS.h.

References _bxsmin.

Referenced by EcalUncalibRecHitMultiFitAlgo::makeRecHit().

{ return _bxsmin; }

ChiSq()

double PulseChiSqSNNLS::ChiSq ( ) const

inline

Definition at line 32 of file PulseChiSqSNNLS.h.

References _chisq.

Referenced by EcalUncalibRecHitMultiFitAlgo::makeRecHit().

{ return _chisq; }

ComputeApproxUncertainty()

double PulseChiSqSNNLS::ComputeApproxUncertainty ( unsigned int  ipulse )

protected

Definition at line 289 of file PulseChiSqSNNLS.cc.

References _covdecomp, and _pulsemat.

Referenced by DoFit().

                                                                    {
  //compute approximate uncertainties
  //(using 1/second derivative since full Hessian is not meaningful in
  //presence of positive amplitude boundaries.)

  return 1. / _covdecomp.matrixL().solve(_pulsemat.col(ipulse)).norm();
}

ComputeChiSq()

double PulseChiSqSNNLS::ComputeChiSq ( )

protected

Definition at line 282 of file PulseChiSqSNNLS.cc.

References _ampvec, _covdecomp, _pulsemat, and _sampvec.

Referenced by DoFit(), and Minimize().

                                     {
  //   SampleVector resvec = _pulsemat*_ampvec - _sampvec;
  //   return resvec.transpose()*_covdecomp.solve(resvec);

  return _covdecomp.matrixL().solve(_pulsemat * _ampvec - _sampvec).squaredNorm();
}

disableErrorCalculation()

void PulseChiSqSNNLS::disableErrorCalculation ( )

inline

Definition at line 33 of file PulseChiSqSNNLS.h.

References _computeErrors.

Referenced by EcalUncalibRecHitMultiFitAlgo::EcalUncalibRecHitMultiFitAlgo(), and EcalUncalibRecHitMultiFitAlgo::makeRecHit().

{ _computeErrors = false; }

DoFit()

bool PulseChiSqSNNLS::DoFit	(	const SampleVector &	samples,
		const SampleMatrix &	samplecov,
		const BXVector &	bxs,
		const FullSampleVector &	fullpulse,
		const FullSampleMatrix &	fullpulsecov,
		const SampleGainVector &	gains = -1 * SampleGainVector::Ones(),
		const SampleGainVector &	badSamples = SampleGainVector::Zero()
	)

Definition at line 59 of file PulseChiSqSNNLS.cc.

References _ampvec, _ampvecmin, _bxs, _bxsmin, _chisq, _computeErrors, _errvec, _nP, _npulsetot, _pulsemat, _sampvec, funct::abs(), aTamat, nano_mu_digi_cff::bx, ComputeApproxUncertainty(), ComputeChiSq(), ecalph2::gains, mps_fire::i, gpuClustering::pixelStatus::mask, SiStripPI::max, Minimize(), NNLSUnconstrainParameter(), hltrates_dqm_sourceclient-live_cfg::offset, EcalCondDBWriter_cfi::pedestal, BXVector< T >::resize(), EgammaValidation_cff::samples, mathSSE::sqrt(), mps_update::status, and edm::swap().

Referenced by EcalUncalibRecHitMultiFitAlgo::makeRecHit().

                                                                {
  int npulse = bxs.rows();

  _sampvec = samples;
  _bxs = bxs;
  _pulsemat.resize(Eigen::NoChange, npulse);

  //construct dynamic pedestals if applicable
  int ngains = gains.maxCoeff() + 1;
  int nPedestals = 0;
  for (int gainidx = 0; gainidx < ngains; ++gainidx) {
    SampleGainVector mask = gainidx * SampleGainVector::Ones();
    SampleVector pedestal = (gains.array() == mask.array()).cast<SampleVector::value_type>();
    if (pedestal.maxCoeff() > 0.) {
      ++nPedestals;
      _bxs.resize(npulse + nPedestals);
      _bxs[npulse + nPedestals - 1] = 100 + gainidx;  //bx values >=100 indicate dynamic pedestals
      _pulsemat.resize(Eigen::NoChange, npulse + nPedestals);
      _pulsemat.col(npulse + nPedestals - 1) = pedestal;
    }
  }

  //construct negative step functions for saturated or potentially slew-rate-limited samples
  for (int isample = 0; isample < SampleVector::RowsAtCompileTime; ++isample) {
    if (badSamples.coeff(isample) > 0) {
      SampleVector step = SampleVector::Zero();
      //step correction has negative sign for saturated or slew-limited samples which have been forced to zero
      step[isample] = -1.;

      ++nPedestals;
      _bxs.resize(npulse + nPedestals);
      _bxs[npulse + nPedestals - 1] =
          -100 - isample;  //bx values <=-100 indicate step corrections for saturated or slew-limited samples
      _pulsemat.resize(Eigen::NoChange, npulse + nPedestals);
      _pulsemat.col(npulse + nPedestals - 1) = step;
    }
  }

  _npulsetot = npulse + nPedestals;

  _ampvec = PulseVector::Zero(_npulsetot);
  _errvec = PulseVector::Zero(_npulsetot);
  _nP = 0;
  _chisq = 0.;

  if (_bxs.rows() == 1 && std::abs(_bxs.coeff(0)) < 100) {
    _ampvec.coeffRef(0) = _sampvec.coeff(_bxs.coeff(0) + 5);
  }

  aTamat.resize(_npulsetot, _npulsetot);

  //initialize pulse template matrix
  for (int ipulse = 0; ipulse < npulse; ++ipulse) {
    int bx = _bxs.coeff(ipulse);
    int offset = 7 - 3 - bx;
    _pulsemat.col(ipulse) = fullpulse.segment<SampleVector::RowsAtCompileTime>(offset);
  }

  //unconstrain pedestals already for first iteration since they should always be non-zero
  if (nPedestals > 0) {
    for (int i = 0; i < _bxs.rows(); ++i) {
      int bx = _bxs.coeff(i);
      if (bx >= 100) {
        NNLSUnconstrainParameter(i);
      }
    }
  }

  //do the actual fit
  bool status = Minimize(samplecov, fullpulsecov);
  _ampvecmin = _ampvec;
  _bxsmin = _bxs;

  if (!status)
    return status;

  if (!_computeErrors)
    return status;

  //compute MINOS-like uncertainties for in-time amplitude
  bool foundintime = false;
  unsigned int ipulseintime = 0;
  for (unsigned int ipulse = 0; ipulse < _npulsetot; ++ipulse) {
    if (_bxs.coeff(ipulse) == 0) {
      ipulseintime = ipulse;
      foundintime = true;
      break;
    }
  }
  if (!foundintime)
    return status;

  const unsigned int ipulseintimemin = ipulseintime;

  double approxerr = ComputeApproxUncertainty(ipulseintime);
  double chisq0 = _chisq;
  double x0 = _ampvecmin[ipulseintime];

  //move in time pulse first to active set if necessary
  if (ipulseintime < _nP) {
    _pulsemat.col(_nP - 1).swap(_pulsemat.col(ipulseintime));
    std::swap(_ampvec.coeffRef(_nP - 1), _ampvec.coeffRef(ipulseintime));
    std::swap(_bxs.coeffRef(_nP - 1), _bxs.coeffRef(ipulseintime));
    ipulseintime = _nP - 1;
    --_nP;
  }

  SampleVector pulseintime = _pulsemat.col(ipulseintime);
  _pulsemat.col(ipulseintime).setZero();

  //two point interpolation for upper uncertainty when amplitude is away from boundary
  double xplus100 = x0 + approxerr;
  _ampvec.coeffRef(ipulseintime) = xplus100;
  _sampvec = samples - _ampvec.coeff(ipulseintime) * pulseintime;
  status &= Minimize(samplecov, fullpulsecov);
  if (!status)
    return status;
  double chisqplus100 = ComputeChiSq();

  double sigmaplus = std::abs(xplus100 - x0) / sqrt(chisqplus100 - chisq0);

  //if amplitude is sufficiently far from the boundary, compute also the lower uncertainty and average them
  if ((x0 / sigmaplus) > 0.5) {
    for (unsigned int ipulse = 0; ipulse < _npulsetot; ++ipulse) {
      if (_bxs.coeff(ipulse) == 0) {
        ipulseintime = ipulse;
        break;
      }
    }
    double xminus100 = std::max(0., x0 - approxerr);
    _ampvec.coeffRef(ipulseintime) = xminus100;
    _sampvec = samples - _ampvec.coeff(ipulseintime) * pulseintime;
    status &= Minimize(samplecov, fullpulsecov);
    if (!status)
      return status;
    double chisqminus100 = ComputeChiSq();

    double sigmaminus = std::abs(xminus100 - x0) / sqrt(chisqminus100 - chisq0);
    _errvec[ipulseintimemin] = 0.5 * (sigmaplus + sigmaminus);

  } else {
    _errvec[ipulseintimemin] = sigmaplus;
  }

  _chisq = chisq0;

  return status;
}

Errors()

const PulseVector& PulseChiSqSNNLS::Errors ( ) const

inline

Definition at line 29 of file PulseChiSqSNNLS.h.

References _errvec.

Referenced by EcalUncalibRecHitMultiFitAlgo::makeRecHit().

{ return _errvec; }

invcov()

const SampleMatrix& PulseChiSqSNNLS::invcov ( ) const

inline

Definition at line 26 of file PulseChiSqSNNLS.h.

References _invcov.

{ return _invcov; }

Minimize()

bool PulseChiSqSNNLS::Minimize ( const SampleMatrix &  samplecov, const FullSampleMatrix &  fullpulsecov  )

protected

Definition at line 214 of file PulseChiSqSNNLS.cc.

References _bxs, _chisq, _maxiters, _maxiterwarnings, funct::abs(), ComputeChiSq(), MillePedeFileConverter_cfg::e, LogDebug, NNLS(), OnePulseMinimize(), mps_update::status, and updateCov().

Referenced by DoFit().

                                                                                                  {
  const unsigned int npulse = _bxs.rows();

  int iter = 0;
  bool status = false;
  while (true) {
    if (iter >= _maxiters) {
      if (_maxiterwarnings) {
        LogDebug("PulseChiSqSNNLS::Minimize") << "Max Iterations reached at iter " << iter;
      }
      break;
    }

    status = updateCov(samplecov, fullpulsecov);
    if (!status)
      break;
    if (npulse > 1) {
      status = NNLS();
    } else {
      //special case for one pulse fit (performance optimized)
      status = OnePulseMinimize();
    }
    if (!status)
      break;

    double chisqnow = ComputeChiSq();
    double deltachisq = chisqnow - _chisq;

    _chisq = chisqnow;
    if (std::abs(deltachisq) < 1e-3) {
      break;
    }
    ++iter;
  }

  return status;
}

NNLS()

bool PulseChiSqSNNLS::NNLS ( )

protected

Definition at line 297 of file PulseChiSqSNNLS.cc.

References _ampvec, _bxs, _covdecomp, _nP, _pulsemat, _sampvec, ampvecpermtest, aTamat, aTbvec, MillePedeFileConverter_cfg::e, eigen_solve_submatrix(), mps_fire::i, invcovp, LogDebug, SiStripPI::max, SiStripPI::min, NNLSConstrainParameter(), NNLSUnconstrainParameter(), particleFlowDisplacedVertex_cfi::ratio, remoteMonitoring_LASER_era2018_cfg::threshold, updatework, and cms::cuda::wmax.

Referenced by Minimize().

                           {
  //Fast NNLS (fnnls) algorithm as per http://citeseerx.ist.psu.edu/viewdoc/download?doi=10.1.1.157.9203&rep=rep1&type=pdf

  const unsigned int npulse = _bxs.rows();
  constexpr unsigned int nsamples = SampleVector::RowsAtCompileTime;

  invcovp = _covdecomp.matrixL().solve(_pulsemat);
  aTamat.noalias() = invcovp.transpose().lazyProduct(invcovp);
  aTbvec.noalias() = invcovp.transpose().lazyProduct(_covdecomp.matrixL().solve(_sampvec));

  int iter = 0;
  Index idxwmax = 0;
  double wmax = 0.0;
  double threshold = 1e-11;
  while (true) {
    //can only perform this step if solution is guaranteed viable
    if (iter > 0 || _nP == 0) {
      if (_nP == std::min(npulse, nsamples))
        break;

      const unsigned int nActive = npulse - _nP;

      updatework = aTbvec - aTamat * _ampvec;
      Index idxwmaxprev = idxwmax;
      double wmaxprev = wmax;
      wmax = updatework.tail(nActive).maxCoeff(&idxwmax);

      //convergence
      if (wmax < threshold || (idxwmax == idxwmaxprev && wmax == wmaxprev))
        break;

      //worst case protection
      if (iter >= 500) {
        LogDebug("PulseChiSqSNNLS::NNLS()") << "Max Iterations reached at iter " << iter;
        break;
      }

      //unconstrain parameter
      Index idxp = _nP + idxwmax;
      NNLSUnconstrainParameter(idxp);
    }

    while (true) {
      //printf("iter in, idxsP = %i\n",int(_idxsP.size()));

      if (_nP == 0)
        break;

      ampvecpermtest = _ampvec;

      //solve for unconstrained parameters
      //need to have specialized function to call optimized versions
      // of matrix solver... this is truly amazing...
      eigen_solve_submatrix(aTamat, aTbvec, ampvecpermtest, _nP);

      //check solution
      bool positive = true;
      for (unsigned int i = 0; i < _nP; ++i)
        positive &= (ampvecpermtest(i) > 0);
      if (positive) {
        _ampvec.head(_nP) = ampvecpermtest.head(_nP);
        break;
      }

      //update parameter vector
      Index minratioidx = 0;

      // no realizable optimization here (because it autovectorizes!)
      double minratio = std::numeric_limits<double>::max();
      for (unsigned int ipulse = 0; ipulse < _nP; ++ipulse) {
        if (ampvecpermtest.coeff(ipulse) <= 0.) {
          const double c_ampvec = _ampvec.coeff(ipulse);
          const double ratio = c_ampvec / (c_ampvec - ampvecpermtest.coeff(ipulse));
          if (ratio < minratio) {
            minratio = ratio;
            minratioidx = ipulse;
          }
        }
      }

      _ampvec.head(_nP) += minratio * (ampvecpermtest.head(_nP) - _ampvec.head(_nP));

      //avoid numerical problems with later ==0. check
      _ampvec.coeffRef(minratioidx) = 0.;

      //printf("removing index %i, orig idx %i\n",int(minratioidx),int(_bxs.coeff(minratioidx)));
      NNLSConstrainParameter(minratioidx);
    }
    ++iter;

    //adaptive convergence threshold to avoid infinite loops but still
    //ensure best value is used
    if (iter % 16 == 0) {
      threshold *= 2;
    }
  }

  return true;
}

NNLSConstrainParameter()

void PulseChiSqSNNLS::NNLSConstrainParameter ( Index  minratioidx )

protected

Definition at line 407 of file PulseChiSqSNNLS.cc.

References _ampvec, _bxs, _nP, _pulsemat, aTamat, aTbvec, and edm::swap().

Referenced by NNLS().

                                                              {
  aTamat.col(_nP - 1).swap(aTamat.col(minratioidx));
  aTamat.row(_nP - 1).swap(aTamat.row(minratioidx));
  _pulsemat.col(_nP - 1).swap(_pulsemat.col(minratioidx));
  std::swap(aTbvec.coeffRef(_nP - 1), aTbvec.coeffRef(minratioidx));
  std::swap(_ampvec.coeffRef(_nP - 1), _ampvec.coeffRef(minratioidx));
  std::swap(_bxs.coeffRef(_nP - 1), _bxs.coeffRef(minratioidx));
  --_nP;
}

NNLSUnconstrainParameter()

void PulseChiSqSNNLS::NNLSUnconstrainParameter ( Index  idxp )

protected

Definition at line 397 of file PulseChiSqSNNLS.cc.

References _ampvec, _bxs, _nP, _pulsemat, aTamat, aTbvec, and edm::swap().

Referenced by DoFit(), and NNLS().

                                                         {
  aTamat.col(_nP).swap(aTamat.col(idxp));
  aTamat.row(_nP).swap(aTamat.row(idxp));
  _pulsemat.col(_nP).swap(_pulsemat.col(idxp));
  std::swap(aTbvec.coeffRef(_nP), aTbvec.coeffRef(idxp));
  std::swap(_ampvec.coeffRef(_nP), _ampvec.coeffRef(idxp));
  std::swap(_bxs.coeffRef(_nP), _bxs.coeffRef(idxp));
  ++_nP;
}

OnePulseMinimize()

bool PulseChiSqSNNLS::OnePulseMinimize ( )

protected

Definition at line 417 of file PulseChiSqSNNLS.cc.

References _ampvec, _covdecomp, _pulsemat, _sampvec, invcovp, and SiStripPI::max.

Referenced by Minimize().

                                       {
  //Fast NNLS (fnnls) algorithm as per http://citeseerx.ist.psu.edu/viewdoc/download?doi=10.1.1.157.9203&rep=rep1&type=pdf

  //   const unsigned int npulse = 1;

  invcovp = _covdecomp.matrixL().solve(_pulsemat);
  //   aTamat = invcovp.transpose()*invcovp;
  //   aTbvec = invcovp.transpose()*_covdecomp.matrixL().solve(_sampvec);

  SingleMatrix aTamatval = invcovp.transpose() * invcovp;
  SingleVector aTbvecval = invcovp.transpose() * _covdecomp.matrixL().solve(_sampvec);
  _ampvec.coeffRef(0) = std::max(0., aTbvecval.coeff(0) / aTamatval.coeff(0));

  return true;
}

pulsemat()

const SamplePulseMatrix& PulseChiSqSNNLS::pulsemat ( ) const

inline

Definition at line 25 of file PulseChiSqSNNLS.h.

References _pulsemat.

{ return _pulsemat; }

setMaxIters()

void PulseChiSqSNNLS::setMaxIters ( int  n )

inline

Definition at line 34 of file PulseChiSqSNNLS.h.

References _maxiters, and dqmiodumpmetadata::n.

Referenced by EcalUncalibRecHitMultiFitAlgo::EcalUncalibRecHitMultiFitAlgo().

{ _maxiters = n; }

setMaxIterWarnings()

void PulseChiSqSNNLS::setMaxIterWarnings ( bool  b )

inline

Definition at line 35 of file PulseChiSqSNNLS.h.

References _maxiterwarnings, and b.

Referenced by EcalUncalibRecHitMultiFitAlgo::EcalUncalibRecHitMultiFitAlgo().

{ _maxiterwarnings = b; }

updateCov()

bool PulseChiSqSNNLS::updateCov ( const SampleMatrix &  samplecov, const FullSampleMatrix &  fullpulsecov  )

protected

Definition at line 252 of file PulseChiSqSNNLS.cc.

References _ampvec, _bxs, _covdecomp, _invcov, funct::abs(), nano_mu_digi_cff::bx, SiStripPI::max, hltrates_dqm_sourceclient-live_cfg::offset, and mps_update::status.

Referenced by Minimize().

                                                                                                   {
  const unsigned int nsample = SampleVector::RowsAtCompileTime;
  const unsigned int npulse = _bxs.rows();

  _invcov = samplecov;  //

  for (unsigned int ipulse = 0; ipulse < npulse; ++ipulse) {
    if (_ampvec.coeff(ipulse) == 0.)
      continue;
    int bx = _bxs.coeff(ipulse);
    if (std::abs(bx) >= 100)
      continue;  //no contribution to covariance from pedestal or saturation/slew step correction

    int firstsamplet = std::max(0, bx + 3);
    int offset = 7 - 3 - bx;

    const double ampveccoef = _ampvec.coeff(ipulse);
    const double ampsq = ampveccoef * ampveccoef;

    const unsigned int nsamplepulse = nsample - firstsamplet;
    _invcov.block(firstsamplet, firstsamplet, nsamplepulse, nsamplepulse) +=
        ampsq * fullpulsecov.block(firstsamplet + offset, firstsamplet + offset, nsamplepulse, nsamplepulse);
  }

  _covdecomp.compute(_invcov);

  bool status = true;
  return status;
}

X()

const PulseVector& PulseChiSqSNNLS::X ( ) const

inline

Definition at line 28 of file PulseChiSqSNNLS.h.

References _ampvecmin.

Referenced by svgfig.Curve.Sample::__repr__(), and EcalUncalibRecHitMultiFitAlgo::makeRecHit().

{ return _ampvecmin; }

Member Data Documentation

_ampvec

PulseVector PulseChiSqSNNLS::_ampvec

protected

Definition at line 50 of file PulseChiSqSNNLS.h.

Referenced by ComputeChiSq(), DoFit(), NNLS(), NNLSConstrainParameter(), NNLSUnconstrainParameter(), OnePulseMinimize(), and updateCov().

_ampvecmin

PulseVector PulseChiSqSNNLS::_ampvecmin

protected

Definition at line 52 of file PulseChiSqSNNLS.h.

Referenced by DoFit(), and X().

_bxs

BXVector PulseChiSqSNNLS::_bxs

protected

Definition at line 59 of file PulseChiSqSNNLS.h.

Referenced by DoFit(), Minimize(), NNLS(), NNLSConstrainParameter(), NNLSUnconstrainParameter(), and updateCov().

_bxsmin

BXVector PulseChiSqSNNLS::_bxsmin

protected

Definition at line 60 of file PulseChiSqSNNLS.h.

Referenced by BXs(), and DoFit().

_chisq

double PulseChiSqSNNLS::_chisq

protected

Definition at line 71 of file PulseChiSqSNNLS.h.

Referenced by ChiSq(), DoFit(), and Minimize().

_computeErrors

bool PulseChiSqSNNLS::_computeErrors

protected

Definition at line 72 of file PulseChiSqSNNLS.h.

Referenced by disableErrorCalculation(), and DoFit().

_covdecomp

SampleDecompLLT PulseChiSqSNNLS::_covdecomp

protected

Definition at line 54 of file PulseChiSqSNNLS.h.

Referenced by ComputeApproxUncertainty(), ComputeChiSq(), NNLS(), OnePulseMinimize(), and updateCov().

_covdecompLinv

SampleMatrix PulseChiSqSNNLS::_covdecompLinv

protected

Definition at line 55 of file PulseChiSqSNNLS.h.

_errvec

PulseVector PulseChiSqSNNLS::_errvec

protected

Definition at line 51 of file PulseChiSqSNNLS.h.

Referenced by DoFit(), and Errors().

_invcov

SampleMatrix PulseChiSqSNNLS::_invcov

protected

Definition at line 48 of file PulseChiSqSNNLS.h.

Referenced by invcov(), and updateCov().

_maxiters

int PulseChiSqSNNLS::_maxiters

protected

Definition at line 73 of file PulseChiSqSNNLS.h.

Referenced by Minimize(), and setMaxIters().

_maxiterwarnings

bool PulseChiSqSNNLS::_maxiterwarnings

protected

Definition at line 74 of file PulseChiSqSNNLS.h.

Referenced by Minimize(), and setMaxIterWarnings().

_nP

unsigned int PulseChiSqSNNLS::_nP

protected

Definition at line 62 of file PulseChiSqSNNLS.h.

Referenced by DoFit(), NNLS(), NNLSConstrainParameter(), and NNLSUnconstrainParameter().

_npulsetot

unsigned int PulseChiSqSNNLS::_npulsetot

protected

Definition at line 61 of file PulseChiSqSNNLS.h.

Referenced by DoFit().

_pulsedecomp

PulseDecompLDLT PulseChiSqSNNLS::_pulsedecomp

protected

Definition at line 57 of file PulseChiSqSNNLS.h.

_pulsemat

SamplePulseMatrix PulseChiSqSNNLS::_pulsemat

protected

Definition at line 49 of file PulseChiSqSNNLS.h.

Referenced by ComputeApproxUncertainty(), ComputeChiSq(), DoFit(), NNLS(), NNLSConstrainParameter(), NNLSUnconstrainParameter(), OnePulseMinimize(), and pulsemat().

_sampvec

SampleVector PulseChiSqSNNLS::_sampvec

protected

Definition at line 47 of file PulseChiSqSNNLS.h.

Referenced by ComputeChiSq(), DoFit(), NNLS(), and OnePulseMinimize().

_topleft_work

PulseMatrix PulseChiSqSNNLS::_topleft_work

protected

Definition at line 56 of file PulseChiSqSNNLS.h.

ampvecpermtest

PulseVector PulseChiSqSNNLS::ampvecpermtest

protected

Definition at line 69 of file PulseChiSqSNNLS.h.

Referenced by NNLS().

aTamat

PulseMatrix PulseChiSqSNNLS::aTamat

protected

Definition at line 65 of file PulseChiSqSNNLS.h.

Referenced by DoFit(), NNLS(), NNLSConstrainParameter(), and NNLSUnconstrainParameter().

aTbvec

PulseVector PulseChiSqSNNLS::aTbvec

protected

Definition at line 66 of file PulseChiSqSNNLS.h.

Referenced by NNLS(), NNLSConstrainParameter(), and NNLSUnconstrainParameter().

invcovp

SamplePulseMatrix PulseChiSqSNNLS::invcovp

protected

Definition at line 64 of file PulseChiSqSNNLS.h.

Referenced by NNLS(), and OnePulseMinimize().

updatework

PulseVector PulseChiSqSNNLS::updatework

protected

Definition at line 67 of file PulseChiSqSNNLS.h.

Referenced by NNLS().