Added integral class.
[ginac.git] / ginac / indexed.cpp
index d693373cae62c6cd622495d05722aa9ce7065d99..29dd86ae5e001f7c8903ce9861c702d95f90d9ec 100644 (file)
@@ -36,6 +36,7 @@
 #include "lst.h"
 #include "archive.h"
 #include "utils.h"
+#include "integral.h"
 
 namespace GiNaC {
 
@@ -522,6 +523,13 @@ exvector power::get_free_indices() const
                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