C++ : mettre un tableau 2D dans un fichier

[julien_4230]

Bonjour,

Je n'arrive absolument pas à stocker un tableau 2D dans un fichier en C++. Voici mon sous-script : mon tableau (" donnees ") est un tableau 100*100

pff = fopen ("data/matrice.dat", "wt");
for (i=0;i<100;i++) {
fwrite (donnees[i],sizeof(int),100,pff);
}
fclose (pff);

Sinon j'ai essayé cela :

pff = fopen ("data/matrice.dat", "wt");
for (i=0;i<100;i++) {
for (j=0;j<100;j++) {
fwrite (&donnees[i][j],sizeof(int),1,pff);
}
}
fclose (pff);

Mais une fois le fichier téléchargé dans R, j'obtiens n'importe quoi : autant de colonnes que de chiffres...

Je ne sais vraiment pas comment faire...

Merci de votre aide !
-----

[julien_4230]

Pardon pour le deuxième essaie c'est plutôt cela que j'ai essayé :

pff = fopen ("data/matrice.dat", "wt");
for (i=0;i<nrow(donnees);i++) {
for (j=0;j<ncol(donnees);j++) {
if(j!=(ncol(donnees)-1)) {fprintf (pff, "%d ", donnees[i][j]);}
else {fprintf (pff, "%d\n", donnees[i][j]);}
}
}
fclose (pff);

[julien_4230]

Et je me rends compte que ce dernier essaie est bon... mais n'y a-t-il pas plus optimal ???

[kwariz]

Bonjour,

ça dépend ce que tu entends pas optimisé. La portion que tu donnes (PENSE AUX BALISES CODE !!!) est un code des plus classiques en C, un code C++ utiliserait plus les iostream sans doute.
Si par optimisé tu entends «obtenir la même chose en tapant moins de code» alors la réponse est non en C. Tu peux juste éviter ton if pour le saut à ligne :

pff = fopen ("data/matrice.dat", "wt");
for (i=0;i<nrow(donnees);i++) {
 for (j=0;j<ncol(donnees);j++)
   fprintf (pff, "%d ", donnees[i][j]);
 fprintf(pff,"\n");
}
fclose (pff);

[A voir en vidéo sur Futura]

[julien_4230]

Ah en effet ! Merci bien !
J'aurais d'autres questions en C++... Par exemple :

J'ai trouvé sur gsl library de trouver les vecteurs propres et valeurs propres d'une matrice, mais je ne sais même pas vérifier la décomposition en vecteurs propres... Je me demande si c'est possible ? Voici une tentative de code pour faire cela :

#include <iostream>
#include <ctime>
#include <cstdlib>
#include <cmath>

#include <gsl/gsl_rng.h>
#include <gsl/gsl_math.h>
#include <gsl/gsl_randist.h>
#include <gsl/gsl_eigen.h>
#include <gsl/gsl_blas.h>

using namespace std ;

int
main (void)
{
double data[] = { 1.0 , 1/2.0, 1/3.0, 1/4.0,
1/2.0, 1/3.0, 1/4.0, 1/5.0,
1/3.0, 1/4.0, 1/5.0, 1/6.0,
1/4.0, 1/5.0, 1/6.0, 1/7.0 };

gsl_matrix_view m
= gsl_matrix_view_array (data, 4, 4);

gsl_vector *eval = gsl_vector_alloc (4);
gsl_matrix *evec = gsl_matrix_alloc (4, 4);

gsl_eigen_symmv_workspace * w =
gsl_eigen_symmv_alloc (4);

gsl_eigen_symmv (&m.matrix, eval, evec, w);

gsl_eigen_symmv_free (w);

gsl_eigen_symmv_sort (eval, evec,
GSL_EIGEN_SORT_ABS_ASC);

{
int i;

for (i = 0; i < 4; i++)
{
double eval_i = gsl_vector_get (eval, i);
gsl_vector_view evec_i = gsl_matrix_column (evec, i);

// cout << "eigenvalue = " << eval_i << "\n" ;
cout << evec_i.vector << "\n";
// gsl_vector_fprintf (stdout, &evec_i.vector, "%g");
}
}



// double data2[4][4];

// for (i=0; i<4 ; i++) {
// for (j=0 ; j<4 ; j++) {
// for (k=0 ; k<4 ; k++) {
// data2[i][j] += eval_i *

gsl_vector_free (eval);
gsl_matrix_free (evec);

return 0;
}


Bon sang, qu'est-ce qui ne va pas, damn it ?

Merci à vous.

inspiré de :
http://www.gnu.org/software/gsl/manu...-Examples.html