un problème dans une équation c++

un problème dans une équation c++ - C++ - Programmation

Marsh Posté le 18-01-2007 à 14:39:06    

j'ai en fait un problème lors de debuger mon programme et c'est dans cette fonction ou il se trouve le problème, mais quand je fais un Build pour mon programme il y 0 erreur, donc si quelqu'un pourra m'aider à trouver ce problème dans cet méthode:

Code :
  1. void AccumulateLaplacian()//X+Lambdadt*Laplacian
  2.  {
  4.   double X[3],Lambdadt=10;
  5.   this->AllocateTable();//ici je fais un appel à la méthode qui me fait l'allocation des tableaux TablePoints,TableLaplacian,Accumulate
  7. for (int i=0; i<this->Input->GetNumberOfPoints(); i++)
  8.   {
  10.    X[0]=TablePoints[i][0];
  11.    X[1]=TablePoints[i][1];
  12.    X[2]=TablePoints[i][2];
  13.    this->ComputeLaplacian();//je fais un appel à la méthode qui me fais le calcul de laplacienne
  15.   }
  16.    Accumulate[i][0]=X[0]+Lambdadt*TableLaplacian[i][0];//et c'est ici ou il se trouve mon problème  
  17.    Accumulate[i][1]=X[1]+Lambdadt*TableLaplacian[i][1];//c'est une équation qui fait X[0]+Lambdadt*L(X)
  18.    Accumulate[i][2]=X[2]+Lambdadt*TableLaplacian[i][2];//et je mets le résultat dans un tableau
  19.  }


merci en avance pour vos réponses

Reply

Marsh Posté le 18-01-2007 à 14:39:06   

Reply

Marsh Posté le 18-01-2007 à 15:01:28    

y a plein de trucs bizarres dans ton code :

  • ligne 13: comment ComputeLaplacian() est-elle au courant de la valeur de X. Si ComputeLaplacian() n'utilise pas X, pourquoi l'appeler à l'intérieur de la boucle ?
  • ligne 16-18: i n'est plus défini à cet endroit là


je ne comprends pas grand chose à ce que tu veux faire, mais il me semble qu'un truc dans ce genre là aurait plus de sens :

Code :
  1. this->AllocateTable();
  2. this->ComputeLaplacian(); // gros point d'interrogation ici: quelles sont les données utilisée pour le calcul du laplacien ?
  4. for( int i=0 ; ... )
  5. {
  6.   Accumulate[i][0] = Points[i][0] + Lambda*Laplacian[i][0];
  7.   Accumulate[i][1] = Points[i][1] + Lambda*Laplacian[i][1];
  8.   Accumulate[i][2] = Points[i][2] + Lambda*Laplacian[i][2];
  9. }


 
essaie de nous expliquer plus clairement ce que tu veux faire, et quelles sont les données que tu manipules...


---------------
TriScale innov
Reply

Marsh Posté le 18-01-2007 à 15:22:19    

voilà je te mets tout le programme c'est mieux  
 

Code :
  1. class vtkImplicitFairing
  2. {
  3. public:
  5. void SetInput(vtkSurface *Input) {this->Input=Input;};
  7. double X[3];
  8. double** TableLaplacian;
  9. double** TablePoints;
  10. double** Accumulate;
  12. void AllocateTable()
  13. {
  15. TablePoints = new double* [this->Input->GetNumberOfPoints()];
  16. Accumulate=new double* [this->Input->GetNumberOfPoints()];
  17. TableLaplacian = new double* [this->Input->GetNumberOfPoints()];
  19. for (int i=0;i<this->Input->GetNumberOfPoints();i++)
  20. {
  21.     TablePoints[i]=new double[3] ;
  22. TableLaplacian[i]=new double[3] ;
  23. Accumulate[i]=new double[3] ;
  24. }
  27. }
  29. void TableCoordinates()
  30. {
  31.  //déclaration + allocation d'un tableau de 2 dimensaions et de i chaines
  33.  for (int i=0;i<this->Input->GetNumberOfPoints();i++ )
  34.  {
  35.    this->Input->GetPointCoordinates(i,X);//Xi
  36.    TablePoints[i][0]=X[0];
  37.    TablePoints[i][1]=X[1];//j'ai mit une copie des mes sommets dans un tableau
  38.    TablePoints[i][2]=X[2];
  40.  }
  42. }
  43. void ComputeLaplacian()
  44.  {
  45.   int m=6;
  46.   double** TableLaplacian = new double* [this->Input->GetNumberOfPoints()] ;
  47.   double Y[3],A[3],L[3];
  48.   vtkIdList *MyList=vtkIdList::New();
  49.   for (int i=0; i<this->Input->GetNumberOfPoints(); i++)
  50.   {
  52.    this->Input->GetVertexNeighbours(i,MyList);//elle me passe les sommets qui ont des voisinages avec le sommet choisi
  54.    TablePoints[i][0]=X[0];
  55.    TablePoints[i][1]=X[1];
  56.    TablePoints[i][2]=X[2];
  58.    double somme[3];
  59.    somme[0]=0; somme[1]=0;somme[2]=0;
  62.   for (int j=0;j<MyList->GetNumberOfIds();j++) //le nbre d'ID correspond au nombre des points de voisinage
  63.   {
  64.                 this->Input->GetPointCoordinates(MyList->GetId(j),Y);//prendre les voisins Y(j) de i
  65.    TablePoints[j][0]=Y[0];
  66.    TablePoints[j][1]=Y[1];
  67.    TablePoints[j][2]=Y[2];
  69.    A[0]=TablePoints[j][0]-TablePoints[i][0];
  70.    A[1]=TablePoints[j][1]-TablePoints[i][1];//on soustraire Y(j) de X(i) pour chaque axe
  71.    A[2]=TablePoints[j][2]-TablePoints[i][2];
  73.    somme[0]=somme[0]+A[0];
  74.    somme[1]=somme[1]+A[1];//pour le calcul de Laplacian Discret= somme (y(j)-X(i))  
  75.    somme[2]=somme[2]+A[2];
  77.    L[0]=1/m*somme[0];
  78.    L[1]=1/m*somme[1];//Laplacian d'ordre 1 L(X)=1/6*somme(Y(j)-X(i))
  79.    L[2]=1/m*somme[2];
  80.   }
  81.   TableLaplacian[i][0]=L[0];
  82.   TableLaplacian[i][1]=L[1];//j'ai mit une copie des mes sommets dans un tableau
  83.   TableLaplacian[i][2]=L[2];
  84.  }
  86.  MyList->Delete();
  87. }
  90.  void AccumulateLaplacian()//X+Lambdadt*Laplacian
  91.  {
  93.   double X[3],Lambdadt=10;
  94.   for (int i=0; i<this->Input->GetNumberOfPoints(); i++)
  95.   {
  97.    X[0]=TablePoints[i][0];
  98.    X[1]=TablePoints[i][1];
  99.    X[2]=TablePoints[i][2];
  100.    this->ComputeLaplacian();
  102.   }
  103.    Accumulate[i][0]=X[0]+Lambdadt*TableLaplacian[i][0];
  104.    Accumulate[i][1]=X[1]+Lambdadt*TableLaplacian[i][1];
  105.    Accumulate[i][2]=X[2]+Lambdadt*TableLaplacian[i][2];
  106.  }
  108.  void ApplyImplicitFilter(double Lambdadt, int n)
  109.  {
  110.   this->AllocateTable();
  111.   for (int n=0;n<4;n++)
  112.   {
  113.    double Lambdadt=0, nombre=10;
  114.    Lambdadt=pow(nombre,n);//lambdadt=10^n;
  115.    this->ComputeLaplacian();
  116.    this->AccumulateLaplacian();
  118.   }
  119.   for (int i=0;i<this->Input->GetNumberOfPoints();i++)
  120.    {
  121.     this->Input->SetPointCoordinates(i, Accumulate[i]);//je fais réagir les nouveaux positions de mes sommets  
  122.    }
  123.  }
  124.  int vtkImplicitFairing::free(double **TablePoints,double **Accumulate,double **TableLaplacian);
  127. protected:
  129. vtkSurface *Input;
  131. };
  135. int main( int argc, char* argv[] )//argc contient le nombre de parametre passé au programme +1  
  136. //argv est un tableau contenant les chaines de caractere transmises
  137. {
  139. int i;
  140. for(i=0;i<argc;i++)
  141. {
  142.  printf("Argument %1i : %s \n", i+1, argv[i]);
  143. }
  146. //Lire le fichier
  148. vtkSurface *Create=vtkSurface::New();
  149. Create->CreateFromFile(argv[1]);
  151. RenderWindow *Window=RenderWindow::New();
  152. Window->SetInput(Create);
  154. //Filtrage
  155. vtkImplicitFairing *Fairing=new (vtkImplicitFairing);
  156. Fairing->SetInput(Create);
  157. Fairing->ApplyImplicitFilter(atof(argv[2]),atoi(argv[3]));
  159. Window->Render();
  160. Window->Interact();
  162. //sauvgarde de maillage
  163. vtkPolyDataWriter *Data=vtkPolyDataWriter::New();
  164. Data->SetInput(Create);
  165. Data->SetFileName("Implicit.vtk" );
  166. Data->Write();
  167. return(0);
  169. }


