Changeset 2583 in Sophya for trunk/SophyaLib/TArray/basarr.cc

Timestamp:

Jul 30, 2004, 12:24:12 PM (21 years ago)

Author:

cmv

Message:

• Intro decision auto d'optimisation produit de matrices
• Possibilite a l'utilisateur pour choisir l'optimisation
• cas FxC optimise par copie

cmv 30/07/04

File:

• trunk/SophyaLib/TArray/basarr.cc (modified) (1 diff)

trunk/SophyaLib/TArray/basarr.cc

@@ -29 +29 @@
 short  BaseArray::default_vector_type = ColumnVector;
 sa_size_t BaseArray::openmp_size_threshold = 200000;
+uint_2 BaseArray::matrix_product_optim = 1;
 
 // ------ Methodes statiques globales --------
+
+//! Set optimization strategy for matrix product
+/*!
+  \param opt : bit coded optimization strategy
+  \warning Default is opt=1
+  \verbatim
+  bit 0 : choose matrix product optimization or not
+        0=no optimization, 1=optimization
+  bit 1 : force optimization in any case (only if bit0=1)
+          (if not the TMatrix::Multiply method decide or
+           not if the product should be optimized)
+        0=do not force and let the method decide, 1=force
+  bit 2 : optimize the product A * B when A-columns and
+          B-rows are not packed (for example if the product
+          is C = ("A Fortran-like" * "B C-like").
+        . That will do a copy of one of the two matrices,
+          so that will result in an increase of the memory
+          space needed.
+        . For big matrices that decrease the computing time.
+        . Do not use this optimisation for small matrices
+          because that would increase the computing time.
+        0=do not optimze that way, 1=optimze that way
+  \endverbatim
+  \verbatim
+   Sumary of the allowed values for "opt"
+           0 = no optimization at all (whatever the other bits are)
+           1 = optimize but let the method decides if optimization
+               is needed regarding the sizes of the matrices.
+           3 = force optimization whatever the sizes of the matrices are.
+           5 = optimisation with method decision ("1")
+               AND optimize by copying when A-columns and B-rows
+               are not packed
+           7 = force optimization ("3")
+               AND optimize by copying when A-columns and B-rows
+               are not packed
+  \endverbatim
+*/
+void BaseArray::SetMatProdOpt(uint_2 opt)
+{
+  matrix_product_optim = opt;
+}
 
 //! Set maximum number of printed elements and print level