* Implementation of GiNaC's special tensors. */
/*
- * GiNaC Copyright (C) 1999-2007 Johannes Gutenberg University Mainz, Germany
+ * GiNaC Copyright (C) 1999-2019 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
* Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA
*/
-#include <iostream>
-#include <stdexcept>
-#include <vector>
-
#include "tensor.h"
#include "idx.h"
#include "indexed.h"
#include "archive.h"
#include "utils.h"
+#include <iostream>
+#include <stdexcept>
+#include <vector>
+
namespace GiNaC {
GINAC_IMPLEMENT_REGISTERED_CLASS(tensor, basic)
// constructors
//////////
-tensor::tensor() : inherited(&tensor::tinfo_static)
+tensor::tensor()
{
setflag(status_flags::evaluated | status_flags::expanded);
}
minkmetric::minkmetric() : pos_sig(false)
{
- tinfo_key = &minkmetric::tinfo_static;
}
spinmetric::spinmetric()
{
- tinfo_key = &spinmetric::tinfo_static;
}
minkmetric::minkmetric(bool ps) : pos_sig(ps)
{
- tinfo_key = &minkmetric::tinfo_static;
}
tensepsilon::tensepsilon() : minkowski(false), pos_sig(false)
{
- tinfo_key = &tensepsilon::tinfo_static;
}
tensepsilon::tensepsilon(bool mink, bool ps) : minkowski(mink), pos_sig(ps)
{
- tinfo_key = &tensepsilon::tinfo_static;
}
//////////
// archiving
//////////
-DEFAULT_ARCHIVING(tensor)
-DEFAULT_ARCHIVING(tensdelta)
-DEFAULT_ARCHIVING(tensmetric)
-DEFAULT_ARCHIVING(spinmetric)
-DEFAULT_UNARCHIVE(minkmetric)
-DEFAULT_UNARCHIVE(tensepsilon)
-
-minkmetric::minkmetric(const archive_node &n, lst &sym_lst) : inherited(n, sym_lst)
+void minkmetric::read_archive(const archive_node& n, lst& sym_lst)
{
+ inherited::read_archive(n, sym_lst);
n.find_bool("pos_sig", pos_sig);
}
+GINAC_BIND_UNARCHIVER(minkmetric);
void minkmetric::archive(archive_node &n) const
{
n.add_bool("pos_sig", pos_sig);
}
-tensepsilon::tensepsilon(const archive_node &n, lst &sym_lst) : inherited(n, sym_lst)
+void tensepsilon::read_archive(const archive_node& n, lst& sym_lst)
{
+ inherited::read_archive(n, sym_lst);
n.find_bool("minkowski", minkowski);
n.find_bool("pos_sig", pos_sig);
}
+GINAC_BIND_UNARCHIVER(tensepsilon);
void tensepsilon::archive(archive_node &n) const
{
n.add_bool("pos_sig", pos_sig);
}
+GINAC_BIND_UNARCHIVER(tensdelta);
+GINAC_BIND_UNARCHIVER(tensmetric);
+GINAC_BIND_UNARCHIVER(spinmetric);
+
//////////
// functions overriding virtual functions from base classes
//////////
if (minkowski) {
for (size_t j=1; j<i.nops(); j++) {
const ex & x = i.op(j);
- if (!is_a<varidx>(x))
+ if (!is_a<varidx>(x)) {
throw(std::runtime_error("indices of epsilon tensor in Minkowski space must be of type varidx"));
- if (ex_to<varidx>(x).is_covariant())
- if (ex_to<idx>(x).get_value().is_zero())
+ }
+ if (ex_to<varidx>(x).is_covariant()) {
+ if (ex_to<idx>(x).get_value().is_zero()) {
sign = (pos_sig ? -sign : sign);
- else
+ }
+ else {
sign = (pos_sig ? sign : -sign);
+ }
+ }
}
}
ex delta_tensor(const ex & i1, const ex & i2)
{
- static ex delta = (new tensdelta)->setflag(status_flags::dynallocated);
+ static ex delta = dynallocate<tensdelta>();
if (!is_a<idx>(i1) || !is_a<idx>(i2))
throw(std::invalid_argument("indices of delta tensor must be of type idx"));
ex metric_tensor(const ex & i1, const ex & i2)
{
- static ex metric = (new tensmetric)->setflag(status_flags::dynallocated);
+ static ex metric = dynallocate<tensmetric>();
if (!is_a<varidx>(i1) || !is_a<varidx>(i2))
throw(std::invalid_argument("indices of metric tensor must be of type varidx"));
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);
+ static ex metric_neg = dynallocate<minkmetric>(false);
+ static ex metric_pos = dynallocate<minkmetric>(true);
if (!is_a<varidx>(i1) || !is_a<varidx>(i2))
throw(std::invalid_argument("indices of metric tensor must be of type varidx"));
ex spinor_metric(const ex & i1, const ex & i2)
{
- static ex metric = (new spinmetric)->setflag(status_flags::dynallocated);
+ static ex metric = dynallocate<spinmetric>();
if (!is_a<spinidx>(i1) || !is_a<spinidx>(i2))
throw(std::invalid_argument("indices of spinor metric must be of type spinidx"));
ex epsilon_tensor(const ex & i1, const ex & i2)
{
- static ex epsilon = (new tensepsilon)->setflag(status_flags::dynallocated);
+ static ex epsilon = dynallocate<tensepsilon>();
if (!is_a<idx>(i1) || !is_a<idx>(i2))
throw(std::invalid_argument("indices of epsilon tensor must be of type idx"));
ex epsilon_tensor(const ex & i1, const ex & i2, const ex & i3)
{
- static ex epsilon = (new tensepsilon)->setflag(status_flags::dynallocated);
+ static ex epsilon = dynallocate<tensepsilon>();
if (!is_a<idx>(i1) || !is_a<idx>(i2) || !is_a<idx>(i3))
throw(std::invalid_argument("indices of epsilon tensor must be of type idx"));
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);
+ static ex epsilon_neg = dynallocate<tensepsilon>(true, false);
+ static ex epsilon_pos = dynallocate<tensepsilon>(true, true);
if (!is_a<varidx>(i1) || !is_a<varidx>(i2) || !is_a<varidx>(i3) || !is_a<varidx>(i4))
throw(std::invalid_argument("indices of Lorentz epsilon tensor must be of type varidx"));