| 129 | **********************************************************************/ |
| 130 | |
| 131 | double LLSQ::rms(double m, double c) const { // get error |
| 132 | double error; // total error |
| 133 | |
| 134 | if (total_weight > 0) { |
| 135 | error = sigyy + m * (m * sigxx + 2 * (c * sigx - sigxy)) + c * |
| 136 | (total_weight * c - 2 * sigy); |
| 137 | if (error >= 0) |
| 138 | error = sqrt(error / total_weight); // sqrt of mean |
| 139 | else |
| 140 | error = 0; |
| 141 | } else { |
| 142 | error = 0; // too little |
| 143 | } |
| 144 | return error; |
| 145 | } |
| 146 | |
| 147 | |
| 148 | /********************************************************************** |
no test coverage detected