martes, 4 de julio de 2023

Ejemplo de una operación (sencilla) de contracción tensorial de índices empleando la herramienta GNU Maxima

Se trata de reproducir la siguiente operación de contracción de índices $b^{i\,\ell}\,c_{ij}\,a^{jk}$, que, cláramente, da como resultado $c^{\ell\,k}$. Lo que sigue es el código que he utilizado para pedirle a MAXIMA que reproduzca automáticamente este sencillo cálculo:

  (%i1)	load("itensor")$ /*carga de la librería necesaria */

  (%i2)	imetric(b); 
  (%o2)	done
  
  (%i3)	ishow(contract(a([-j,-k],[])*b([-i,-l],[])*c([i,j],[])))$
  (%t3)	c^l k

$\diamond$

No hay comentarios:

Publicar un comentario