ProgImageRotationalPCA Class Reference

#include <image_rotational_pca.h>

Inheritance diagram for ProgImageRotationalPCA:

Collaboration diagram for ProgImageRotationalPCA:

Public Member Functions
	ProgImageRotationalPCA ()
	Empty constructor. More...
	~ProgImageRotationalPCA ()
	Destructor. More...
void	readParams ()
	Read argument from command line. More...
void	show ()
	Show. More...
void	defineParams ()
	Usage. More...
void	produceSideInfo ()
	Produce side info. More...
void	writeToHBuffer (int idx, double *dest)
void	flushHBuffer ()
void	clearHbuffer ()
void	applyT ()
void	applyTt ()
int	QR ()
void	run ()
virtual void	selectPartFromMd (MetaData &MDin)
virtual void	comunicateMatrix (Matrix2D< double > &W)
virtual void	createMutexes (size_t Nimgs)
virtual void	allReduceApplyT (Matrix2D< double > &Wnode_0)
virtual void	comunicateQrDim (int &qrDim)
virtual void	mapMatrix (int qrDim)
virtual void	applySVD ()
virtual void	copyHtoF (int block)
Public Member Functions inherited from XmippProgram
const char *	getParam (const char *param, int arg=0)
const char *	getParam (const char *param, const char *subparam, int arg=0)
int	getIntParam (const char *param, int arg=0)
int	getIntParam (const char *param, const char *subparam, int arg=0)
double	getDoubleParam (const char *param, int arg=0)
double	getDoubleParam (const char *param, const char *subparam, int arg=0)
float	getFloatParam (const char *param, int arg=0)
float	getFloatParam (const char *param, const char *subparam, int arg=0)
void	getListParam (const char *param, StringVector &list)
int	getCountParam (const char *param)
bool	checkParam (const char *param)
bool	existsParam (const char *param)
void	addParamsLine (const String &line)
void	addParamsLine (const char *line)
ParamDef *	getParamDef (const char *param) const
virtual void	quit (int exit_code=0) const
virtual int	tryRun ()
void	initProgress (size_t total, size_t stepBin=60)
void	setProgress (size_t value=0)
void	endProgress ()
void	processDefaultComment (const char *param, const char *left)
void	setDefaultComment (const char *param, const char *comment)
virtual void	initComments ()
void	setProgramName (const char *name)
void	addUsageLine (const char *line, bool verbatim=false)
void	clearUsage ()
void	addExampleLine (const char *example, bool verbatim=true)
void	addSeeAlsoLine (const char *seeAlso)
void	addKeywords (const char *keywords)
const char *	name () const
virtual void	usage (int verb=0) const
virtual void	usage (const String &param, int verb=2)
int	version () const
virtual void	show () const
virtual void	read (int argc, const char **argv, bool reportErrors=true)
virtual void	read (int argc, char **argv, bool reportErrors=true)
void	read (const String &argumentsLine)
	XmippProgram ()
	XmippProgram (int argc, const char **argv)
virtual	~XmippProgram ()

Public Attributes
FileName	fnIn
FileName	fnRoot
int	Neigen
int	Nits
double	psi_step
double	max_shift_change
double	shift_step
int	maxNimgs
int	Nthreads
int	rank
std::vector< MetaDataVec >	MD
size_t	Nimg
int	Nangles
int	Nshifts
size_t	Xdim
int	Npixels
std::vector< double * >	Hbuffer
std::vector< double * >	HbufferDestination
std::unique_ptr< Mutex >	fileMutex
std::unique_ptr< Mutex >	threadMutex
Matrix2D< double >	H
Matrix2D< double >	F
std::vector< Image< double > >	I
std::vector< MultidimArray< double > >	Iaux
std::vector< Matrix2D< double > >	A
std::vector< Matrix2D< double > >	Hblock
std::vector< Matrix2D< double > >	Wnode
Matrix2D< double >	Wtranspose
MultidimArray< unsigned char >	mask
std::unique_ptr< ThreadTaskDistributor >	taskDistributor
std::unique_ptr< ThreadManager >	thMgr
std::vector< size_t >	objId
Public Attributes inherited from XmippProgram
bool	doRun
bool	runWithoutArgs
int	verbose
	Verbosity level. More...
int	debug

Static Public Attributes
static const int	HbufferMax =20

