Spread total immune pressure strength uniformly over targeted epitopes.

diff --git a/wf.cpp b/wf.cpp
index 30a9364..0064dea 100644 (file)
--- a/wf.cpp
+++ b/wf.cpp
@@ -100,7 +100,7 @@ void run(RunParameters &r, unsigned seed) {
             
                 EpitopeHamiltonian H(r.couplingsInfile);
                 for (unsigned ep=0; ep<r.numEpitopes; ++ep)
-                    H.set_epitope(r.eWT[ep], r.eMut[ep], r.penalty);
+                    H.set_epitope(r.eWT[ep], r.eMut[ep], r.penalty/(double)r.numEpitopes);
                 H.set_temp(r.bh, r.bJ);
                 Population P(H, r.n, r.mu, r.initPop, r.initFrac);