78spinmetric::spinmetric()
100 inherited::read_archive(n, sym_lst);
107 inherited::archive(n);
113 inherited::read_archive(n, sym_lst);
121 inherited::archive(n);
163 return inherited::compare_same_type(other);
184 return inherited::compare_same_type(other);
233 }
catch (std::exception &e) {
306 return inherited::eval_indexed(i);
322 if (!(
static_cast<const indexed &
>(i).get_dummy_indices().empty()))
348 if (!(
static_cast<const indexed &
>(i).get_dummy_indices().empty()))
357 v.reserve(i.
nops() - 1);
358 for (
size_t j=1; j<i.
nops(); j++)
364 for (
size_t j=1; j<i.
nops(); j++) {
365 const ex & x = i.
op(j);
367 throw(std::runtime_error(
"indices of epsilon tensor in Minkowski space must be of type varidx"));
371 sign = (
pos_sig ? -sign : sign);
374 sign = (
pos_sig ? sign : -sign);
392 bool first_index_tried =
false;
396 for (
size_t i=1; i<other->nops(); i++) {
410 }
catch (std::exception &e) {
417 if (!first_index_tried) {
422 first_index_tried =
true;
508 bool first_index_tried =
false;
513 for (
size_t i=1; i<other->nops(); i++) {
520 *other = other->subs(other_idx == *free_idx);
527 if (!first_index_tried) {
532 first_index_tried =
true;
546 size_t num = self->nops() - 1;
553 for (
size_t i=0; i<num; i++) {
554 for (
size_t j=0; j<num; j++) {
581 throw(std::invalid_argument(
"indices of delta tensor must be of type idx"));
591 throw(std::invalid_argument(
"indices of metric tensor must be of type varidx"));
602 throw(std::invalid_argument(
"indices of metric tensor must be of type varidx"));
612 throw(std::invalid_argument(
"indices of spinor metric must be of type spinidx"));
614 throw(std::runtime_error(
"index dimension for spinor metric must be 2"));
624 throw(std::invalid_argument(
"indices of epsilon tensor must be of type idx"));
628 throw(std::invalid_argument(
"all indices of epsilon tensor must have the same dimension"));
630 throw(std::runtime_error(
"index dimension of epsilon tensor must match number of indices"));
643 throw(std::invalid_argument(
"indices of epsilon tensor must be of type idx"));
647 throw(std::invalid_argument(
"all indices of epsilon tensor must have the same dimension"));
649 throw(std::runtime_error(
"index dimension of epsilon tensor must match number of indices"));
663 throw(std::invalid_argument(
"indices of Lorentz epsilon tensor must be of type varidx"));
667 throw(std::invalid_argument(
"all indices of epsilon tensor must have the same dimension"));
669 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...
bool find_bool(const std::string &name, bool &ret, unsigned index=0) const
Retrieve property of type "bool" from node.
void add_bool(const std::string &name, bool value)
Add property of type "bool" to node.
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.
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.
bool info(unsigned inf) const override
Information about the object.
ex eval_indexed(const basic &i) const override
Automatic symbolic evaluation of an indexed delta tensor.
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.
container< std::list > lst
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
B & dynallocate(Args &&... args)
Constructs a new (class basic or derived) B object on the heap.
const symmetry & antisymmetric3()
attribute_pure const T & ex_to(const ex &e)
Return a reference to the basic-derived class T object embedded in an expression.
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.
GINAC_IMPLEMENT_REGISTERED_CLASS_OPT(add, expairseq, print_func< print_context >(&add::do_print). print_func< print_latex >(&add::do_print_latex). print_func< print_csrc >(&add::do_print_csrc). print_func< print_tree >(&add::do_print_tree). print_func< print_python_repr >(&add::do_print_python_repr)) add
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(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)