EigenVector.cpp File Reference

#include "FirstIncludes.h"
#include <stdlib.h>
#include <stdio.h>
#include <cmath>
#include <iostream>
#include <vector>
#include "MemoryDebug.h"
#include "KKBaseTypes.h"
#include "EigenVector.h"
#include "KKException.h"
#include "OSservices.h"

Go to the source code of this file.

Macros
#define	SQR(X)   ((X) * (X))

Functions
double	pythag (const double a, const double b)
template<class T >
const T	SIGN (const T &a, const T &b)

Macro Definition Documentation

#define SQR

(

X

)

((X) * (X))

Definition at line 166 of file EigenVector.cpp.

Function Documentation

double pythag	(	const double	a,
		const double	b
	)

Definition at line 170 of file EigenVector.cpp.

Referenced by KKB::tqli().

{
  double  absa, absb;

  absa = fabs (a);
  absb = fabs (b);

  if  (absa > absb)
  {
    return  absa * sqrt (1.0 + SQR (absb / absa));
  }
  else
  {
    if  (absb == 0.0)
      return 0.0;
    else 
      return absb * sqrt (1.0 + SQR (absa / absb));
  }
}  /* pythag */

template<class T >

const T SIGN ( const T &  a, const T &  b  )

inline

Definition at line 144 of file EigenVector.cpp.

{
  if  (b >= 0)
  {
    if  (a >= 0)
      return  a;
    else
      return -a;
  }

  else
  {
    if  (a >- 0)
      return -a;
    else
      return  a;
  }
}  /* SIGN */