Additional Inherited Members
Protected Member Functions inherited from XmippProgram
void	defineCommons ()
Protected Attributes inherited from XmippProgram
int	errorCode
ProgramDef *	progDef
	Program definition and arguments parser. More...
std::map< String, CommentList >	defaultComments
int	argc
	Original command line arguments. More...
const char **	argv

Detailed Description

Rotational invariant PCA parameters.

Definition at line 41 of file image_rotational_pca.h.

Constructor & Destructor Documentation

ProgImageRotationalPCA()

ProgImageRotationalPCA::ProgImageRotationalPCA

(

)

Empty constructor.

Definition at line 34 of file image_rotational_pca.cpp.

{
  rank = 0;
  verbose = 1;
}

~ProgImageRotationalPCA()

ProgImageRotationalPCA::~ProgImageRotationalPCA

(

)

Destructor.

Definition at line 41 of file image_rotational_pca.cpp.

{
  clearHbuffer();
}

Member Function Documentation

allReduceApplyT()

void ProgImageRotationalPCA::allReduceApplyT ( Matrix2D< double > &  Wnode_0 )

virtual

Last part of function applyT

Reimplemented in MpiProgImageRotationalPCA.

Definition at line 322 of file image_rotational_pca.cpp.

{
}

applySVD()

void ProgImageRotationalPCA::applySVD ( )

virtual

Apply SVD

Reimplemented in MpiProgImageRotationalPCA.

Definition at line 531 of file image_rotational_pca.cpp.

{
  // Apply SVD and extract the basis
    std::cout << "Performing SVD decomposition ..." << std::endl;
    // SVD of W
    Matrix2D<double> U,V;
    Matrix1D<double> S;
    svdcmp(Wnode[0],U,S,V);

    // Keep the first Neigen images from U
    Image<double> I;
    I().resizeNoCopy(Xdim,Xdim);
    const MultidimArray<double> &mI=I();
    FileName fnImg;
    MetaDataVec MD;
    for (int eig=0; eig<Neigen; eig++)
    {
      int Un=0;
      FOR_ALL_DIRECT_ELEMENTS_IN_MULTIDIMARRAY(mI)
      if (DIRECT_MULTIDIM_ELEM(mask,n))
      DIRECT_MULTIDIM_ELEM(mI,n)=MAT_ELEM(U,Un++,eig);
      fnImg.compose(eig+1,fnRoot,"stk");
      I.write(fnImg);
      size_t id=MD.addObject();
      MD.setValue(MDL_IMAGE,fnImg,id);
      MD.setValue(MDL_WEIGHT,VEC_ELEM(S,eig),id);
    }
    MD.write(fnRoot+".xmd");
}

applyT()

void ProgImageRotationalPCA::applyT

(

)

Apply T. W=T(H).

Definition at line 326 of file image_rotational_pca.cpp.

{
  Matrix2D<double> &Wnode_0=Wnode[0];
  Wnode_0.initZeros(Npixels,MAT_XSIZE(H));
  taskDistributor->reset();
  thMgr->run(threadApplyT);

  // Gather all Wnodes from all threads
  for (int n=1; n<Nthreads; ++n)
    Wnode_0 += Wnode[n];
  allReduceApplyT(Wnode_0);
  if (IS_MASTER)
    progress_bar(objId.size());
}

applyTt()

void ProgImageRotationalPCA::applyTt

(

)

Apply Tt. H=Tt(W).

Definition at line 445 of file image_rotational_pca.cpp.

{
  // Compute W transpose to accelerate memory access
  Matrix2D<double> &W=Wnode[0];
  Wtranspose.resizeNoCopy(MAT_XSIZE(W),MAT_YSIZE(W));
  FOR_ALL_ELEMENTS_IN_MATRIX2D(Wtranspose)
  MAT_ELEM(Wtranspose,i,j) = MAT_ELEM(W,j,i);

  taskDistributor->reset();
  thMgr->run(threadApplyTt);
  flushHBuffer();
  if (IS_MASTER)
    progress_bar(objId.size());
}

clearHbuffer()

void ProgImageRotationalPCA::clearHbuffer

(

)

Clear H buffer

Definition at line 47 of file image_rotational_pca.cpp.

{
  int nmax = Hbuffer.size();
  for (int n = 0; n < nmax; n++)
    delete[] Hbuffer[n];
  Hbuffer.clear();
}

comunicateMatrix()

void ProgImageRotationalPCA::comunicateMatrix ( Matrix2D< double > &  W )

virtual

