1c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath//===================================================== 2c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// File : init_matrix.hh 3c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// Author : L. Plagne <laurent.plagne@edf.fr)> 4c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// Copyright (C) EDF R&D, lun sep 30 14:23:19 CEST 2002 5c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath//===================================================== 6c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// 7c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// This program is free software; you can redistribute it and/or 8c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// modify it under the terms of the GNU General Public License 9c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// as published by the Free Software Foundation; either version 2 10c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// of the License, or (at your option) any later version. 11c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// 12c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// This program is distributed in the hope that it will be useful, 13c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// but WITHOUT ANY WARRANTY; without even the implied warranty of 14c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the 15c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// GNU General Public License for more details. 16c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// You should have received a copy of the GNU General Public License 17c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// along with this program; if not, write to the Free Software 18c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA. 19c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// 20c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath#ifndef INIT_MATRIX_HH 21c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath#define INIT_MATRIX_HH 22c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 23c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// The Vector class must satisfy the following part of STL vector concept : 24c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// resize() method 25c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// [] operator for setting element 26c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// value_type defined 27c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamathtemplate<double init_function(int,int), class Vector> 28c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan KamathBTL_DONT_INLINE void init_row(Vector & X, int size, int row){ 29c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 30c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath X.resize(size); 31c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 32c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath for (int j=0;j<X.size();j++){ 33c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath X[j]=typename Vector::value_type(init_function(row,j)); 34c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 35c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath} 36c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 37c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 38c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// Matrix is a Vector of Vector 39c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// The Matrix class must satisfy the following part of STL vector concept : 40c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// resize() method 41c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// [] operator for setting rows 42c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamathtemplate<double init_function(int,int),class Vector> 43c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan KamathBTL_DONT_INLINE void init_matrix(Vector & A, int size){ 44c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath A.resize(size); 45c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath for (int row=0; row<A.size() ; row++){ 46c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath init_row<init_function>(A[row],size,row); 47c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 48c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath} 49c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 50c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamathtemplate<double init_function(int,int),class Matrix> 51c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan KamathBTL_DONT_INLINE void init_matrix_symm(Matrix& A, int size){ 52c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath A.resize(size); 53c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath for (int row=0; row<A.size() ; row++) 54c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath A[row].resize(size); 55c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath for (int row=0; row<A.size() ; row++){ 56c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath A[row][row] = init_function(row,row); 57c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath for (int col=0; col<row ; col++){ 58c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath double x = init_function(row,col); 59c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath A[row][col] = A[col][row] = x; 60c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 61c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 62c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath} 63c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 64c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath#endif 65