X-Git-Url: https://www.ginac.de/ginac.git//ginac.git?p=ginac.git;a=blobdiff_plain;f=ginac%2Ftensor.cpp;h=1a8f8e4041e542661e83b7190bfb2bd34fb052dd;hp=f8e5635a64af3caeb67c6557ca8e21cff6b50815;hb=54c1eca199508d821795ca08506f6e7788ebee91;hpb=5ef801553eb39aed7bd2df9dd1aff9d752c3ea9d diff --git a/ginac/tensor.cpp b/ginac/tensor.cpp index f8e5635a..1a8f8e40 100644 --- a/ginac/tensor.cpp +++ b/ginac/tensor.cpp @@ -3,7 +3,7 @@ * Implementation of GiNaC's special tensors. */ /* - * GiNaC Copyright (C) 1999-2003 Johannes Gutenberg University Mainz, Germany + * GiNaC Copyright (C) 1999-2005 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 @@ -17,7 +17,7 @@ * * You should have received a copy of the GNU General Public License * along with this program; if not, write to the Free Software - * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA + * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA */ #include @@ -33,18 +33,32 @@ #include "lst.h" #include "numeric.h" #include "matrix.h" -#include "print.h" #include "archive.h" #include "utils.h" namespace GiNaC { GINAC_IMPLEMENT_REGISTERED_CLASS(tensor, basic) -GINAC_IMPLEMENT_REGISTERED_CLASS(tensdelta, tensor) -GINAC_IMPLEMENT_REGISTERED_CLASS(tensmetric, tensor) -GINAC_IMPLEMENT_REGISTERED_CLASS(minkmetric, tensmetric) -GINAC_IMPLEMENT_REGISTERED_CLASS(spinmetric, tensmetric) -GINAC_IMPLEMENT_REGISTERED_CLASS(tensepsilon, tensor) + +GINAC_IMPLEMENT_REGISTERED_CLASS_OPT(tensdelta, tensor, + print_func(&tensdelta::do_print). + print_func(&tensdelta::do_print_latex)) + +GINAC_IMPLEMENT_REGISTERED_CLASS_OPT(tensmetric, tensor, + print_func(&tensmetric::do_print). + print_func(&tensmetric::do_print)) + +GINAC_IMPLEMENT_REGISTERED_CLASS_OPT(minkmetric, tensmetric, + print_func(&minkmetric::do_print). + print_func(&minkmetric::do_print_latex)) + +GINAC_IMPLEMENT_REGISTERED_CLASS_OPT(spinmetric, tensmetric, + print_func(&spinmetric::do_print). + print_func(&spinmetric::do_print_latex)) + +GINAC_IMPLEMENT_REGISTERED_CLASS_OPT(tensepsilon, tensor, + print_func(&tensepsilon::do_print). + print_func(&tensepsilon::do_print_latex)) ////////// // constructors @@ -174,7 +188,7 @@ ex tensdelta::eval_indexed(const basic & i) const exmap m; m[i1] = i1.replace_dim(min_dim); m[i2] = i2.replace_dim(min_dim); - return i.subs(m); + return i.subs(m, subs_options::no_pattern); } // Trace of delta tensor is the (effective) dimension of the space @@ -218,7 +232,7 @@ ex tensmetric::eval_indexed(const basic & i) const exmap m; m[i1] = i1.replace_dim(min_dim); m[i2] = i2.replace_dim(min_dim); - return i.subs(m); + return i.subs(m, subs_options::no_pattern); } // A metric tensor with one covariant and one contravariant index gets @@ -520,40 +534,51 @@ bool tensepsilon::contract_with(exvector::iterator self, exvector::iterator othe ex delta_tensor(const ex & i1, const ex & i2) { + static ex delta = (new tensdelta)->setflag(status_flags::dynallocated); + if (!is_a(i1) || !is_a(i2)) throw(std::invalid_argument("indices of delta tensor must be of type idx")); - return indexed(tensdelta(), sy_symm(), i1, i2); + return indexed(delta, symmetric2(), i1, i2); } ex metric_tensor(const ex & i1, const ex & i2) { + static ex metric = (new tensmetric)->setflag(status_flags::dynallocated); + if (!is_a(i1) || !is_a(i2)) throw(std::invalid_argument("indices of metric tensor must be of type varidx")); - return indexed(tensmetric(), sy_symm(), i1, i2); + return indexed(metric, symmetric2(), i1, i2); } ex lorentz_g(const ex & i1, const ex & i2, bool pos_sig) { + static ex metric_neg = (new minkmetric(false))->setflag(status_flags::dynallocated); + static ex metric_pos = (new minkmetric(true))->setflag(status_flags::dynallocated); + if (!is_a(i1) || !is_a(i2)) throw(std::invalid_argument("indices of metric tensor must be of type varidx")); - return indexed(minkmetric(pos_sig), sy_symm(), i1, i2); + return indexed(pos_sig ? metric_pos : metric_neg, symmetric2(), i1, i2); } ex spinor_metric(const ex & i1, const ex & i2) { + static ex metric = (new spinmetric)->setflag(status_flags::dynallocated); + if (!is_a(i1) || !is_a(i2)) throw(std::invalid_argument("indices of spinor metric must be of type spinidx")); if (!ex_to(i1).get_dim().is_equal(2) || !ex_to(i2).get_dim().is_equal(2)) throw(std::runtime_error("index dimension for spinor metric must be 2")); - return indexed(spinmetric(), sy_anti(), i1, i2); + return indexed(metric, antisymmetric2(), i1, i2); } ex epsilon_tensor(const ex & i1, const ex & i2) { + static ex epsilon = (new tensepsilon)->setflag(status_flags::dynallocated); + if (!is_a(i1) || !is_a(i2)) throw(std::invalid_argument("indices of epsilon tensor must be of type idx")); @@ -563,11 +588,13 @@ ex epsilon_tensor(const ex & i1, const ex & i2) if (!ex_to(i1).get_dim().is_equal(_ex2)) throw(std::runtime_error("index dimension of epsilon tensor must match number of indices")); - return indexed(tensepsilon(), sy_anti(), i1, i2); + return indexed(epsilon, antisymmetric2(), i1, i2); } ex epsilon_tensor(const ex & i1, const ex & i2, const ex & i3) { + static ex epsilon = (new tensepsilon)->setflag(status_flags::dynallocated); + if (!is_a(i1) || !is_a(i2) || !is_a(i3)) throw(std::invalid_argument("indices of epsilon tensor must be of type idx")); @@ -577,11 +604,14 @@ ex epsilon_tensor(const ex & i1, const ex & i2, const ex & i3) if (!ex_to(i1).get_dim().is_equal(_ex3)) throw(std::runtime_error("index dimension of epsilon tensor must match number of indices")); - return indexed(tensepsilon(), sy_anti(), i1, i2, i3); + return indexed(epsilon, antisymmetric3(), i1, i2, i3); } ex lorentz_eps(const ex & i1, const ex & i2, const ex & i3, const ex & i4, bool pos_sig) { + static ex epsilon_neg = (new tensepsilon(true, false))->setflag(status_flags::dynallocated); + static ex epsilon_pos = (new tensepsilon(true, true))->setflag(status_flags::dynallocated); + if (!is_a(i1) || !is_a(i2) || !is_a(i3) || !is_a(i4)) throw(std::invalid_argument("indices of Lorentz epsilon tensor must be of type varidx")); @@ -591,7 +621,7 @@ ex lorentz_eps(const ex & i1, const ex & i2, const ex & i3, const ex & i4, bool if (!ex_to(i1).get_dim().is_equal(_ex4)) throw(std::runtime_error("index dimension of epsilon tensor must match number of indices")); - return indexed(tensepsilon(true, pos_sig), sy_anti(), i1, i2, i3, i4); + return indexed(pos_sig ? epsilon_pos : epsilon_neg, antisymmetric4(), i1, i2, i3, i4); } } // namespace GiNaC