ElecDiagonalTensor.h File Reference

#include "NumType.h"
#include "Vector.h"
#include <sstream>
#include "ElecOffDiagonalTensor.h"

Include dependency graph for ElecDiagonalTensor.h:

This graph shows which files directly or indirectly include this file:

Go to the source code of this file.

Classes
class	ElecDiagonalTensor< T, N >

Functions
template<class T , int N>
ostream &	operator<< (ostream &os, const ElecDiagonalTensor< T, N > &v)
template<class T , int N>
ElecDiagonalTensor< T, N >	operator+ (const ElecDiagonalTensor< T, N > &a, const ElecDiagonalTensor< T, N > &b)
template<class T , int N>
ElecDiagonalTensor< T, N >	operator- (const ElecDiagonalTensor< T, N > &a, const ElecDiagonalTensor< T, N > &b)
template<class T , int N>
ElecDiagonalTensor< T, N >	operator- (const ElecDiagonalTensor< T, N > &a)
template<class T , int N>
ElecDiagonalTensor< T, N >	operator* (const ElecDiagonalTensor< T, N > &a, const ElecDiagonalTensor< T, N > &b)
template<class T , int N>
ElecDiagonalTensor< T, N >	operator* (const T s, const ElecDiagonalTensor< T, N > &a)
template<class T , int N>
ElecDiagonalTensor< T, N >	operator* (const ElecDiagonalTensor< T, N > &a, const T s)
template<class T , int N>
Vector< T, N+1 >	operator* (const ElecDiagonalTensor< T, N > &a, const Vector< T, N+1 > &b)
template<class T , int N>
ElecDiagonalTensor< T, N >	operator/ (const ElecDiagonalTensor< T, N > &a, const T s)
template<class T , int N>
Vector< T, N+1 >	operator/ (const Vector< T, N+1 > &a, const ElecDiagonalTensor< T, N > &b)
template<class T , int N>
ElecDiagonalTensor< T, N >	operator/ (const ElecDiagonalTensor< T, N > &a, const ElecDiagonalTensor< T, N > &b)
template<class T , int N>
ElecDiagonalTensor< T, N >	operator/ (const T s, const ElecDiagonalTensor< T, N > &a)
template<class T , int N>
T	DiagToOffDiag (const ElecDiagonalTensor< T, N > &x)

Function Documentation

template<class T , int N>

T DiagToOffDiag

(

const ElecDiagonalTensor< T, N > &

x

)

Definition at line 409 of file ElecDiagonalTensor.h.

{
  throw CException("diag to offdiag not defined");
}

template<class T , int N>

ElecDiagonalTensor<T,N> operator*	(	const ElecDiagonalTensor< T, N > &	a,
		const ElecDiagonalTensor< T, N > &	b
	)

Definition at line 320 of file ElecDiagonalTensor.h.

{
  throw CException("no operator defined for diag * diag");
}

template<class T , int N>

ElecDiagonalTensor<T,N> operator*	(	const T	s,
		const ElecDiagonalTensor< T, N > &	a
	)

Definition at line 327 of file ElecDiagonalTensor.h.

{
  return ElecDiagonalTensor<T,N>(a) *= s;
}

template<class T , int N>

ElecDiagonalTensor<T,N> operator*	(	const ElecDiagonalTensor< T, N > &	a,
		const T	s
	)

Definition at line 334 of file ElecDiagonalTensor.h.

{
  return ElecDiagonalTensor<T,N>(a) *= s;
}

template<class T , int N>

Vector<T,N+1> operator*	(	const ElecDiagonalTensor< T, N > &	a,
		const Vector< T, N+1 > &	b
	)

Definition at line 341 of file ElecDiagonalTensor.h.

{
  Vector<T,N+1> r;
  r = 0;
  for(int i=0; i<N; i++){
    r[i] = a[i]*b[i]+a[N+i]*b[N];
  }
  for(int i=0; i<=N; i++)
    r[N] += a[2*N+i]*b[i]; 
  
  return r;
}

template<class T , int N>

ElecDiagonalTensor<T,N> operator+	(	const ElecDiagonalTensor< T, N > &	a,
		const ElecDiagonalTensor< T, N > &	b
	)

Definition at line 299 of file ElecDiagonalTensor.h.

{
  return ElecDiagonalTensor<T,N>(a) += b;
}

template<class T , int N>

ElecDiagonalTensor<T,N> operator-	(	const ElecDiagonalTensor< T, N > &	a,
		const ElecDiagonalTensor< T, N > &	b
	)

Definition at line 306 of file ElecDiagonalTensor.h.

{
  return ElecDiagonalTensor<T,N>(a) -= b;
}

template<class T , int N>

ElecDiagonalTensor<T,N> operator-

(

const ElecDiagonalTensor< T, N > &

a

)

Definition at line 313 of file ElecDiagonalTensor.h.

{
  return -ElecDiagonalTensor<T,N>(a);
}

template<class T , int N>

ElecDiagonalTensor<T,N> operator/	(	const ElecDiagonalTensor< T, N > &	a,
		const T	s
	)

Definition at line 356 of file ElecDiagonalTensor.h.

{
  return ElecDiagonalTensor<T,N>(a) /= s;
}

template<class T , int N>

Vector<T,N+1> operator/	(	const Vector< T, N+1 > &	a,
		const ElecDiagonalTensor< T, N > &	b
	)

Definition at line 363 of file ElecDiagonalTensor.h.

{
  //solve bx = a and return x
  //rename Mx = v where M=b v=a
  Vector<T,N+1> x;
  Vector<T,N+1> v;
  ElecDiagonalTensor<T,N> M;
  x = 0;
  M = b;
  v = a;
  
  //first step: eleminate M to upper diagonal  
  for(int i=0; i<N;i++){
    if(M[2*N+i] !=0.0 && M[i] != 0.0){
      M[3*N] -= M[N+i] * (M[2*N+i]/M[i]);
      v[N] -= v[i] * (M[2*N+i]/M[i]);
    }
  }
  
  //second step: go backward to solve x
  if (M[3*N] != 0.0)
    x[N] = v[N] / M[3*N];
  for(int i=0; i<N; i++){
    if (M[i] != 0.0)
      x[i] = (v[i] - M[N+i]*x[N]) / M[i];
  }
    
  return x;
}

template<class T , int N>

ElecDiagonalTensor<T,N> operator/	(	const ElecDiagonalTensor< T, N > &	a,
		const ElecDiagonalTensor< T, N > &	b
	)

Definition at line 396 of file ElecDiagonalTensor.h.

{
  throw CException("operator not defined for diag/diag"); 
}

template<class T , int N>

ElecDiagonalTensor<T,N> operator/	(	const T	s,
		const ElecDiagonalTensor< T, N > &	a
	)

Definition at line 403 of file ElecDiagonalTensor.h.

{
  throw CException("operator not defined for s/diag");
}

template<class T , int N>

ostream& operator<< ( ostream &  os, const ElecDiagonalTensor< T, N > &  v  )

inline

Definition at line 289 of file ElecDiagonalTensor.h.

References ElecDiagonalTensor< T, N >::printFromC().

{
  v.printFromC(os);
  return os;
}