Bug nicht reproduzierbar.

Evtl. Kompilation ausser Takt.

TODO

 == BUGS ==
 
-  oi x[,,ni-1]-x[,,0]
-index j=102 out of bounds
-  oi x[,,0]-x[,,ni-1]
-i out of range
-  oi x[,j,ni-1]-x[,j,0]
-index j=4294967295 out of bounds
-
-noch einfacher:
-fm 9 1 \
-x [,,ni-1]
-index j=4294967295 out of bounds
-
 - coord name algebra
 - convolutand must not be defined over full energy range
 - we need kww with double precision, period.

pub/src/expr.cpp

@@ -22,6 +22,7 @@
 
 extern bool allow_slow_conv; // TODO unify parameter treatment
 
+#define DEBUG 0
 
 //***************************************************************************//
 //* class CContext                                                          *//
@@ -221,8 +222,10 @@ void CRef::ref_val( CResult& ret, const CContext& ctx ) const
     } else {
         j = ctx.j;
     }
-    if ( j>=f->nJ() )
-        throw string( "index j=" + strg(j) + " out of bounds" );
+    if ( j==-1 )
+        throw "index j not defined";
+    else if ( j>=f->nJ() )
+        throw "index j=" + strg(j) + " out of bounds";
 
     // Return k-j-dependent reference (unless done above):
     if        ( typ == _J ) {
@@ -252,7 +255,8 @@ void CRef::ref_val( CResult& ret, const CContext& ctx ) const
                         throw string( "missing index i" );
                 }
                 if( i>= sj->size() )
-                    throw string( "i out of range" );
+                    throw "i=" + strg(i) + " exceeds scan size" +
+                        strg(sj->size());
                 if       ( typ == _I ) {
                     ret.set_r( i );
                 } else if ( typ == _X ) {
@@ -276,7 +280,8 @@ void CRef::ref_val( CResult& ret, const CContext& ctx ) const
                         uint i;
                         ti->tree_uival( &i, myctx );
                         if( i >= sj->size() )
-                            throw string( "i out of range" );
+                            throw "i=" + strg(i) + " exceeds scan size" +
+                                strg(sj->size());
                         if        ( typ == _I ) {
                             ret.v[ii] = i;
                         } else if ( typ == _X ) {
@@ -741,6 +746,12 @@ void CTree::tree_op( CResult& ret, const CContext& ctx ) const
             }
         }
     }
+    if( DEBUG ) {
+        cout << "tree_op:\n";
+        for ( uint iarg=0; iarg<narg; ++iarg )
+            cout << "  arg: " << a[iarg].result_info() << "\n";
+        cout << "  res: " << ret.result_info() << "\n";
+    }
 }
 
 
@@ -761,7 +772,7 @@ void CTree::tree_cev( CResult& ret, const CContext& ctx ) const
     }
     if ( k==(uint)-1 )
         throw "missing or invalid k";
-    if ( k>=NOlm::MOM.size())
+    else if ( k>=NOlm::MOM.size())
         throw "index k=" + strg(k) + " out of bounds";
 
     // k points to a curve: