calcul de min et de max

calcul de min et de max - C++ - Programmation

Marsh Posté le 17-05-2007 à 03:08:58    

bonjour à tous,
 
J'ai un petit problème de faire le calcul de min et de max en c++.
En fait j'au un double He, tel que He est comprise entre 0.5<He<0.1e-020
mais on fait je vous dis ça mais moi je connais comment ça se pourra être puisque je calcul He en fonction de chaque maillage que j'ai donc je ne peux pas dire qu'elle est comprise entre telle et telle, puisque ça dépend de maillage est de son calcul. donc comment pourrai-je faire??????

Reply

Marsh Posté le 17-05-2007 à 03:08:58   

Reply

Marsh Posté le 17-05-2007 à 04:59:29    

on fait j'ai arrivé à trouver mon valeur de max mais j'arrive pas à trouver le bon min :(
 
voilà mon bou de code.

Code :
  1. #include <stdio.h>
  2. #include <stdlib.h> 
  3. #include <conio.h> 
  4. #include "vtkSurface.h"
  5. #include "RenderWindow.h"
  6. #include <vtkCommand.h>
  7. #include <vtkMath.h>
  8. #include <math.h>
  9. #include "vtkSurfaceBase.h"
  10. #include <fstream>
  12. vtkIdType Edge,T1,T2,i,j;
  13. double X[3],Y[3],Z[3],Vi[3],Vj[3],Normal[3],Normalize,CourbureMoyenne,Min=0,Max=0;
  15. int main( int argc, char* argv[] )
  17. {
  18. ofstream myfile;
  19. myfile.open ("Min et Max.txt" );
  20. vtkSurface *Surface=vtkSurface::New();
  21. Surface->CreateFromFile(argv[1]);---------------------------->c'est mon maillage en 3D
  22. int NbEdge=Surface->GetNumberOfEdges();---------------------------->ça me passe le nombre des arêtes
  24. for (int Edge=0;Edge<NbEdge;Edge++)
  25. {
  26. vtkIdType T1,T2,i,j;
  27. Surface->GetEdgeFaces(Edge,T1,T2);//elle me passe les 2 faces de chaque arête
  28. Surface->GetTriangleNormal(T1,Vi);//elle me passe le normal
  29. Surface->GetTriangleNormal(T2,Vj);
  30. Normal[0]=Vi[0]+Vj[0];
  31. Normal[1]=Vi[1]+Vj[1];
  32. Normal[2]=Vi[2]+Vj[2];//additionner les 2 vecteurs Vi et Vj
  33. vtkMath::Normalize(Normal);
  34. Surface->GetEdgeVertices(Edge,i,j);
  35. double Alpha[3];
  36. vtkMath::Cross(Normal,Vi,Alpha);
  37. double SinAlpha=vtkMath::Norm(Alpha);
  38. CourbureMoyenne=2*(Surface->GetDistanceBetweenVertices(i,j))*SinAlpha;
  39. if (CourbureMoyenne!=0)
  40. {
  41.  if(Min<=pow(10,-30) || Max>=1)
  42.           {
  43.      Min=CourbureMoyenne;
  44.             Max=CourbureMoyenne;
  45.   myfile<<"min="<<Min<<"  "<<"max="<<Max<<endl;
  47.            }
  48.         if(CourbureMoyenne<Min) 
  49.             Min=CourbureMoyenne; 
  50.         if(CourbureMoyenne>Max) 
  51.             Max=CourbureMoyenne;
  52.     }
  54. }
  55. cout<<"min="<<Min<<"  "<<"max="<<Max<<endl;
  56. myfile.close();
  57. return(0); 
  58. }

 
Merci pour vos aides :jap:

Reply

Sujets relatifs:

Leave a Replay

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