Discussion :

[erloroloco]

Salut !

Si je calcule le cross product de deux vectors comme A=[1 3 2] et B=[5 1 2] a main, mon résultat est C[4 -8 -14]. Mais si j'utilise la formule en Fortran
C=cshift(A,1)*cshift(B,-1)-cshift(B,1)cshift(A,-1) le résultat est différent : C=[4 8 -14]. Le deux vectors résultats ne sont pas linear dependant.

Est-ce que vous savez si la formule est une autre, ou utilisez-vous un autre méthode ??

Merci d'avance...

[FR119492]

Salut !

Je ne comprends pas très bien ta formule, mais pourquoi ne pas faire simplement:

Code :

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

1
2
3
      C(1)=A(2)*B(3)-A(3)*B(2)
      C(2)=A(3)*B(1)-A(1)*B(3)
      C(3)=A(1)*B(2)-A(2)*B(1)

Jean-Marc Blanc

[erloroloco]

Ça marche naturalment, mais je ai cherché un formul general au lieu de 3. La formul montrée appartiens le repertoire de "La boˆıte `a outils de la PHYSIQUE NUM´ERIQUE".

MerÇi pour l´essai

[FR119492]

Salut !

Je lis dans l'article "Produit vectoriel" de Wikipedia:

Le produit vectoriel est une opération vectorielle effectuée dans les espaces euclidiens orientés de dimension trois.

Jean-Marc Blanc