j'espére que tu vas comprendre ce que je voudrai

Reply

Marsh Posté le 18-01-2007 à 15:23:19    

Reply

Marsh Posté le 18-01-2007 à 15:27:09    

je n'ai pas compris tout à fais ce que tu voudras me dire j'ai un problème avec la compilation et ce problème vient d'une méthode et j'ai mis tout le programme pour que la personne comprend le déroulement du programme voilà. Si tu pourras m'aider en expliquant ou il se trouve le problème je serai ravis merci. :D

Reply

Marsh Posté le 18-01-2007 à 15:33:36    

Quand francesco te dit que "i" n'a pas la bonne portée ca te parle ?


---------------
Töp of the plöp
Reply

Marsh Posté le 18-01-2007 à 15:48:39    

c'est bon je l'ai trouvé c'est exactement ce qu'il m'a dit je n'ai pas fais gaffe mais je toujours un problème dans la premiére argument il arrive pas à lire mon fichier pourquoi??

Reply

Marsh Posté le 18-01-2007 à 15:53:57    

:??:  
 
tu peux parler en français (accords, ponctuation, ...) stp ?  Parce que là, j'ai vraiment du mal à comprendre...


---------------
TriScale innov
Reply

Marsh Posté le 18-01-2007 à 15:59:51    

ok , on Fait j'arrive à l'exécuter mais, j'ai toujours un problème dans l'argument [1], parce que le programme n'arrive pas à le lire pourtant, je l'ai mis en argument (propriétés->debugging->commandArgument)ici je mets mes arguments . mais il ne le lis pas . Pourqoui????
j'espére avoir être clair cette fois.

Reply

Marsh Posté le 18-01-2007 à 16:31:53    

qu'est-ce que les lignes 139-143 t'affichent à l'écran ?
 
Je ne sais pas comment marche ton environnement de développement, mais tu peux aussi tester en mettant tes arguments directement en ligne de commande dans une console:

myprog.exe arg1 arg2


 
(en supposant que tu es sous Windows)


---------------
TriScale innov
Reply

Marsh Posté le 18-01-2007 à 16:31:53   

Reply

Marsh Posté le 18-01-2007 à 16:45:16    

Code :
  1. double** TableLaplacian;
  2. double** TablePoints;
  3. double** Accumulate;
  5. void AllocateTable()
  6. {
  7.    TablePoints = new double* [this->Input->GetNumberOfPoints()];
  8.    Accumulate=new double* [this->Input->GetNumberOfPoints()];       
  9.    TableLaplacian = new double* [this->Input->GetNumberOfPoints()];
  11.    for (int i=0;i<this->Input->GetNumberOfPoints();i++)
  12.    {
  13.       TablePoints[i]=new double[3] ;
  14.       TableLaplacian[i]=new double[3] ;
  15.       Accumulate[i]=new double[3] ;
  16.    }


 
wééééééééééééééé :D
 
bon plus sérieusement: std::vector<>
 

Reply

Sujets relatifs:

Leave a Replay

Make sure you enter the(*)required information where indicate.HTML code is not allowed