Comunicate matrix, only meanful for MPI

Reimplemented in MpiProgImageRotationalPCA.

Definition at line 114 of file image_rotational_pca.cpp.

{
}

comunicateQrDim()

void ProgImageRotationalPCA::comunicateQrDim ( int &  qrDim )

virtual

Comunicate int param

Reimplemented in MpiProgImageRotationalPCA.

Definition at line 517 of file image_rotational_pca.cpp.

{
    std::cout << "Performing QR decomposition ..." << std::endl;
    qrDim = QR();
}

copyHtoF()

void ProgImageRotationalPCA::copyHtoF ( int  block )

virtual

Copy H to F.

Reimplemented in MpiProgImageRotationalPCA.

Definition at line 510 of file image_rotational_pca.cpp.

{
    size_t Hidx=block*MAT_XSIZE(H);
    FOR_ALL_ELEMENTS_IN_MATRIX2D(H)
    MAT_ELEM(F,Hidx+j,i)=MAT_ELEM(H,i,j);
}

createMutexes()

void ProgImageRotationalPCA::createMutexes ( size_t  Nimgs )

virtual

Create mutexes and distributor

Reimplemented in MpiProgImageRotationalPCA.

Definition at line 118 of file image_rotational_pca.cpp.

{
  fileMutex = std::make_unique<Mutex>();
  threadMutex = std::make_unique<Mutex>();
  taskDistributor = std::make_unique<ThreadTaskDistributor>(Nimgs, XMIPP_MAX(1,Nimgs/5));
}

defineParams()

void ProgImageRotationalPCA::defineParams ( )

virtual

Usage.

Reimplemented from XmippProgram.

Definition at line 88 of file image_rotational_pca.cpp.

{
  addUsageLine(
      "Makes a rotational invariant representation of the image collection");
  addParamsLine(
      "    -i <selfile>               : Selfile with experimental images");
  addParamsLine("   --oroot <rootname>          : Rootname for output");
  addParamsLine("  [--eigenvectors <N=200>]     : Number of eigenvectors");
  addParamsLine("  [--iterations <N=2>]         : Number of iterations");
  addParamsLine(
      "  [--max_shift_change <r=0>]   : Maximum change allowed in shift");
  addParamsLine("  [--psi_step <ang=1>]         : Step in psi in degrees");
  addParamsLine("  [--shift_step <r=1>]         : Step in shift in pixels");
  addParamsLine("  [--maxImages <N=-1>]         : Maximum number of images");
  addParamsLine("  [--thr <N=1>]                : Number of threads");
  addExampleLine("Typical use (4 nodes with 4 processors):", false);
  addExampleLine(
      "mpirun -np 4 `which xmipp_mpi_image_rotational_pca` -i images.stk --oroot images_eigen --thr 4");
}

flushHBuffer()

void ProgImageRotationalPCA::flushHBuffer

(

)

Flush buffer

Definition at line 215 of file image_rotational_pca.cpp.

{
  int nmax=HbufferDestination.size();
  const Matrix2D<double> &Hblock_0=Hblock[0];
  fileMutex->lock();
  for (int n=0; n<nmax; ++n)
  memcpy(HbufferDestination[n],Hbuffer[n],MAT_XSIZE(Hblock_0)*MAT_YSIZE(Hblock_0)*sizeof(double));
  fileMutex->unlock();
  HbufferDestination.clear();
}

mapMatrix()

void ProgImageRotationalPCA::mapMatrix ( int  qrDim )

virtual

Map matrix

Reimplemented in MpiProgImageRotationalPCA.

Definition at line 523 of file image_rotational_pca.cpp.

{
    FileName(fnRoot+"_matrixH.raw").createEmptyFileWithGivenLength(Nimg*2*Nangles*Nshifts*qrDim*sizeof(double));
    H.mapToFile(fnRoot+"_matrixH.raw",MAT_XSIZE(F),qrDim);
    FOR_ALL_ELEMENTS_IN_MATRIX2D(H)
      MAT_ELEM(H,i,j)=MAT_ELEM(F,j,i);
}

produceSideInfo()

void ProgImageRotationalPCA::produceSideInfo

(

)

Produce side info.

Definition at line 126 of file image_rotational_pca.cpp.

