1c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath//===================================================== 2c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// File : STL_interface.hh 3c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// Author : L. Plagne <laurent.plagne@edf.fr)> 4c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// Copyright (C) EDF R&D, lun sep 30 14:23:24 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 STL_INTERFACE_HH 21c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath#define STL_INTERFACE_HH 22c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath#include <string> 23c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath#include <vector> 24c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath#include "utilities.h" 25c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 26c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamathusing namespace std; 27c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 28c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamathtemplate<class real> 29c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamathclass STL_interface{ 30c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 31c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamathpublic : 32c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 33c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath typedef real real_type ; 34c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 35c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath typedef std::vector<real> stl_vector; 36c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath typedef std::vector<stl_vector > stl_matrix; 37c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 38c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath typedef stl_matrix gene_matrix; 39c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 40c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath typedef stl_vector gene_vector; 41c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 42c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath static inline std::string name( void ) 43c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath { 44c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath return "STL"; 45c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 46c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 47c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath static void free_matrix(gene_matrix & A, int N){} 48c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 49c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath static void free_vector(gene_vector & B){} 50c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 51c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath static inline void matrix_from_stl(gene_matrix & A, stl_matrix & A_stl){ 52c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath A = A_stl; 53c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 54c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 55c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath static inline void vector_from_stl(gene_vector & B, stl_vector & B_stl){ 56c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath B = B_stl; 57c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 58c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 59c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath static inline void vector_to_stl(gene_vector & B, stl_vector & B_stl){ 60c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath B_stl = B ; 61c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 62c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 63c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 64c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath static inline void matrix_to_stl(gene_matrix & A, stl_matrix & A_stl){ 65c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath A_stl = A ; 66c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 67c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 68c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath static inline void copy_vector(const gene_vector & source, gene_vector & cible, int N){ 69c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath for (int i=0;i<N;i++){ 70c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath cible[i]=source[i]; 71c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 72c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 73c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 74c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 75c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath static inline void copy_matrix(const gene_matrix & source, gene_matrix & cible, int N){ 76c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath for (int i=0;i<N;i++) 77c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath for (int j=0;j<N;j++) 78c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath cible[i][j]=source[i][j]; 79c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 80c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 81c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// static inline void ata_product(const gene_matrix & A, gene_matrix & X, int N) 82c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// { 83c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// real somme; 84c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// for (int j=0;j<N;j++){ 85c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// for (int i=0;i<N;i++){ 86c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// somme=0.0; 87c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// for (int k=0;k<N;k++) 88c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// somme += A[i][k]*A[j][k]; 89c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// X[j][i]=somme; 90c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// } 91c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// } 92c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath// } 93c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 94c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath static inline void aat_product(const gene_matrix & A, gene_matrix & X, int N) 95c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath { 96c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath real somme; 97c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath for (int j=0;j<N;j++){ 98c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath for (int i=0;i<N;i++){ 99c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath somme=0.0; 100c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath if(i>=j) 101c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath { 102c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath for (int k=0;k<N;k++){ 103c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath somme+=A[k][i]*A[k][j]; 104c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 105c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath X[j][i]=somme; 106c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 107c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 108c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 109c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 110c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 111c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 112c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath static inline void matrix_matrix_product(const gene_matrix & A, const gene_matrix & B, gene_matrix & X, int N) 113c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath { 114c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath real somme; 115c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath for (int j=0;j<N;j++){ 116c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath for (int i=0;i<N;i++){ 117c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath somme=0.0; 118c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath for (int k=0;k<N;k++) 119c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath somme+=A[k][i]*B[j][k]; 120c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath X[j][i]=somme; 121c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 122c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 123c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 124c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 125c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath static inline void matrix_vector_product(gene_matrix & A, gene_vector & B, gene_vector & X, int N) 126c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath { 127c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath real somme; 128c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath for (int i=0;i<N;i++){ 129c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath somme=0.0; 130c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath for (int j=0;j<N;j++) 131c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath somme+=A[j][i]*B[j]; 132c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath X[i]=somme; 133c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 134c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 135c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 136c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath static inline void symv(gene_matrix & A, gene_vector & B, gene_vector & X, int N) 137c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath { 138c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath for (int j=0; j<N; ++j) 139c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath X[j] = 0; 140c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath for (int j=0; j<N; ++j) 141c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath { 142c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath real t1 = B[j]; 143c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath real t2 = 0; 144c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath X[j] += t1 * A[j][j]; 145c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath for (int i=j+1; i<N; ++i) { 146c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath X[i] += t1 * A[j][i]; 147c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath t2 += A[j][i] * B[i]; 148c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 149c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath X[j] += t2; 150c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 151c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 152c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 153c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath static inline void syr2(gene_matrix & A, gene_vector & B, gene_vector & X, int N) 154c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath { 155c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath for (int j=0; j<N; ++j) 156c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath { 157c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath for (int i=j; i<N; ++i) 158c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath A[j][i] += B[i]*X[j] + B[j]*X[i]; 159c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 160c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 161c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 162c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath static inline void ger(gene_matrix & A, gene_vector & X, gene_vector & Y, int N) 163c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath { 164c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath for (int j=0; j<N; ++j) 165c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath { 166c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath for (int i=j; i<N; ++i) 167c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath A[j][i] += X[i]*Y[j]; 168c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 169c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 170c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 171c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath static inline void atv_product(gene_matrix & A, gene_vector & B, gene_vector & X, int N) 172c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath { 173c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath real somme; 174c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath for (int i=0;i<N;i++){ 175c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath somme = 0.0; 176c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath for (int j=0;j<N;j++) 177c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath somme += A[i][j]*B[j]; 178c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath X[i] = somme; 179c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 180c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 181c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 182c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath static inline void axpy(real coef, const gene_vector & X, gene_vector & Y, int N){ 183c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath for (int i=0;i<N;i++) 184c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath Y[i]+=coef*X[i]; 185c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 186c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 187c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath static inline void axpby(real a, const gene_vector & X, real b, gene_vector & Y, int N){ 188c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath for (int i=0;i<N;i++) 189c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath Y[i] = a*X[i] + b*Y[i]; 190c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 191c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 192c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath static inline void trisolve_lower(const gene_matrix & L, const gene_vector & B, gene_vector & X, int N){ 193c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath copy_vector(B,X,N); 194c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath for(int i=0; i<N; ++i) 195c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath { 196c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath X[i] /= L[i][i]; 197c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath real tmp = X[i]; 198c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath for (int j=i+1; j<N; ++j) 199c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath X[j] -= tmp * L[i][j]; 200c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 201c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 202c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 203c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath static inline real norm_diff(const stl_vector & A, const stl_vector & B) 204c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath { 205c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath int N=A.size(); 206c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath real somme=0.0; 207c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath real somme2=0.0; 208c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 209c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath for (int i=0;i<N;i++){ 210c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath real diff=A[i]-B[i]; 211c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath somme+=diff*diff; 212c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath somme2+=A[i]*A[i]; 213c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 214c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath return somme/somme2; 215c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 216c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 217c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath static inline real norm_diff(const stl_matrix & A, const stl_matrix & B) 218c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath { 219c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath int N=A[0].size(); 220c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath real somme=0.0; 221c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath real somme2=0.0; 222c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 223c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath for (int i=0;i<N;i++){ 224c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath for (int j=0;j<N;j++){ 225c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath real diff=A[i][j] - B[i][j]; 226c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath somme += diff*diff; 227c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath somme2 += A[i][j]*A[i][j]; 228c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 229c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 230c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 231c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath return somme/somme2; 232c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 233c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 234c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath static inline void display_vector(const stl_vector & A) 235c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath { 236c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath int N=A.size(); 237c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath for (int i=0;i<N;i++){ 238c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath INFOS("A["<<i<<"]="<<A[i]<<endl); 239c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 240c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath } 241c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 242c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath}; 243c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath 244c981c48f5bc9aefeffc0bcb0cc3934c2fae179ddNarayan Kamath#endif 245