Accueil > > > CLASSE DE VECTEURS EN 2D / 3D PORTABLE, UTILISATION POUR SIMULATIONS PHYSIQUES
CLASSE DE VECTEURS EN 2D / 3D PORTABLE, UTILISATION POUR SIMULATIONS PHYSIQUES
Information sur la source
Description
J'ai utilisé cette classe de vecteurs dans plusieurs simulations physiques (cinétiques surtout, physique classique de Newton et Co). L'aspect graphique était géré en OpenGL, mais la classe en elle-même est absolument indépendante de TOUT. Je n'ai rien utilisé de spécial: je n'ai utilisé que iostream pour la fonction de débuggage console, si on l'enlève on peut aussi enlever l'utilisation de la STD. Et puis j'ai aussi inclus math.h, en écrivant à côté les fonctions que j'ai utilisées, donc si vous ne voulez pas utiliser math.h, il vous suffit de réécrire ces fonctions, et tout marchera! NOTE: le screenshot provient d'une simulation de collisions non-élastiques, les segments colorés sont les vecteurs vitesse des billes ATTENTION GRANDE MISE A JOUR DU 5 MAI 2005: -------------------------------------------- ------------------ je remplace le code (dégueulasse et informe) que j'avais posté ici (quelle daube :/) et je le remplace par deux nouvelles classes écrites récemment, et au passage, vous avez la version 2D et la version 3D. c'est doublement templaté, mais je vous suggères de garder les types proposés dans les 3 typedefs de bas de fichier: V3Di V3Df V3Dd et V2Di V2Df V2Dd. i = int, f = float, d = double. pas mal de fonctions ont été ajoutées, et puis surtout c'est plus propre ^^ Edit du 30 septembre: correction d'un bug pour GetAngle et ajout des "const" pour les méthodes concernées
Source
- //------------------------------------------------------------------------------
- #ifndef CVECTEUR2D_H
- #define CVECTEUR2D_H
- //écrit par Kirua - nicolasboumal@hotmail.com
- //------------------------------------------------------------------------------
- #include <cmath> //pour la fonction hypot() (renvoie la racine carrée de la somme des carrés des 2 paramètres)
- //------------------------------------------------------------------------------
- //Explication: le premier paramètre (T) est le type des coordonnées. Pour un point de l'écran, on voudra
- //des int, mais pour un point de l'espace, on préférera float ou même double. Le second paramètre (S)
- //indique le type à utiliser pour effectuer les calculs intermédiaires. Il est recommandé d'utilisé un
- //type plus précis (par défaut, double) afin d'éviter les pertes trop importantes liées aux arrondis.
- template<typename T, typename S = double>
- class CVecteur2D
- {
- public:
- CVecteur2D() : x(0), y(0) {}
- CVecteur2D(T _x, T _y) : x(_x), y(_y) {}
- CVecteur2D(const CVecteur2D& V) { x = V.x; y = V.y; };
-
- T x, y;
-
- S GetModule() const {return S(hypot(y,x));}
- void SetModule(S m) { S M = GetModule(); x = T(x*m/M); y = T(y*m/M); }
-
- S GetAngle() const { return S(atan2(y, x)); }
-
- CVecteur2D operator+ (const CVecteur2D& V) const {return CVecteur2D(x+V.x, y+V.y);}
- CVecteur2D operator- (const CVecteur2D& V) const {return CVecteur2D(x-V.x, y-V.y);}
- CVecteur2D operator* (const S R) const {return CVecteur2D(T(x*R), T(y*R));}
- CVecteur2D operator/ (const S R) const {return CVecteur2D(T(x/R), T(y/R));}
- S operator* (const CVecteur2D& V) const {return x*V.x + y*V.y;}
- CVecteur2D& operator+=(const CVecteur2D& V) {x += V.x; y += V.y; return *this;}
- CVecteur2D& operator-=(const CVecteur2D& V) {x -= V.x; y -= V.y; return *this;}
- CVecteur2D& operator*=(const S R) {x = T(S(x) * R); y = T(S(y) * R); return *this;}
- CVecteur2D& operator/=(const S R) {x = T(S(x) / R); y = T(S(y) / R); return *this;}
- bool operator==(const CVecteur2D& V) const {return (x == V.x && y == V.y);}
- bool operator!=(const CVecteur2D& V) const {return (x != V.x || y != V.y);}
- };
- //------------------------------------------------------------------------------
- typedef CVecteur2D<int> V2Di;
- typedef CVecteur2D<float> V2Df;
- typedef CVecteur2D<double> V2Dd;
- //------------------------------------------------------------------------------
- #endif
- //------------------------------------------------------------------------------
-
-
-
-
- Et la classe 3D:
-
-
-
- //------------------------------------------------------------------------------
- #ifndef CVECTEUR3D_H
- #define CVECTEUR3D_H
- //écrit par Kirua - nicolasboumal@hotmail.com
- //------------------------------------------------------------------------------
- #include <cmath>
- //------------------------------------------------------------------------------
- template<typename T, typename S = double>
- class CVecteur3D
- {
- public:
- CVecteur3D() : x(0), y(0), z(0) {}
- CVecteur3D(T _x, T _y, T _z) : x(_x), y(_y), z(_z) {}
- CVecteur3D(const CVecteur3D& V) { x = V.x; y = V.y; z = V.z; }
-
- T x, y, z;
-
- S GetModule() const { return S(sqrt(x*x + y*y + z*z)); }
- void SetModule(S m) { S M = GetModule(); x = T(x*m/M); y = T(y*m/M); z = T(z*m/M); }
-
- CVecteur3D operator+ (const CVecteur3D& V) const {return CVecteur3D(x+V.x, y+V.y, z+V.z);}
- CVecteur3D operator- (const CVecteur3D& V) const {return CVecteur3D(x-V.x, y-V.y, z-V.z);}
- CVecteur3D operator* (const S R) const {return CVecteur3D(T(x*R), T(y*R), T(z*R));}
- CVecteur3D operator/ (const S R) const {return CVecteur3D(T(x/R), T(y/R), T(z/R));}
- S operator* (const CVecteur3D& V) const {return x*V.x + y*V.y + z*V.z;}
- CVecteur3D& operator+=(const CVecteur3D& V) {x += V.x; y += V.y; z += V.z; return *this;}
- CVecteur3D& operator-=(const CVecteur3D& V) {x -= V.x; y -= V.y; z -= V.z; return *this;}
- CVecteur3D& operator*=(const S R) {x = T(S(x) * R); y = T(S(y) * R); z = T(S(z) * R); return *this;}
- CVecteur3D& operator/=(const S R) {x = T(S(x) / R); y = T(S(y) / R); z = T(S(z) / R); return *this;}
- bool operator==(const CVecteur3D& V) const {return (x == V.x && y == V.y && z == V.z);}
- bool operator!=(const CVecteur3D& V) const {return (x != V.x || y != V.y || z != V.z);}
-
- //produit croisé (cross product, produit vectoriel)
- CVecteur3D operator^ (const CVecteur3D& V) const {return CVecteur3D(y*V.z - z*V.y, z*V.x - x*V.z, x*V.y - y*V.x);}
- CVecteur3D& operator^= (const CVecteur3D& V) {x = y*V.z - z*V.y; y = z*V.x - x*V.z; z = x*V.y - y*V.x; return *this;}
- };
- //------------------------------------------------------------------------------
- typedef CVecteur3D<int> V3Di;
- typedef CVecteur3D<float> V3Df;
- typedef CVecteur3D<double> V3Dd;
- //------------------------------------------------------------------------------
- #endif
- //------------------------------------------------------------------------------
//------------------------------------------------------------------------------
#ifndef CVECTEUR2D_H
#define CVECTEUR2D_H
//écrit par Kirua - nicolasboumal@hotmail.com
//------------------------------------------------------------------------------
#include <cmath> //pour la fonction hypot() (renvoie la racine carrée de la somme des carrés des 2 paramètres)
//------------------------------------------------------------------------------
//Explication: le premier paramètre (T) est le type des coordonnées. Pour un point de l'écran, on voudra
//des int, mais pour un point de l'espace, on préférera float ou même double. Le second paramètre (S)
//indique le type à utiliser pour effectuer les calculs intermédiaires. Il est recommandé d'utilisé un
//type plus précis (par défaut, double) afin d'éviter les pertes trop importantes liées aux arrondis.
template<typename T, typename S = double>
class CVecteur2D
{
public:
CVecteur2D() : x(0), y(0) {}
CVecteur2D(T _x, T _y) : x(_x), y(_y) {}
CVecteur2D(const CVecteur2D& V) { x = V.x; y = V.y; };
T x, y;
S GetModule() const {return S(hypot(y,x));}
void SetModule(S m) { S M = GetModule(); x = T(x*m/M); y = T(y*m/M); }
S GetAngle() const { return S(atan2(y, x)); }
CVecteur2D operator+ (const CVecteur2D& V) const {return CVecteur2D(x+V.x, y+V.y);}
CVecteur2D operator- (const CVecteur2D& V) const {return CVecteur2D(x-V.x, y-V.y);}
CVecteur2D operator* (const S R) const {return CVecteur2D(T(x*R), T(y*R));}
CVecteur2D operator/ (const S R) const {return CVecteur2D(T(x/R), T(y/R));}
S operator* (const CVecteur2D& V) const {return x*V.x + y*V.y;}
CVecteur2D& operator+=(const CVecteur2D& V) {x += V.x; y += V.y; return *this;}
CVecteur2D& operator-=(const CVecteur2D& V) {x -= V.x; y -= V.y; return *this;}
CVecteur2D& operator*=(const S R) {x = T(S(x) * R); y = T(S(y) * R); return *this;}
CVecteur2D& operator/=(const S R) {x = T(S(x) / R); y = T(S(y) / R); return *this;}
bool operator==(const CVecteur2D& V) const {return (x == V.x && y == V.y);}
bool operator!=(const CVecteur2D& V) const {return (x != V.x || y != V.y);}
};
//------------------------------------------------------------------------------
typedef CVecteur2D<int> V2Di;
typedef CVecteur2D<float> V2Df;
typedef CVecteur2D<double> V2Dd;
//------------------------------------------------------------------------------
#endif
//------------------------------------------------------------------------------
Et la classe 3D:
//------------------------------------------------------------------------------
#ifndef CVECTEUR3D_H
#define CVECTEUR3D_H
//écrit par Kirua - nicolasboumal@hotmail.com
//------------------------------------------------------------------------------
#include <cmath>
//------------------------------------------------------------------------------
template<typename T, typename S = double>
class CVecteur3D
{
public:
CVecteur3D() : x(0), y(0), z(0) {}
CVecteur3D(T _x, T _y, T _z) : x(_x), y(_y), z(_z) {}
CVecteur3D(const CVecteur3D& V) { x = V.x; y = V.y; z = V.z; }
T x, y, z;
S GetModule() const { return S(sqrt(x*x + y*y + z*z)); }
void SetModule(S m) { S M = GetModule(); x = T(x*m/M); y = T(y*m/M); z = T(z*m/M); }
CVecteur3D operator+ (const CVecteur3D& V) const {return CVecteur3D(x+V.x, y+V.y, z+V.z);}
CVecteur3D operator- (const CVecteur3D& V) const {return CVecteur3D(x-V.x, y-V.y, z-V.z);}
CVecteur3D operator* (const S R) const {return CVecteur3D(T(x*R), T(y*R), T(z*R));}
CVecteur3D operator/ (const S R) const {return CVecteur3D(T(x/R), T(y/R), T(z/R));}
S operator* (const CVecteur3D& V) const {return x*V.x + y*V.y + z*V.z;}
CVecteur3D& operator+=(const CVecteur3D& V) {x += V.x; y += V.y; z += V.z; return *this;}
CVecteur3D& operator-=(const CVecteur3D& V) {x -= V.x; y -= V.y; z -= V.z; return *this;}
CVecteur3D& operator*=(const S R) {x = T(S(x) * R); y = T(S(y) * R); z = T(S(z) * R); return *this;}
CVecteur3D& operator/=(const S R) {x = T(S(x) / R); y = T(S(y) / R); z = T(S(z) / R); return *this;}
bool operator==(const CVecteur3D& V) const {return (x == V.x && y == V.y && z == V.z);}
bool operator!=(const CVecteur3D& V) const {return (x != V.x || y != V.y || z != V.z);}
//produit croisé (cross product, produit vectoriel)
CVecteur3D operator^ (const CVecteur3D& V) const {return CVecteur3D(y*V.z - z*V.y, z*V.x - x*V.z, x*V.y - y*V.x);}
CVecteur3D& operator^= (const CVecteur3D& V) {x = y*V.z - z*V.y; y = z*V.x - x*V.z; z = x*V.y - y*V.x; return *this;}
};
//------------------------------------------------------------------------------
typedef CVecteur3D<int> V3Di;
typedef CVecteur3D<float> V3Df;
typedef CVecteur3D<double> V3Dd;
//------------------------------------------------------------------------------
#endif
//------------------------------------------------------------------------------
Conclusion
Si vous trouvez des erreurs / incohérences ... ou si vous avez des suggestions: postez joyeusement.
Historique
- 05 mai 2005 17:58:39 :
- Remplacement total du code par deux nouvelles classes templatées de bien meilleure qualité (ça je peux le dire, maintenant, reste à savoir si le nouveau code vaut qq ch ^^).
- 05 mai 2005 18:02:03 :
- Remplacement du zip par le nouveau code.
- 30 septembre 2006 19:38:46 :
- Modification d'une ERREUR: la fonction GetAngle de la classe CVecteur2D appelait atan2 avec les mauvais arguments: c'est y puis x, et non l'inverse. Mettez vos codes à jour si vous avez utilisé cette classe ...
De plus, j'ai rajouté les "const" pour les méthodes utiles (+ - * / == != et ^ + GetAngle et GetModule).
- 30 septembre 2006 19:48:33 :
- Correction d'une bêtise dans le texte de présentation ^^.
Sources du même auteur
Sources de la même categorie
Commentaires et avis
Discussions en rapport avec ce code source dans le forum
Coord 2D to 3D [ par bat67000 ]
Comment optenir d'un point 2D sur l'app les coordonnees du point 3D associé avec la matrice de projection ?(je pige pas bien comment fonctionne la mat
de la 3d a la 2d [ par azkab ]
bonjour tt le monde,dans mon projet j ai un nombre de primitives (rectangles, cercles, polygones, lignes..) que j ai modelisé en 2d. je stock leur con
allegroGL ... [ par screamX ]
Je me tate en se moment pour développer un apps en 2D + 3D...je sais pas trop qu'elle lib utiliser... je 'maitrise' allegro mais la 3D c'est pas son f
Probleme vecteur 3D et projection orthogonale [ par supergrey ]
Salut, voila le probleme:J'ai deux vecteurs 3D, V(xv,yv,zv) et P(xp,yp,zp), et je voudrais calculer les vecteurs v0 et v1 de sorte que: v0+v1=V (
matrice*vecteur=vecteur? [ par Arnaud16022 ]
petit pb d'operator sur prografix on trouve: X' = A * X X' étant les nouvelles coordonnées du point, X les anciennes et A la matrice de t
Objet COM [ par 0xYg3n3 ]
Bonsoir, Voila je ne sais pas si je vais arriver a me faire comprendre. Je me demandais si un objet qui n'a pas initialiser d'interfaces COM peut qu
Comment fonctionne un moteur 2D ? [ par psykocrash ]
Bonjour,J'ai pour projet de développer un rpg à la Final Fantasy (en 2D, 4 héros, un inventaire pour les techniques et les armes, des c
une file Template dans un Vecteur Template [ par grand_duck ]
voici le problème j'ai une classe de vecteur template et une classe de file template je voudrais mettre ma file template dans mon vecteur est-ce
Image 2D sur fond 3D [ par CyberP ]
J'ai une question peut être stupide mais comment faire pour dessiner une image par dessus un rendu OpenGL ? Je pensais utiliser les trucs avec bi
[VC++] Affichage d'une en courbe 2D et 3D [ par crocgunter ]
Bonjour à tous..Je dispose d'une image BMP que je transforme en tableau de points.Dans un premier temps, je veux réaliser une coupe de l'ima
|
Derniers Blogs
IMAGINE CUP 2012, MAKE A SIGN EN FINALEIMAGINE CUP 2012, MAKE A SIGN EN FINALE par junarnoalg
Voilà qui est fait, la nouvelle est officielle ! L'équipe belge "Make a Sign" va au pays des kangourous défendre son projet dans la catégorie Software Design. http://www.imaginecup.com/CompetitionsContent/Competition/WorldwideFinalists.aspx V...
Cliquez pour lire la suite de l'article par junarnoalg KINECT 1.5 IS OUT !KINECT 1.5 IS OUT ! par Vko
La version 1.5 du Kinect For Microsoft vient tout juste de sortir ! Plein de nouveautés: Tracking de squelette en Near Mode Détection en position assise Détection faciale avec un SDK dédié Documentation et des guideline (enfin) Un out...
Cliquez pour lire la suite de l'article par Vko LES ACTUALITéS DE LA SEMAINE SUR C2I.FR (14 MAI - 20 MAI) LES ACTUALITéS DE LA SEMAINE SUR C2I.FR (14 MAI - 20 MAI) par richardc
Mise à jour des Web API du 14 Mai
Réservez dès maintenant votre journée du 20 juin pour le Windows Azure Dev Camp 2012 à Paris
Mise à jour de Team Foundation Service
MechCommander 2 sur Windows 8
Entity Framework 5 Release Candidate e...
Cliquez pour lire la suite de l'article par richardc REACTIVE EXTENSIONS : CONSOMMER DES SERVICES AVEC RX PARTIE 3, LES PIèGES à éVITERREACTIVE EXTENSIONS : CONSOMMER DES SERVICES AVEC RX PARTIE 3, LES PIèGES à éVITER par Groc
Une mauvaise utilisation de rx lors de l'écriture d'une couche d'accès à des services peut conduire à des cas embarassants avec des erreurs mal gérées, des appels qui ne partent lorsqu'ils le devraient, et même des résultats incorrects . le tout nuis...
Cliquez pour lire la suite de l'article par Groc SHAREPOINT BLOG SITE, PROBLèME D'ARCHIVESSHAREPOINT BLOG SITE, PROBLèME D'ARCHIVES par junarnoalg
Dernièrement, nous avons migré le site
myTIC
vers un nouveau serveur SharePoint 2010. Dans les contenus que nous vouloins récupérer, nous avions un certain nombre de blogs.
Nous avons utilisé les commandes Power...
Cliquez pour lire la suite de l'article par junarnoalg
Forum
MATRICE TEMPLATEMATRICE TEMPLATE par hjr2610
Cliquez pour lire la suite par hjr2610 RE : SAC A DOS RE : SAC A DOS par hadjkaddour
Cliquez pour lire la suite par hadjkaddour
Logiciels
sDEVIS-FACTURES vlPRO (8.1.0.3)SDEVIS-FACTURES VLPRO (8.1.0.3)sDEVIS-FACTURES vlPRO a été mis au point pour les particuliers, créateurs, entrepreneurs, artisa... Cliquez pour télécharger sDEVIS-FACTURES vlPRO 974 Application Server (12.2.4.6)974 APPLICATION SERVER (12.2.4.6)Développez de puissantes applications dans un environnement de 'cloud computing', clusterisé, séc... Cliquez pour télécharger 974 Application Server vPicture (1.4.2.1)VPICTURE (1.4.2.1)Avec vPicture, hébergez vos images facilement et rapidement.
vPicture est un utilitaire simple, ... Cliquez pour télécharger vPicture Easy-Planning (2.2.1.6)EASY-PLANNING (2.2.1.6)Easy-Planning permet de créer des plannings sous la représentation de diagrammes et est adapté au... Cliquez pour télécharger Easy-Planning COM-BACKUP (2.0)COM-BACKUP (2.0)
COM-BACKUP est un logiciel de sauvegarde qui permet de planifier les sauvegardes de vos dossiers ...
Cliquez pour télécharger COM-BACKUP
|