Changeset 1064 in Sophya for trunk/SophyaLib/HiStats/histos2.cc

Timestamp:

Jul 11, 2000, 8:57:25 PM (25 years ago)

Author:

ansari

Message:

pour PutValues,Error etc ... on prend maintenant le mini

du nbre d'elements du vecteur/matrice et de l'histo.

cmv 11/07/00

File:

• trunk/SophyaLib/HiStats/histos2.cc (modified) (5 diffs)

trunk/SophyaLib/HiStats/histos2.cc

@@ -559 +559 @@
 {
 int i,j;
-if(v.NRows()!=(uint_4)nx || v.NCol()!=(uint_4)ny) THROW(sizeMismatchErr);
-for(i=0;i<nx;i++) for(j=0;j<ny;j++) {
+//if(v.NRows()!=(uint_4)nx || v.NCol()!=(uint_4)ny) THROW(sizeMismatchErr);
+uint_4 nnx = (v.NRows()<(uint_4)nx)? v.NRows(): (uint_4)nx;
+uint_4 nny = (v.NCol() <(uint_4)ny)? v.NCol() : (uint_4)ny;
+if(nnx>0 && nny>0) for(i=0;i<nnx;i++) for(j=0;j<nny;j++) {
   (*this)(i,j) = v(i,j);
   if(err2 && ierr) Error2(i,j) = fabs(v(i,j));
@@ -573 +575 @@
 {
 int i,j;
-if(v.NRows()!=(uint_4)nx || v.NCol()!=(uint_4)ny) THROW(sizeMismatchErr);
-for(i=0;i<nx;i++) for(j=0;j<ny;j++) {
+//if(v.NRows()!=(uint_4)nx || v.NCol()!=(uint_4)ny) THROW(sizeMismatchErr);
+uint_4 nnx = (v.NRows()<(uint_4)nx)? v.NRows(): (uint_4)nx;
+uint_4 nny = (v.NCol() <(uint_4)ny)? v.NCol() : (uint_4)ny;
+if(nnx>0 && nny>0) for(i=0;i<nnx;i++) for(j=0;j<nny;j++) {
   (*this)(i,j) += v(i,j);
   if(err2 && ierr) Error2(i,j) += fabs(v(i,j));
@@ -588 +592 @@
 {
 int i,j;
-if(v.NRows()!=(uint_4)nx || v.NCol()!=(uint_4)ny) THROW(sizeMismatchErr);
-if(!err2) Errors();
-for(i=0;i<nx;i++) for(j=0;j<ny;j++) Error2(i,j) = v(i,j);
+//if(v.NRows()!=(uint_4)nx || v.NCol()!=(uint_4)ny) THROW(sizeMismatchErr);
+uint_4 nnx = (v.NRows()<(uint_4)nx)? v.NRows(): (uint_4)nx;
+uint_4 nny = (v.NCol() <(uint_4)ny)? v.NCol() : (uint_4)ny;
+if(nnx>0 && nny>0) {
+  if(!err2) Errors();
+  for(i=0;i<nnx;i++) for(j=0;j<nny;j++) Error2(i,j) = v(i,j);
+}
 return;
 }
@@ -601 +609 @@
 {
 int i,j;
-if(v.NRows()!=(uint_4)nx || v.NCol()!=(uint_4)ny) THROW(sizeMismatchErr);
-if(!err2) Errors();
-for(i=0;i<nx;i++) for(j=0;j<ny;j++)
+//if(v.NRows()!=(uint_4)nx || v.NCol()!=(uint_4)ny) THROW(sizeMismatchErr);
+uint_4 nnx = (v.NRows()<(uint_4)nx)? v.NRows(): (uint_4)nx;
+uint_4 nny = (v.NCol() <(uint_4)ny)? v.NCol() : (uint_4)ny;
+if(nnx>0 && nny>0) {
+  if(!err2) Errors();
+  for(i=0;i<nnx;i++) for(j=0;j<nny;j++)
           if(v(i,j)>0.) Error2(i,j) += v(i,j);
+}
 return;
 }
@@ -614 +626 @@
 {
 int i,j;
-if(v.NRows()!=(uint_4)nx || v.NCol()!=(uint_4)ny) THROW(sizeMismatchErr);
-if(!err2) Errors();
-for(i=0;i<nx;i++) for(j=0;j<ny;j++)
-  if(v(i,j)>0.) Error2(i,j)=v(i,j)*v(i,j); else Error2(i,j)= -v(i,j)*v(i,j);
+//if(v.NRows()!=(uint_4)nx || v.NCol()!=(uint_4)ny) THROW(sizeMismatchErr);
+uint_4 nnx = (v.NRows()<(uint_4)nx)? v.NRows(): (uint_4)nx;
+uint_4 nny = (v.NCol() <(uint_4)ny)? v.NCol() : (uint_4)ny;
+if(nnx>0 && nny>0) {
+  if(!err2) Errors();
+  for(i=0;i<nnx;i++) for(j=0;j<nny;j++)
+    if(v(i,j)>0.) Error2(i,j)=v(i,j)*v(i,j); else Error2(i,j)= -v(i,j)*v(i,j);
+}
 return;
 }