debmsg

pub/src/plot.cpp

@@ -290,6 +290,7 @@ void NPlot::Plot( class CPlot *plot, bool add )
 
                     // evaluate curve for grid xc:
                     fc->curve_val_vec( &yc, xc, k, j );
+                    cout << "DEB evaluation -> " << xc.size() << "\n";
 
                     // merge xc,yc and xn,yn:
                     if( iref==0 ) {
@@ -311,6 +312,7 @@ void NPlot::Plot( class CPlot *plot, bool add )
                         }
                         xn = xa; yn = ya;
                     }
+                    cout << "DEB merger -> " << xn.size() << "\n";
 
                     // end of loop
                     if( iref>=10 || xn.size()>=plot->maxpoints )
@@ -333,6 +335,7 @@ void NPlot::Plot( class CPlot *plot, bool add )
                             ++i;
                     }
                     xn = xa; yn = ya;
+                    cout << "DEB within frame -> " << xn.size() << "\n";
 
                     // set additional base points:
                     xc.clear();
@@ -364,8 +367,9 @@ void NPlot::Plot( class CPlot *plot, bool add )
                         xa.push_back( xn[i] );
                         ya.push_back( yn[i] );
                     } 
-                    if( plot->Y.contains(yc[i]) || i==npts-1 ) {
+                    if( !plot->Y.contains(yc[i]) || i==npts-1 ) {
                         if( xa.size() ){
+                            cout << "DEB segment -> " << xa.size() << "\n";
                             xpv.push_back(xa);
                             ypv.push_back(ya);
                             xa.clear();
@@ -373,6 +377,7 @@ void NPlot::Plot( class CPlot *plot, bool add )
                         }
                     }
                 }
+                cout << "DEB #segments " << xpv.size() << "\n";
                 
                 // plot subcurves: TODO merge with above, eliminate xpv,ypv
                 for ( uint ipv=0; ipv<xpv.size();++ipv) {