MatrixTriangle Class Reference

#include <MatrixTriangle.h>

Collaboration diagram for MatrixTriangle:

Classes
struct	MatrixTriangleDataTag

Public Member Functions
	MatrixTriangle (int _n=0)
void	setSize (int _n)
	~MatrixTriangle ()
	MatrixTriangle (const MatrixTriangle &A)
MatrixTriangle &	operator= (const MatrixTriangle &A)
MatrixTriangle	clone ()
void	copyFrom (MatrixTriangle r)
bool	operator== (const MatrixTriangle &A)
int	nLine ()
double *	operator[] (int i)
	operator double ** () const
void	solveInPlace (Vector b)
void	solveTransposInPlace (Vector y)
void	LINPACK (Vector &u)

Static Public Attributes
static MatrixTriangle	emptyMatrixTriangle

Protected Types
typedef struct MatrixTriangle::MatrixTriangleDataTag	MatrixTriangleData

Protected Member Functions
void	destroyCurrentBuffer ()

Protected Attributes
MatrixTriangleData *	d

Friends
class	Matrix

Detailed Description

Definition at line 34 of file MatrixTriangle.h.

Member Typedef Documentation

MatrixTriangleData

typedef struct MatrixTriangle::MatrixTriangleDataTag MatrixTriangle::MatrixTriangleData

protected

Constructor & Destructor Documentation

MatrixTriangle() [1/2]

MatrixTriangle::MatrixTriangle

(

int

_n = 0

)

Definition at line 31 of file MatrixTriangle.cpp.

{
    d=(MatrixTriangleData*)malloc(sizeof(MatrixTriangleData));
    d->n=_n; d->ext=_n;
    d->ref_count=1;
    if (_n>0)
    {
        double **t,*t2;
        int i=1;
        t=d->p=(double**)malloc(_n*sizeof(double*));
        t2=(double*)malloc((_n+1)*_n/2*sizeof(double));
        while(_n--)
        {
            *(t++)=t2; t2+=i; i++;
        }
    } else d->p=NULL;
}

~MatrixTriangle()

MatrixTriangle::~MatrixTriangle

(

)

Definition at line 137 of file MatrixTriangle.cpp.

{
    destroyCurrentBuffer();
}

MatrixTriangle() [2/2]

MatrixTriangle::MatrixTriangle

(

const MatrixTriangle &

A

)

Definition at line 153 of file MatrixTriangle.cpp.

{
    // shallow copy
    d=A.d;
    (d->ref_count)++ ;
}

Member Function Documentation

clone()

MatrixTriangle MatrixTriangle::clone

(

)

Definition at line 172 of file MatrixTriangle.cpp.

{
    // a deep copy
    MatrixTriangle r(nLine());
    r.copyFrom(*this);
    return r;
}

copyFrom()

void MatrixTriangle::copyFrom

(

MatrixTriangle

r

)

Definition at line 180 of file MatrixTriangle.cpp.

{
    int n=r.nLine();
    setSize(n);
    if (n==0) return;
    memcpy(*d->p,*(r.d->p),(n+1)*n/2*sizeof(double));
}

destroyCurrentBuffer()

void MatrixTriangle::destroyCurrentBuffer ( )

protected

Definition at line 142 of file MatrixTriangle.cpp.

{
    if (!d) return;
    (d->ref_count) --;
    if (d->ref_count==0)
    {
        if (d->p) { free(*d->p); free(d->p); }
        free(d);
    };
}

LINPACK()

void MatrixTriangle::LINPACK

(

Vector &

u

)

Definition at line 119 of file MatrixTriangle.cpp.

{
    int i,j,n=nLine();
    R.setSize(n);
    double **L=(*this), *w=R, sum;

    for (i=0; i<n; i++)
    {
        if (L[i][i]==0) w[i]=1.0;

        sum=0; j=i-1;
        if (i) while (j--) sum+=L[i][j]*w[j];
        if (((1.0-sum)/L[i][i])>((-1.0-sum)/L[i][i])) w[i]=1.0; else w[i]=-1.0;
    }
    solveTransposInPlace(R);
    R.multiply(1/R.euclidianNorm());
}

nLine()

int MatrixTriangle::nLine ( )

inline

Definition at line 63 of file MatrixTriangle.h.

{ return d->n; };

operator double **()

MatrixTriangle::operator double ** ( ) const

inline

Definition at line 65 of file MatrixTriangle.h.

{ return d->p; };

operator=()

MatrixTriangle & MatrixTriangle::operator=

(

const MatrixTriangle &

A

)

Definition at line 160 of file MatrixTriangle.cpp.

{
    // shallow copy
    if (this != &A)
    {
        destroyCurrentBuffer();
        d=A.d;
        (d->ref_count) ++ ;
    }
    return *this;
}

operator==()

bool MatrixTriangle::operator== ( const MatrixTriangle &  A )

inline

Definition at line 62 of file MatrixTriangle.h.

{ return (A.d==d);}

operator[]()

double* MatrixTriangle::operator[] ( int  i )

inline

Definition at line 64 of file MatrixTriangle.h.

{ return d->p[i]; };

setSize()

void MatrixTriangle::setSize

(

int

_n

)

Definition at line 49 of file MatrixTriangle.cpp.

{
    d->n=_n;
    if (_n>d->ext)
    {
        d->ext=_n;
        double **t,*t2;
        if (!d->p)
        {
            t2=(double*)malloc((_n+1)*_n/2*sizeof(double));
            t=d->p=(double**)malloc(_n*sizeof(double));
        } else
        {
            t2=(double*)realloc(*d->p,(_n+1)*_n/2*sizeof(double));
            t=d->p=(double**)realloc(d->p,_n*sizeof(double));
        }
        int i=1;
        while(_n--)
        {
            *(t++)=t2; t2+=i; i++;
        }
    }
}

solveInPlace()

void MatrixTriangle::solveInPlace

(

Vector

b

)

Definition at line 73 of file MatrixTriangle.cpp.

{
    int i,k,n=nLine();
    double **a=(*this), *x=b, sum;

    if ((int)b.sz()!=n)
    {
        printf("error in matrixtriangle solve.\n"); getchar(); exit(254);
    }
    for (i=0; i<n; i++)
    {
        sum=x[i]; k=i;
        while (k--) sum-=a[i][k]*x[k];
        x[i]=sum/a[i][i];
    }
}

solveTransposInPlace()

void MatrixTriangle::solveTransposInPlace

(

Vector

y

)

Definition at line 90 of file MatrixTriangle.cpp.

{
    int n=nLine(),i=n,k;
    double **a=(*this), *x=y, sum;

    while(i--)
    {
        sum=x[i];
        for (k=i+1; k<n; k++) sum-=a[k][i]*x[k];
        x[i]=sum/a[i][i];
    }
}

Friends And Related Function Documentation

Matrix

friend class Matrix

friend

Definition at line 36 of file MatrixTriangle.h.

Member Data Documentation

d

MatrixTriangleData* MatrixTriangle::d

protected

Definition at line 46 of file MatrixTriangle.h.

emptyMatrixTriangle

MatrixTriangle MatrixTriangle::emptyMatrixTriangle

static

Definition at line 74 of file MatrixTriangle.h.

---

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

• xmipp/external/condor/MatrixTriangle.h
• xmipp/external/condor/Matrix.cpp
• xmipp/external/condor/MatrixTriangle.cpp