{
  time_config();
  MetaDataVec MDin(fnIn);

  if (maxNimgs > 0)
    selectPartFromMd(MDin);

    Nimg = MDin.size();
    size_t Ydim, Zdim, Ndim;
    getImageSize(MDin, Xdim, Ydim, Zdim, Ndim);
    Nangles = (int)floor(360.0 / psi_step);
    Nshifts = (int)((2 * max_shift_change + 1) / shift_step);
    Nshifts *= Nshifts;

    // Construct mask
    mask.resizeNoCopy(Xdim, Xdim);
    mask.setXmippOrigin();
    double R2 = 0.25 * Xdim * Xdim;
    FOR_ALL_ELEMENTS_IN_ARRAY2D(mask)
      A2D_ELEM(mask,i,j)=(i*i+j*j<R2);
    Npixels = (int) mask.sum();

    // Thread Manager
    thMgr = std::make_unique<ThreadManager>(Nthreads, this);
    Image<double> dummy;
    Matrix2D<double> dummyMatrix, dummyHblock, dummyW;
    dummyW.resizeNoCopy(Npixels, Neigen + 2);
    dummyHblock.resizeNoCopy(2 * Nangles * Nshifts, Neigen + 2);
    for (int n = 0; n < Nthreads; ++n)
    {
      Wnode.push_back(dummyW);
      Hblock.push_back(dummyHblock);
      A.push_back(dummyMatrix);
      I.push_back(dummy);
      Iaux.push_back(dummy());
      MD.push_back(MDin);
    }

    Matrix2D<double> &W = Wnode[0];
    FileName fnMatrixF(fnRoot + "_matrixF.raw");
    FileName fnMatrixH(fnRoot + "_matrixH.raw");

    if (IS_MASTER)
    {
      // F (#images*#shifts*#angles) x (#eigenvectors+2)*(its+1)
      fnMatrixF.createEmptyFileWithGivenLength(
          Nimg * 2 * Nangles * Nshifts * (Neigen + 2) * (Nits + 1)
              * sizeof(double));
      // H (#images*#shifts*#angles) x (#eigenvectors+2)
      fnMatrixH.createEmptyFileWithGivenLength(
          Nimg * 2 * Nangles * Nshifts * (Neigen + 2) * sizeof(double));

      // Initialize with random numbers between -1 and 1
      FOR_ALL_ELEMENTS_IN_MATRIX2D(W)
        MAT_ELEM(W,i,j)=rnd_unif(-1.0,1.0);
    }

    comunicateMatrix(W);

    F.mapToFile(fnMatrixF, (Neigen + 2) * (Nits + 1),
        Nimg * 2 * Nangles * Nshifts);
    H.mapToFile(fnMatrixH, Nimg * 2 * Nangles * Nshifts, Neigen + 2);

    // Prepare buffer
    for (int n = 0; n < HbufferMax; n++)
      Hbuffer.push_back(
          new double[MAT_XSIZE(dummyHblock) * MAT_YSIZE(dummyHblock)]);

    // Construct a FileTaskDistributor
    MDin.findObjects(objId);
    size_t Nimgs = objId.size();
    createMutexes(Nimgs);
}

QR()

int ProgImageRotationalPCA::QR

(

)

QR decomposition. In fact, only Q is computed. It returns the number of columns of Q different from 0.

Definition at line 461 of file image_rotational_pca.cpp.

{
  size_t jQ=0;
  Matrix1D<double> qj1, qj2;
  int iBlockMax=MAT_XSIZE(F)/4;

  for (size_t j1=0; j1<MAT_YSIZE(F); j1++)
  {
    F.getRow(j1,qj1);
    // Project twice in the already established subspace
    // One projection should be enough but Gram-Schmidt suffers
    // from numerical problems
    for (int it=0; it<2; it++)
    {
      for (size_t j2=0; j2<jQ; j2++)
      {
        F.getRow(j2,qj2);

        // Compute dot product
        double s12=qj1.dotProduct(qj2);

        // Subtract the part of qj2 from qj1
        double *ptr1=&VEC_ELEM(qj1,0);
        const double *ptr2=&VEC_ELEM(qj2,0);
        for (int i=0; i<iBlockMax; i++)
        {
          (*ptr1++)-=s12*(*ptr2++);
          (*ptr1++)-=s12*(*ptr2++);
          (*ptr1++)-=s12*(*ptr2++);
          (*ptr1++)-=s12*(*ptr2++);
        }
        for (size_t i=iBlockMax*4; i<MAT_XSIZE(F); ++i)
        (*ptr1++)-=s12*(*ptr2++);
      }
    }

    // Keep qj1 in Q if it has enough norm
    double Rii=qj1.module();
    if (Rii>1e-14)
    {
      // Make qj1 to be unitary and store in Q
      qj1/=Rii;
      F.setRow(jQ++,qj1);
    }
  }
  return jQ;
}

