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.