Discussion :

[bisontin]

Bonjour,
voila j'ai fait une petite modif ds un test d'un code de calcul, ce test utlise une viscosité _mu et le but de ma modification est d'utiliser deux viscosités _mu1 et _mu2;
à la ligne 288 de NavierStokesSolverPC.hpp le code utlise la fonction stiff qui va la chercher dans elemOper.hpp (ligne 953) et elemOper.cpp (ligne 60), dans ces deux fichiers j'ai modifié les arguments de la fonction stiff ainsi dans NavierStokesSolverPC.hpp. quand je lance le test j'ai une erreur:

Code :

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

1
2
3
4
main.o(.gnu.linkonce.t._ZN5LifeV20NavierStokesSolverPCINS_12RegionMesh3DINS_11LinearTetraENS_12MarkerCommonINS_17MarkerTraits_BaseEEEEEEC1ERK6GetPotRKNS_5RefFEESD_RKNS_8QuadRuleESG_SG_SG_RNS_9BCHandlerERKi+0xcb1): In function `LifeV::NavierStokesSolverPC<LifeV::RegionMesh3D<LifeV::LinearTetra, LifeV::MarkerCommon<LifeV::MarkerTraits_Base> > >::NavierStokesSolverPC(GetPot const&, LifeV::RefFE const&, LifeV::RefFE const&, LifeV::QuadRule const&, LifeV::QuadRule const&, LifeV::QuadRule const&, LifeV::QuadRule const&, LifeV::BCHandler&, int const&)':
: undefined reference to `LifeV::stiff(double, double, LifeV::ElemMat&, LifeV::CurrentFE const&, int, int, int)'
collect2: ld returned 1 exit status
make[1]: *** [test_ns_cyl] Erreur 1

elemOper.cpp : http://cjoint.com/?kznXnogUDJ
Avez vous une ideé sur la nature de l'erreur

[zooro]

Bonjour,

Oui, aucun prototype compatible de la fonction LifeV::stiff n'a été trouvé.
Soit la fonction n'existe pas ou n'est pas visible de l'endroit où elle est appelée, soit les arguments utilisés à l'appel ne correspondent pas aux arguments des prototypes déclarés pour la fonction.

PS: en utilisant la balise CODE lors de la rédaction de ton post, tu aurais simplifié sa lecture. Peut-être que plus de personnes auraient pris le temps de lire ton message

[bisontin]

voila j'ai pas vu la balise code
sinon je suppose que c'est les arguments utilisés à l'appel qui ne correspondent pas aux arguments déclarés pour stiff. peux tu m'eclairer un peu plus sur ce point parceque jene vois pas du tout parceque a l'appel de stiff il y a bien deux coefficients mu1 et mu2 et idem pour ceux declarés dans stiff.

[zooro]

Définition de la fonction :

Code :

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

1
2
void stiff( Real coef1, Real coef2, ElemMat& elmat, const CurrentFE& fe,
            int iblock, int jblock, int nb )

Appel de la fonction avec :

Code :

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

1
2
LifeV::stiff(double, double, LifeV::ElemMat&, LifeV::CurrentFE const&,
            int, int, int)

Donc déjà, il veut des "Real" et tu lui donnes des "double".

[bisontin]

Oui c'est ce que je me suis dit au début ca saute aux yuex mais avant ma modif,
c'était à l'appel :

Code :

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

1
2
void stiff( Real coef, ElemMat& elmat, const CurrentFE& fe,
 int iblock, int jblock, int nb )

et ca marchait parfaitement

[zooro]

Oui c'est ce que je me suis dit au début ca saute aux yuex mais avant ma modif,
c'était à l'appel :

Code :

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

1
2
void stiff( Real coef, ElemMat& elmat, const CurrentFE& fe,
 int iblock, int jblock, int nb )

et ca marchait parfaitement

Bizarre.
Deux solutions simples :
- soit tu forces le cast des double en Real (si ça marchait avant ta modif, ça doit être possible);
- soit tu ajoutes un prototype avec des doubles à la place des Real.