Discussion :

[moooona]

Bonjour,
J'ai besoin d'utiliser la fonction suivante

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
 
extern "C" {
int dsyev_(char *jobz, char *uplo, 
	   int *n, double *a, int *lda, 
	   double *w, double *work, int *lwork, 
	   int *info);
}
 
void softassign::eigenvectorOfN(double *N, float* q){
 
  static float q_pre[4]; // previous result
 
  int dimN = 4;
  double w[4]; // eigenvalues
  double *work = new double; // workspace
  int info;
  int lwork = -1;
 
  dsyev_((char*)"V", (char*)"U",
	 &dimN, N, &dimN,
	 w, work, &lwork, &info);
  if(info != 0){
    fprintf(stderr, "info = %d\n", info);
    exit(1);
  }
  lwork = (int)work[0];
  delete work;
 
  work = new double [lwork];
 
  dsyev_((char*)"V", (char*)"U",
	 &dimN, N, &dimN,
	 w, work, &lwork, &info);
 
  delete [] work;
 
 
  if(info != 0){
    fprintf(stderr, "computing eigenvector FAIL! info = %d\n", info);
    //exit(1);
 
    // if fail, put back the previous result
    for(int i=0; i<4; i++){
      q[i] = q_pre[i];
    }
 
 
  }else{
 
    // last column of N is the eigenvector of the largest eigenvalue 
    // and N is stored column-major
    for(int i=0; i<4; i++){
      q[i] = N[4*3 + i];
      q_pre[i] = q[i];
    }
 
  }
 
 
}

Cette fonction m'a crée beaucoup de problèmes à cause de l'utilisation de dsyev_ de la bibliothèque Lapack.
Est ce qu'il y a quelqu'un qui peut me donner l'implémentation en c de cette fonction
Merci beaucoup pour votre aide

[bacelar]

1- C'est quoi le rapport avec C++/CLI ?
2- Google est ton ami : http://www.netlib.org/clapack/old/double/dsyev.c
3- votre code est extrêmement alambiqué et je pense même faux dû à de multiple cast sauvages à la C, des new et des delete complètement inutiles, etc...

[moooona]

Merci bacelar pour votre réponse..

Concernant le premier point : ce code crée un problème dans un projet c++/CLI. EN fait en ajoutant ce code , l'application ne s'exécute pas.
Concernant le deuxième point: D'aprés google j'ai trouvé que la fonction dsyev retourne eigenvalues et eigenvector mais d'aprés ce que j'ai implémenté l'eigenvector est un vecteur avec nombre réel et pas des nombres complexes alors que cette fonction retourne des valeurs complexes
Concernant troisième point: J'ai pas compris qu'est ce que vous voulez dire exactement mais ce que je cherche c'est l'implémentation de dsyev.

[bacelar]

>Concernant le premier point : ce code crée un problème dans un projet c++/CLI. EN fait en ajoutant ce code , l'application ne s'exécute pas.
Oui, peut-être, mais je pense que si le projet serait en C++ standard vous auriez les mêmes problèmes, car votre code utilisateur est une horreur en terme de complexité, de fuite mémoire, de maintenabilité.
En C++ moderne, il n'y quasiment pas de cast, et jamais de cast à la C, quasiment pas de delete et encore moins de delete [].
Franchement, plutôt que de jeter l'opprobre sur l'implémentation de dsyev_ dans une librairie extrêmement utilisées donc extrêmement bien testées, je vérifierais mon code.
Votre code ne devrait pas passer la moindre revue de code de vos paires.
Il devrait être plus proche du code qui suit, qui est thread-safe, exception-safe, plus facile à maintenir et à débuguer car utilisant les exceptions (la création du classe d'exception dédiée serait un plus), initialise automatiquement les data données à une librairie avec des valeurs par défaut (les routines C "systèmes" sont extrêmement susceptible sur l'initialisation de TOUTES les données, même les champs non (encore) utilisés).

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
const int DIM = 4;
int MatrixOder=12;
char V = 'V';
char U = 'U';
 
void softassign::eigenvectorOfN(double *N, float* q)
{
	std::array<double,DIM> w = {};
 
	int info = 0;
	double optimalSize;
	int lwork = -1;
	std::ostringstream msg;
 
	dsyev_(&V, &U,&MatrixOder, N, &MatrixOder,&w[0], &optimalSize, &lwork, &info);
 
	if(info != 0) { msg << "info =" << info << std::endl; throw (msg.str().c_str()); }
 
	std::vector<double> work(optimalSize); // workspace
	lwork = work.size();
 
	dsyev_(&V, &U,&MatrixOder, N, &MatrixOder, &w[0], &work[0], &lwork, &info);
 
	if(info != 0) { msg << "computing eigenvector FAIL! info =" << info << std::endl; throw (msg.str().c_str());}
}

[Médinoc]

Attention, à moins que j'aie encore manqué un truc sur C++11, les additions de chaînes littérales ne marcheront pas en C++ standard (et en C++/CLI, elles devraient être des chaînes larges).

De plus, si les deux premiers paramètres de dsyev_() sont censés être des chaînes C, le passage de &V et &U se base sur des détails d'implémentation sur l'alignement des variables. Ah non j'ai compris, c'est une lib FORTRAN, donc tous les paramètres sont passés par référence.

Edit: Franchement, la première chose que je ferais pour cet appel, c'est faire un wrapper de dsyev() pour "dé-pointeuriser" les arguments non-tableaux déclarés en entrée seule. Et sans doute constifier les tableaux en entrée seule, au passage.
PS: Le lien posté par bacelar est désormais 404. Version du cache Google

[bacelar]

>Attention, à moins que j'aie encore manqué un truc sur C++11
C'est mon poil dans la main, que ta manqué ;-).