79spinmetric::spinmetric()
101 inherited::read_archive(
n, sym_lst);
108 inherited::archive(
n);
114 inherited::read_archive(
n, sym_lst);
122 inherited::archive(
n);
164 return inherited::compare_same_type(other);
185 return inherited::compare_same_type(other);
234 }
catch (std::exception &e) {
307 return inherited::eval_indexed(i);
323 if (!(
static_cast<const indexed &
>(i).get_dummy_indices().empty()))
349 if (!(
static_cast<const indexed &
>(i).get_dummy_indices().empty()))
358 v.reserve(i.
nops() - 1);
359 for (
size_t j=1; j<i.
nops(); j++)
365 for (
size_t j=1; j<i.
nops(); j++) {
366 const ex &
x = i.
op(j);
368 throw(std::runtime_error(
"indices of epsilon tensor in Minkowski space must be of type varidx"));
372 sign = (
pos_sig ? -sign : sign);
375 sign = (
pos_sig ? sign : -sign);
393 bool first_index_tried =
false;
397 for (
size_t i=1; i<other->nops(); i++) {
411 }
catch (std::exception &e) {
418 if (!first_index_tried) {
423 first_index_tried =
true;
509 bool first_index_tried =
false;
514 for (
size_t i=1; i<other->nops(); i++) {
521 *other = other->subs(other_idx == *free_idx);
528 if (!first_index_tried) {
533 first_index_tried =
true;
547 size_t num = self->nops() - 1;
554 for (
size_t i=0; i<num; i++) {
555 for (
size_t j=0; j<num; j++) {
582 throw(std::invalid_argument(
"indices of delta tensor must be of type idx"));
592 throw(std::invalid_argument(
"indices of metric tensor must be of type varidx"));
603 throw(std::invalid_argument(
"indices of metric tensor must be of type varidx"));
613 throw(std::invalid_argument(
"indices of spinor metric must be of type spinidx"));
615 throw(std::runtime_error(
"index dimension for spinor metric must be 2"));
625 throw(std::invalid_argument(
"indices of epsilon tensor must be of type idx"));
629 throw(std::invalid_argument(
"all indices of epsilon tensor must have the same dimension"));
631 throw(std::runtime_error(
"index dimension of epsilon tensor must match number of indices"));
644 throw(std::invalid_argument(
"indices of epsilon tensor must be of type idx"));
648 throw(std::invalid_argument(
"all indices of epsilon tensor must have the same dimension"));
650 throw(std::runtime_error(
"index dimension of epsilon tensor must match number of indices"));
664 throw(std::invalid_argument(
"indices of Lorentz epsilon tensor must be of type varidx"));
668 throw(std::invalid_argument(
"all indices of epsilon tensor must have the same dimension"));
670 throw(std::runtime_error(
"index dimension of epsilon tensor must match number of indices"));
Archiving of GiNaC expressions.
#define GINAC_ASSERT(X)
Assertion macro for checking invariances.
This class stores all properties needed to record/retrieve the state of one object of class basic (or...
This class is the ABC (abstract base class) of GiNaC's class hierarchy.
virtual size_t nops() const
Number of operands/members.
virtual ex op(size_t i) const
Return operand/member at position i.
const basic & hold() const
Stop further evaluation.
virtual ex subs(const exmap &m, unsigned options=0) const
Substitute a set of objects by arbitrary expressions.
virtual int compare_same_type(const basic &other) const
Returns order relation between two objects of same type.
Wrapper template for making GiNaC classes out of STL containers.
Lightweight wrapper for GiNaC's symbolic objects.
bool is_equal(const ex &other) const
ex simplify_indexed(unsigned options=0) const
Simplify/canonicalize expression containing indexed objects.
ex subs(const exmap &m, unsigned options=0) const
This class holds one index of an indexed object.
ex replace_dim(const ex &new_dim) const
Make a new index with the same value but a different dimension.
ex get_dim() const
Get dimension of index space.
bool is_symbolic() const
Check whether the index is symbolic.
ex minimal_dim(const idx &other) const
Return the minimum of the dimensions of this and another index.
ex get_value() const
Get value of index.
This class holds an indexed expression.
ex determinant(unsigned algo=determinant_algo::automatic) const
Determinant of square matrix.
This class represents a Minkowski metric tensor.
minkmetric(bool pos_sig)
Construct Lorentz metric tensor with given signature.
ex eval_indexed(const basic &i) const override
Automatic symbolic evaluation of an indexed Lorentz metric tensor.
bool pos_sig
If true, the metric is diag(-1,1,1...).
void archive(archive_node &n) const override
Save (a.k.a.
void read_archive(const archive_node &n, lst &syms) override
Read (a.k.a.
bool info(unsigned inf) const override
Information about the object.
Context for default (ginsh-parsable) output.
Context for latex-parsable output.
This class holds a spinor index that can be dotted or undotted and that also has a variance.
This class represents an antisymmetric spinor metric tensor which can be used to raise/lower indices ...
bool info(unsigned inf) const override
Information about the object.
ex eval_indexed(const basic &i) const override
Automatic symbolic evaluation of an indexed metric tensor.
bool contract_with(exvector::iterator self, exvector::iterator other, exvector &v) const override
Contraction of an indexed spinor metric with something else.
@ expanded
.expand(0) has already done its job (other expand() options ignore this flag)
@ evaluated
.eval() has already done its job
@ no_pattern
disable pattern matching
This class represents the delta tensor.
bool contract_with(exvector::iterator self, exvector::iterator other, exvector &v) const override
Contraction of an indexed delta tensor with something else.
This class represents the totally antisymmetric epsilon tensor.
bool minkowski
If true, tensor is in Minkowski-type space.
bool pos_sig
If true, the metric is assumed to be diag(-1,1,1...).
void read_archive(const archive_node &n, lst &syms) override
Read (a.k.a.
bool contract_with(exvector::iterator self, exvector::iterator other, exvector &v) const override
Contraction of epsilon tensor with something else.
bool info(unsigned inf) const override
Information about the object.
tensepsilon(bool minkowski, bool pos_sig)
void archive(archive_node &n) const override
Save (a.k.a.
ex eval_indexed(const basic &i) const override
Automatic symbolic evaluation of an indexed epsilon tensor.
This class represents a general metric tensor which can be used to raise/lower indices.
bool info(unsigned inf) const override
Information about the object.
ex eval_indexed(const basic &i) const override
Automatic symbolic evaluation of an indexed metric tensor.
bool contract_with(exvector::iterator self, exvector::iterator other, exvector &v) const override
Contraction of an indexed metric tensor with something else.
This class holds one of GiNaC's predefined special tensors such as the delta and the metric tensors.
bool replace_contr_index(exvector::iterator self, exvector::iterator other) const
Replace dummy index in contracted-with object by the contracting object's second index (used internal...
This class holds an index with a variance (co- or contravariant).
bool is_covariant() const
Check whether the index is covariant.
Interface to GiNaC's indices.
Interface to GiNaC's indexed expressions.
Definition of GiNaC's lst.
Interface to symbolic matrices.
const symmetry & antisymmetric4()
bool is_zero(const ex &thisex)
ex spinor_metric(const ex &i1, const ex &i2)
Create a spinor metric tensor with specified indices.
ex metric_tensor(const ex &i1, const ex &i2)
Create a symmetric metric tensor with specified indices.
std::map< ex, ex, ex_is_less > exmap
const T & ex_to(const ex &e)
Return a reference to the basic-derived class T object embedded in an expression.
B & dynallocate(Args &&... args)
Constructs a new (class basic or derived) B object on the heap.
const symmetry & antisymmetric3()
const symmetry & antisymmetric2()
GINAC_IMPLEMENT_REGISTERED_CLASS_OPT_T(lst, basic, print_func< print_context >(&lst::do_print). print_func< print_tree >(&lst::do_print_tree)) template<> bool ls GINAC_BIND_UNARCHIVER)(lst)
Specialization of container::info() for lst.
const symmetry & symmetric2()
bool is_dummy_pair(const idx &i1, const idx &i2)
Check whether two indices form a dummy pair.
bool is_a(const basic &obj)
Check if obj is a T, including base classes.
ex delta_tensor(const ex &i1, const ex &i2)
Create a delta tensor with specified indices.
ex lorentz_eps(const ex &i1, const ex &i2, const ex &i3, const ex &i4, bool pos_sig)
Create an epsilon tensor in a Minkowski space with four indices.
ex lorentz_g(const ex &i1, const ex &i2, bool pos_sig)
Create a Minkowski metric tensor with specified indices.
bool is_exactly_a(const basic &obj)
Check if obj is a T, not including base classes.
int permutation_sign(It first, It last)
std::vector< ex > exvector
ex epsilon_tensor(const ex &i1, const ex &i2)
Create an epsilon tensor in a Euclidean space with two indices.
Makes the interface to the underlying bignum package available.
Interface to GiNaC's overloaded operators.
#define GINAC_IMPLEMENT_REGISTERED_CLASS_OPT(classname, supername, options)
Macro for inclusion in the implementation of each registered class.
#define GINAC_IMPLEMENT_REGISTERED_CLASS(classname, supername)
Macro for inclusion in the implementation of each registered class.
Interface to relations between expressions.
Interface to GiNaC's symmetry definitions.
Interface to GiNaC's special tensors.
Interface to several small and furry utilities needed within GiNaC but not of any interest to the use...
#define DEFAULT_PRINT_LATEX(classname, text, latex)
#define DEFAULT_PRINT(classname, text)
#define DEFAULT_CTOR(classname)
#define DEFAULT_COMPARE(classname)