renormalizing custom rr profiles as done for built-in profiles (average = 1)

bayesiancook · bayesiancook · commit 3c776560e822 · 2023-02-01T22:05:24.000+01:00
diff --git a/sources/GTRProfileProcess.cpp b/sources/GTRProfileProcess.cpp
@@ -373,6 +373,13 @@ void GTRProfileProcess::SetRR(string type)	{
 					
 				}
 			}
+            double total = 0;
+			for (int i=0; i<Nrr; i++)	{
+				total += rr[i];
+			}
+			for (int i=0; i<Nrr; i++)	{
+				rr[i] /= total / Nrr;
+			}
 		}	
 	}
 }

Original file line number	Diff line number	Diff line change
`@@ -373,6 +373,13 @@ void GTRProfileProcess::SetRR(string type) {`
`373`	`373`
`374`	`374`	`}`
`375`	`375`	`}`
	`376`	`+ double total = 0;`
	`377`	`+ for (int i=0; i<Nrr; i++) {`
	`378`	`+ total += rr[i];`
	`379`	`+ }`
	`380`	`+ for (int i=0; i<Nrr; i++) {`
	`381`	`+ rr[i] /= total / Nrr;`
	`382`	`+ }`
`376`	`383`	`}`
`377`	`384`	`}`
`378`	`385`	`}`