From aae6dcb11162edfe778b093d88a6ef531c107f14 Mon Sep 17 00:00:00 2001
From: "Joachim Wuttke (laptop)" <j.wuttke@fz-juelich.de>
Date: Fri, 18 Jan 2013 10:03:11 +0100
Subject: [PATCH] com
---
pub/src/curve.cpp | 7 ++++---
1 file changed, 4 insertions(+), 3 deletions(-)
diff --git a/pub/src/curve.cpp b/pub/src/curve.cpp
index 13d3dffa..141049d5 100644
--- a/pub/src/curve.cpp
+++ b/pub/src/curve.cpp
@@ -925,11 +925,12 @@ void NCurveFile::Fit( bool _allow_slow_conv )
if( status.info==11 )
throw "exception in fit function";
- // compute fit quality metrics:
+ // Fit quality metrics Q:
+ // Q: status code
C->Quality[0] = (double) status.info;
- // - coefficient of determination R^2 (for linear regression)
+ // Q: coefficient of determination R^2 (for linear regression)
double ysum1 = 0, ysum2 = 0;
for( uint i=0; i<nd; ++i ){
ysum1 += D->y[i];
@@ -938,7 +939,7 @@ void NCurveFile::Fit( bool _allow_slow_conv )
double ydev2 = ysum2 - SQR( ysum1 )/nd; // = nd * <(y-<y>)^2>
C->Quality[1] = ydev2!=0 ? 1 - SQR(status.fnorm) / ydev2 : -1;
- // - mean deviation (chi^2 if fvec is weighed with 1/dy)
+ // Q: mean deviation (chi^2 if fvec is weighed with 1/dy)
uint nfreedom = nd - npfree - 1;
if( nfreedom<1 ) nfreedom = 1;
C->Quality[2] = SQR(status.fnorm) / nfreedom;
--
GitLab