[GiNaC-list] Problem with scalar_products : Scalar product of two vectors in index notation

[aravind b]

Hi,
Hope everyone's doing fine !

Please go through this code :


#include <iostream>
#include <ginac/ginac.h>
using namespace std;
using namespace GiNaC;

int main()
{
    symbol i_sym("i"), j_sym("j"), e("e"), a("a"), b("b") ;
    idx i(i_sym, 3), j( j_sym, 3) ;          //defining indices i and j

    scalar_products sp;
    sp.add( indexed(e, i), indexed(e, j), delta_tensor(i,j) );       // e.i*
e.j = delta_tensor.i.j

    cout << (indexed(e,i)*indexed(e,j)).simplify_indexed(sp) <<endl ;

    return 0;
}

does NOT give delta.j.i instead just gives e.i * e.j - why?

Basically, i want to define an orthonormal basis vector set {e.i}.
If i have vector_a = indexed(a, i) * indexed(e, i) and vector_b = indexed(b,
j) * indexed(e, j) , then i want adotb to reduce to a.i * b.i ( ofcoure
after i apply simplify_indexed(sp) to adotb )

Thanks a lot for going through this mail.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://www.cebix.net/pipermail/ginac-list/attachments/20061109/b47dc4d8/attachment.html