[Metil-lmtpp] probleme avec -O3
Status: Pre-Alpha
Brought to you by:
hugo_lec
From: Felipe B. W. <fel...@gm...> - 2009-07-04 14:16:15
|
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 ------------------------------------------------------------------------ |