[Metil-lmtpp] probleme avec -O3

[Felipe B. W. <fel...@gm...>]

le même code compilé avec -O0 et -O3:

Vec<double> U ;
U.resize(indice_i.size(),0);

b = Get_K()(indice_b, indice_i)*U;

PRINT(Get_K().nb_rows());       ///   2385
PRINT(Get_K().nb_cols());        ///  2385
PRINT(max(indice_b));              /// 2384
PRINT(max(indice_i));              /// 2210
PRINT(b);

avec -O0:

b -> 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ...
correct :)

avec -O3:
b ->  -2.09192e-316 -2.09192e-316 -2.09192e-316 -2.09192e-316  
-2.09192e-316 -2.09192e-316 -2.09192e-316 ...

et à chaque fois il donne une valeur différente.
tres bizarre si U et vide [0 0 0 0 0 0 0...] comment on peut avoir un  
résultat non nul??


---------
avec:

Mat<double,Sym<>,SparseLine<> >& Get_K() { return pb->K;};
et K une matrice élément finit bien définit

et pb :

FormulationAssembly<PB> *pb;

----

l'opérateur Mat * Vec a un bug????

------------------------------------------------------------------------
Felipe Bordeu Weldt
Doctorant - 3° année   (3° year PhD student)
------------------------------------------------------------------------
LMT-Cachan
(ENS Cachan/CNRS/UPMC/PRES UniverSud Paris)
61, av. du Président Wilson
F-94230 Cachan
France
------------------------------------------------------------------------