?

src/curve.cpp

@@ -387,13 +387,17 @@ void globalEvaluate( double* par, int m_dat, double* fvec,
         fc->T->tree_curve_val( &fit, s->x, mydata->k, mydata->j );
         
         for ( uint i=0; i<m_dat; i++) {
+            if( fit[i]==INFINITY ) // TEMPORARY UNTIL exceptions can be
+                                   // thrown everywhere
+                throw string( "function evaluation return INFINITY" );
+
             if (!fc->weight_y_log) {
                 fvec[i] = s->y[i] - fit[i];
-
+/*
                 fvec[i] *= fvec[i]; // TEMPORARY
                 fvec[i] *= fvec[i]; // TEMPORARY
                 fvec[i] *= fvec[i]; // TEMPORARY
-
+*/
             } else {
                 if ( s->y[i] < 1e-99 || fit[i] < 1e-99 ) {
                     fvec[i] = 0; // just don't fit