20 theLabel(aLabel), theValue(aValue), thePrecision(aPrec), theDownWeight(
21 1.), thePrediction(0.), theParameters(), theDerivatives(), globalLabels(), globalDerivatives() {
42 const std::vector<unsigned int> &labDer,
const SMatrix55 &matDer,
43 unsigned int iOff,
const TMatrixD &derLocal,
44 const std::vector<int> &labGlobal,
const TMatrixD &derGlobal,
45 unsigned int extOff,
const TMatrixD &extDer) {
47 unsigned int nParMax = 5 + derLocal.GetNcols() + extDer.GetNcols();
51 for (
int i = 0; i < derLocal.GetNcols(); ++i)
53 if (derLocal(iRow - iOff, i)) {
59 for (
int i = 0; i < extDer.GetNcols(); ++i)
61 if (extDer(iRow - iOff, i)) {
67 for (
unsigned int i = 0; i < 5; ++i)
69 if (labDer[i] and matDer(iRow, i)) {
76 for (
int i = 0; i < derGlobal.GetNcols(); ++i)
90 const std::vector<unsigned int> &labDer,
const SMatrix27 &matDer,
91 unsigned int extOff,
const TMatrixD &extDer) {
93 unsigned int nParMax = 7 + extDer.GetNcols();
97 for (
int i = 0; i < extDer.GetNcols(); ++i)
99 if (extDer(iRow, i)) {
105 for (
unsigned int i = 0; i < 7; ++i)
107 if (labDer[i] and matDer(iRow, i)) {
121 const std::vector<double> &derivatives) {
122 for (
unsigned int i = 0; i < derivatives.size(); ++i)
124 if (derivatives[i]) {
150 if (scaledResidual < 4.6851) {
151 aWeight = (1.0 - 0.045558 * scaledResidual * scaledResidual);
156 }
else if (aMethod == 2)
158 if (scaledResidual >= 1.345) {
159 aWeight = 1.345 / scaledResidual;
161 }
else if (aMethod == 3)
163 aWeight = 1.0 / (1.0 + (scaledResidual * scaledResidual / 5.6877));
187 std::cout << std::endl;
192 std::cout << std::endl;
203 std::vector<unsigned int>* &indLocal, std::vector<double>* &derLocal) {
221 std::vector<unsigned int>* &indLocal, std::vector<double>* &derLocal,
222 std::vector<int>* &labGlobal, std::vector<double>* &derGlobal) {
240 double &aDownWeight, std::vector<unsigned int>* &indLocal,
241 std::vector<double>* &derLocal) {