readParams()

void ProgImageRotationalPCA::readParams ( )

virtual

Read argument from command line.

Reimplemented from XmippProgram.

Definition at line 57 of file image_rotational_pca.cpp.

{
  fnIn = getParam("-i");
  fnRoot = getParam("--oroot");
  Neigen = getIntParam("--eigenvectors");
  Nits = getIntParam("--iterations");
  max_shift_change = getDoubleParam("--max_shift_change");
  psi_step = getDoubleParam("--psi_step");
  shift_step = getDoubleParam("--shift_step");
  maxNimgs = getIntParam("--maxImages");
  Nthreads = getIntParam("--thr");
}

run()

void ProgImageRotationalPCA::run ( )

virtual

Run.

Reimplemented from XmippProgram.

Definition at line 563 of file image_rotational_pca.cpp.

{
  show();
  produceSideInfo();

  // Compute matrix F:
  // Set H pointing to the first block of F
  applyTt();// H=Tt(W)
  copyHtoF(0);
  for (int it=0; it<Nits; it++)
  {
    applyT(); // W=T(H)
    applyTt();// H=Tt(W)
    copyHtoF(it+1);
  }
  H.clear();
  clearHbuffer();

  // QR decomposition of matrix F
  int qrDim;

  comunicateQrDim(qrDim);

  if (qrDim==0)
  REPORT_ERROR(ERR_VALUE_INCORRECT,"No subspace have been found");

  mapMatrix(qrDim);
  F.clear();

  // Apply T
  for (int n=0; n<Nthreads; ++n)
  Hblock[n].resizeNoCopy(2*Nangles*Nshifts,qrDim);
  applyT();

  // Free memory
  H.clear();
  Hblock.clear();

  applySVD();

  // Clean files
  if (IS_MASTER)
  {
    unlink((fnRoot+"_matrixH.raw").c_str());
    unlink((fnRoot+"_matrixF.raw").c_str());
  }
}

selectPartFromMd()

void ProgImageRotationalPCA::selectPartFromMd ( MetaData &  MDin )

virtual

Read input images

Reimplemented in MpiProgImageRotationalPCA.

Definition at line 108 of file image_rotational_pca.cpp.

{
    MetaDataVec MDaux;
    MDaux.randomize(MDin);
    MDin.selectPart(MDaux, 0, maxNimgs);
}

show()

void ProgImageRotationalPCA::show

(

)

Show.

Definition at line 72 of file image_rotational_pca.cpp.

{
  if (!verbose)
    return;
  std::cout << "Input:               " << fnIn << std::endl
      << "Output root:         " << fnRoot << std::endl
      << "Number eigenvectors: " << Neigen << std::endl
      << "Number iterations:   " << Nits << std::endl << "Psi step:            "
      << psi_step << std::endl << "Max shift change:    " << max_shift_change
      << " step: " << shift_step << std::endl << "Max images:          "
      << maxNimgs << std::endl << "Number of threads:   " << Nthreads
      << std::endl;
}

writeToHBuffer()

void ProgImageRotationalPCA::writeToHBuffer	(	int	idx,
		double *	dest
	)

Write to H buffer

Definition at line 202 of file image_rotational_pca.cpp.

{
  threadMutex->lock();
  int n=HbufferDestination.size();
  const Matrix2D<double> &Hblock_idx=Hblock[idx];
  // COSS std::cout << "Copying block " << idx << " into " << n << "(" << Hbuffer.size() << ")" << std::endl;
  memcpy(Hbuffer[n],&MAT_ELEM(Hblock_idx,0,0),MAT_XSIZE(Hblock_idx)*MAT_YSIZE(Hblock_idx)*sizeof(double));
  HbufferDestination.push_back(dest);
  if (n==(HbufferMax-1))
  flushHBuffer();
  threadMutex->unlock();
}

Member Data Documentation

A

std::vector< Matrix2D<double> > ProgImageRotationalPCA::A

Definition at line 100 of file image_rotational_pca.h.

F

Matrix2D<double> ProgImageRotationalPCA::F

Definition at line 93 of file image_rotational_pca.h.

fileMutex

