X-Git-Url: https://www.ginac.de/ginac.git//ginac.git?a=blobdiff_plain;f=check%2Fexam_clifford.cpp;h=6262af8ae07fb63b8cfc873f3f3064f476c0c69f;hb=67467d256b44f5e08498ca81c946d9ffaa25d1e2;hp=51e34f899878bbf0e7ee077c1746697561195609;hpb=619d77d2676f7f1a562fb9fefc0ba6754fe2d750;p=ginac.git diff --git a/check/exam_clifford.cpp b/check/exam_clifford.cpp index 51e34f89..6262af8a 100644 --- a/check/exam_clifford.cpp +++ b/check/exam_clifford.cpp @@ -3,7 +3,7 @@ * Here we test GiNaC's Clifford algebra objects. */ /* - * GiNaC Copyright (C) 1999-2007 Johannes Gutenberg University Mainz, Germany + * GiNaC Copyright (C) 1999-2008 Johannes Gutenberg University Mainz, Germany * * This program is free software; you can redistribute it and/or modify * it under the terms of the GNU General Public License as published by @@ -20,7 +20,10 @@ * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA */ -#include "exams.h" +#include +#include "ginac.h" +using namespace std; +using namespace GiNaC; const numeric half(1, 2); @@ -410,13 +413,23 @@ template unsigned clifford_check6(const matrix &A) result += check_equal(canonicalize_clifford(e), 0); /* lst_to_clifford() and clifford_inverse() check*/ - realsymbol x("x"), y("y"), t("t"), z("z"); + realsymbol s("s"), t("t"), x("x"), y("y"), z("z"); ex c = clifford_unit(nu, A, 1); e = lst_to_clifford(lst(t, x, y, z), mu, A, 1) * lst_to_clifford(lst(1, 2, 3, 4), c); e1 = clifford_inverse(e); result += check_equal_simplify_term2((e*e1).simplify_indexed(), dirac_ONE(1)); +/* lst_to_clifford() and clifford_to_lst() check for vectors*/ + e = lst(t, x, y, z); + result += check_equal_lst(clifford_to_lst(lst_to_clifford(e, c), c, false), e); + result += check_equal_lst(clifford_to_lst(lst_to_clifford(e, c), c, true), e); + +/* lst_to_clifford() and clifford_to_lst() check for pseudovectors*/ + e = lst(s, t, x, y, z); + result += check_equal_lst(clifford_to_lst(lst_to_clifford(e, c), c, false), e); + result += check_equal_lst(clifford_to_lst(lst_to_clifford(e, c), c, true), e); + /* Moebius map (both forms) checks for symmetric metrics only */ matrix M1(2, 2), M2(2, 2); c = clifford_unit(nu, A); @@ -529,7 +542,6 @@ unsigned exam_clifford() unsigned result = 0; cout << "examining clifford objects" << flush; - clog << "----------clifford objects:" << endl; result += clifford_check1(); cout << '.' << flush; result += clifford_check2(); cout << '.' << flush; @@ -590,12 +602,10 @@ unsigned exam_clifford() result += clifford_check7(indexed(-2*minkmetric(), sy_symm(), xi, chi), dim); cout << '.' << flush; result += clifford_check7(-2*delta_tensor(xi, chi), dim); cout << '.' << flush; - if (!result) { - cout << " passed " << endl; - clog << "(no output)" << endl; - } else { - cout << " failed " << endl; - } - return result; } + +int main(int argc, char** argv) +{ + return exam_clifford(); +}