Added integral class.

[ginac.git] / ginac / indexed.cpp

diff --git a/ginac/indexed.cpp b/ginac/indexed.cpp
index d693373cae62c6cd622495d05722aa9ce7065d99..29dd86ae5e001f7c8903ce9861c702d95f90d9ec 100644 (file)
--- a/ginac/indexed.cpp
+++ b/ginac/indexed.cpp
@@ -36,6 +36,7 @@
 #include "lst.h"
 #include "archive.h"
 #include "utils.h"
 #include "lst.h"
 #include "archive.h"
 #include "utils.h"

+#include "integral.h"

 
 namespace GiNaC {
 
 
 namespace GiNaC {
 


@@ -522,6 +523,13 @@ exvector power::get_free_indices() const
                return basis_indices;
 }
 
                return basis_indices;
 }
 

+exvector integral::get_free_indices() const
+{
+       if (a.get_free_indices().size() || b.get_free_indices().size())
+               throw (std::runtime_error("integral::get_free_indices: boundary values should not have free indices"));
+       return f.get_free_indices();
+}
+

 /** Rename dummy indices in an expression.
  *
  *  @param e Expression to work on
 /** Rename dummy indices in an expression.
  *
  *  @param e Expression to work on