aboutsummaryrefslogtreecommitdiffhomepage
path: root/bench/btl/libs/STL/STL_interface.hh
diff options
context:
space:
mode:
authorGravatar Gael Guennebaud <g.gael@free.fr>2009-08-04 11:30:33 +0200
committerGravatar Gael Guennebaud <g.gael@free.fr>2009-08-04 11:30:33 +0200
commitc2a92e92a6e5aa0bcbb465220a8beb0687c1449a (patch)
treeb3d615aa6875c0d2bdc3778babd42f9e9b1456aa /bench/btl/libs/STL/STL_interface.hh
parent4bec101470091aae27e7469025c80e31b889f566 (diff)
add ger and lu with partial pivoting in BTL
Diffstat (limited to 'bench/btl/libs/STL/STL_interface.hh')
-rw-r--r--bench/btl/libs/STL/STL_interface.hh9
1 files changed, 9 insertions, 0 deletions
diff --git a/bench/btl/libs/STL/STL_interface.hh b/bench/btl/libs/STL/STL_interface.hh
index 3958d4af5..0b73382f3 100644
--- a/bench/btl/libs/STL/STL_interface.hh
+++ b/bench/btl/libs/STL/STL_interface.hh
@@ -156,6 +156,15 @@ public :
}
}
+ static inline void ger(gene_matrix & A, gene_vector & X, gene_vector & Y, int N)
+ {
+ for (int j=0; j<N; ++j)
+ {
+ for (int i=j; i<N; ++i)
+ A[j][i] += X[i]*Y[j];
+ }
+ }
+
static inline void atv_product(gene_matrix & A, gene_vector & B, gene_vector & X, int N)
{
real somme;