std::unique_ptr<Mutex> ProgImageRotationalPCA::fileMutex

Definition at line 87 of file image_rotational_pca.h.

fnIn

FileName ProgImageRotationalPCA::fnIn

Input selfile

Definition at line 45 of file image_rotational_pca.h.

fnRoot

FileName ProgImageRotationalPCA::fnRoot

Output root

Definition at line 47 of file image_rotational_pca.h.

H

Matrix2D<double> ProgImageRotationalPCA::H

Definition at line 91 of file image_rotational_pca.h.

Hblock

std::vector< Matrix2D<double> > ProgImageRotationalPCA::Hblock

Definition at line 102 of file image_rotational_pca.h.

Hbuffer

std::vector<double *> ProgImageRotationalPCA::Hbuffer

Definition at line 81 of file image_rotational_pca.h.

HbufferDestination

std::vector<double *> ProgImageRotationalPCA::HbufferDestination

Definition at line 83 of file image_rotational_pca.h.

HbufferMax

const int ProgImageRotationalPCA::HbufferMax =20

static

Definition at line 85 of file image_rotational_pca.h.

I

std::vector< Image<double> > ProgImageRotationalPCA::I

Definition at line 96 of file image_rotational_pca.h.

Iaux

std::vector< MultidimArray<double> > ProgImageRotationalPCA::Iaux

Definition at line 98 of file image_rotational_pca.h.

mask

MultidimArray< unsigned char > ProgImageRotationalPCA::mask

Definition at line 108 of file image_rotational_pca.h.

max_shift_change

double ProgImageRotationalPCA::max_shift_change

Maximum shift change

Definition at line 55 of file image_rotational_pca.h.

maxNimgs

int ProgImageRotationalPCA::maxNimgs

Maximum number of images

Definition at line 59 of file image_rotational_pca.h.

MD

std::vector<MetaDataVec> ProgImageRotationalPCA::MD

Definition at line 67 of file image_rotational_pca.h.

Nangles

int ProgImageRotationalPCA::Nangles

Definition at line 71 of file image_rotational_pca.h.

Neigen

int ProgImageRotationalPCA::Neigen

Number of eigenvectors

Definition at line 49 of file image_rotational_pca.h.

Nimg

size_t ProgImageRotationalPCA::Nimg

Definition at line 69 of file image_rotational_pca.h.

Nits

int ProgImageRotationalPCA::Nits

Number of iterations

Definition at line 51 of file image_rotational_pca.h.

Npixels

int ProgImageRotationalPCA::Npixels

Definition at line 77 of file image_rotational_pca.h.

Nshifts

int ProgImageRotationalPCA::Nshifts

Definition at line 73 of file image_rotational_pca.h.

Nthreads

int ProgImageRotationalPCA::Nthreads

Number of threads

Definition at line 61 of file image_rotational_pca.h.

objId

std::vector<size_t> ProgImageRotationalPCA::objId

Definition at line 114 of file image_rotational_pca.h.

psi_step

double ProgImageRotationalPCA::psi_step

Psi step

Definition at line 53 of file image_rotational_pca.h.

rank

int ProgImageRotationalPCA::rank

Rank, used later for MPI

Definition at line 63 of file image_rotational_pca.h.

shift_step

double ProgImageRotationalPCA::shift_step

Shift step

Definition at line 57 of file image_rotational_pca.h.

taskDistributor

std::unique_ptr<ThreadTaskDistributor> ProgImageRotationalPCA::taskDistributor

Definition at line 110 of file image_rotational_pca.h.

thMgr

std::unique_ptr<ThreadManager> ProgImageRotationalPCA::thMgr

Definition at line 112 of file image_rotational_pca.h.

threadMutex

std::unique_ptr<Mutex> ProgImageRotationalPCA::threadMutex

Definition at line 89 of file image_rotational_pca.h.

Wnode

std::vector< Matrix2D<double> > ProgImageRotationalPCA::Wnode

Definition at line 104 of file image_rotational_pca.h.

Wtranspose

Matrix2D<double> ProgImageRotationalPCA::Wtranspose

Definition at line 106 of file image_rotational_pca.h.

Xdim

size_t ProgImageRotationalPCA::Xdim

Definition at line 75 of file image_rotational_pca.h.

---

The documentation for this class was generated from the following files:

• xmipp/libraries/reconstruction/image_rotational_pca.h
• xmipp/libraries/reconstruction/image_rotational_pca.cpp