escript/src/LocalOps.h

// $Id$
/* 
 ******************************************************************************
 *                                                                            *
 *       COPYRIGHT  ACcESS 2004 -  All Rights Reserved                        *
 *                                                                            *
 * This software is the property of ACcESS. No part of this code              *
 * may be copied in any form or by any means without the expressed written    *
 * consent of ACcESS.  Copying, use or modification of this software          *
 * by any unauthorised person is illegal unless that person has a software    *
 * license agreement with ACcESS.                                             *
 *                                                                            *
 ******************************************************************************
*/

#if !defined escript_LocalOps_H
#define escript_LocalOps_H

namespace escript {


/**
   \brief
   solves a 1x1 eigenvalue A*V=ev*V problem

   \param A00 Input - A_00 
   \param ev0 Output - eigenvalue
*/
inline
void eigenvalues1(const double A00,double* ev0) {

   *ev0=1.;

}
/**
   \brief
   solves a 2x2 eigenvalue A*V=ev*V problem for symmetric A

   \param A00 Input - A_00 
   \param A01 Input - A_01 
   \param A11 Input - A_11
   \param ev0 Output - smallest eigenvalue
   \param ev1 Output - largest eigenvalue
*/
inline
void eigenvalues2(const double A00,const double A01
                                 ,const double A11,
                 double* ev0, double* ev1) {

   *ev0=1.;
   *ev1=2.;
}
/**
   \brief
   solves a 3x3 eigenvalue A*V=ev*V problem for symmetric A

   \param A00 Input - A_00 
   \param A01 Input - A_01 
   \param A02 Input - A_02 
   \param A11 Input - A_11 
   \param A12 Input - A_12 
   \param A22 Input - A_22 
   \param ev0 Output - smallest eigenvalue
   \param ev1 Output - eigenvalue
   \param ev2 Output - largest eigenvalue
*/
inline
void eigenvalues3(const double A00, const double A01, const double A02,
                                   const double A11, const double A12,
                                                     const double A22,
                 double* ev0, double* ev1,double* ev2) {

   *ev0=1.;
   *ev1=2.;
   *ev2=3.;

}

} // end of namespace
#endif
1	gross	576	// $Id$
2			/*
3			******************************************************************************
4			* *
5			* COPYRIGHT ACcESS 2004 - All Rights Reserved *
6			* *
7			* This software is the property of ACcESS. No part of this code *
8			* may be copied in any form or by any means without the expressed written *
9			* consent of ACcESS. Copying, use or modification of this software *
10			* by any unauthorised person is illegal unless that person has a software *
11			* license agreement with ACcESS. *
12			* *
13			******************************************************************************
14			*/
15
16			#if !defined escript_LocalOps_H
17			#define escript_LocalOps_H
18
19			namespace escript {
20
21
22			/**
23			\brief
24			solves a 1x1 eigenvalue AV=evV problem
25
26			\param A00 Input - A_00
27			\param ev0 Output - eigenvalue
28			*/
29			inline
30			void eigenvalues1(const double A00,double* ev0) {
31
32			*ev0=1.;
33
34			}
35			/**
36			\brief
37			solves a 2x2 eigenvalue AV=evV problem for symmetric A
38
39			\param A00 Input - A_00
40			\param A01 Input - A_01
41			\param A11 Input - A_11
42			\param ev0 Output - smallest eigenvalue
43			\param ev1 Output - largest eigenvalue
44			*/
45			inline
46			void eigenvalues2(const double A00,const double A01
47			,const double A11,
48			double* ev0, double* ev1) {
49
50			*ev0=1.;
51			*ev1=2.;
52			}
53			/**
54			\brief
55			solves a 3x3 eigenvalue AV=evV problem for symmetric A
56
57			\param A00 Input - A_00
58			\param A01 Input - A_01
59			\param A02 Input - A_02
60			\param A11 Input - A_11
61			\param A12 Input - A_12
62			\param A22 Input - A_22
63			\param ev0 Output - smallest eigenvalue
64			\param ev1 Output - eigenvalue
65			\param ev2 Output - largest eigenvalue
66			*/
67			inline
68			void eigenvalues3(const double A00, const double A01, const double A02,
69			const double A11, const double A12,
70			const double A22,
71			double* ev0, double* ev1,double* ev2) {
72
73			*ev0=1.;
74			*ev1=2.;
75			*ev2=3.;
76
77			}
78
79			} // end of namespace
80			#endif