models/degrader_constant.py [107:131]:
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
        gr = self.r * torch.sigmoid(4.0 * (t - self.tlag))

        # Specific growth and dilution
        g = 1.0 - x / self.K
        gamma = gr * g

        # Promoter activity
        boundLuxR = luxR * luxR * self.fracLuxR
        boundLasR = lasR * lasR * self.fracLasR
        P76 = (self.e76 + self.KGR_76 * boundLuxR + self.KGS_76 * boundLasR) / (
            1.0 + self.KGR_76 * boundLuxR + self.KGS_76 * boundLasR
        )
        P81 = (self.e81 + self.KGR_81 * boundLuxR + self.KGS_81 * boundLasR) / (
            1.0 + self.KGR_81 * boundLuxR + self.KGS_81 * boundLasR
        )

        # Right-hand sides
        d_x = gamma * x
        d_rfp = self.rc - (gamma + self.drfp) * rfp
        d_yfp = self.rc * self.aYFP * P81 - (gamma + self.dyfp) * yfp
        d_cfp = self.rc * self.aCFP * P76 - (gamma + self.dcfp) * cfp
        d_f530 = self.rc * self.a530 - gamma * f530
        d_f480 = self.rc * self.a480 - gamma * f480
        d_luxR = self.rc * self.aR - (gamma + self.dR) * luxR
        d_lasR = self.rc * self.aS - (gamma + self.dS) * lasR
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -



models/relay_constant.py [95:119]:
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
        gr = self.r * torch.sigmoid(4.0 * (t - self.tlag))

        # Specific growth and dilution
        g = 1.0 - x / self.K
        gamma = gr * g

        # Promoter activity
        boundLuxR = luxR * luxR * self.fracLuxR
        boundLasR = lasR * lasR * self.fracLasR
        P76 = (self.e76 + self.KGR_76 * boundLuxR + self.KGS_76 * boundLasR) / (
            1.0 + self.KGR_76 * boundLuxR + self.KGS_76 * boundLasR
        )
        P81 = (self.e81 + self.KGR_81 * boundLuxR + self.KGS_81 * boundLasR) / (
            1.0 + self.KGR_81 * boundLuxR + self.KGS_81 * boundLasR
        )

        # Right-hand sides
        d_x = gamma * x
        d_rfp = self.rc - (gamma + self.drfp) * rfp
        d_yfp = self.rc * self.aYFP * P81 - (gamma + self.dyfp) * yfp
        d_cfp = self.rc * self.aCFP * P76 - (gamma + self.dcfp) * cfp
        d_f530 = self.rc * self.a530 - gamma * f530
        d_f480 = self.rc * self.a480 - gamma * f480
        d_luxR = self.rc * self.aR - (gamma + self.dR) * luxR
        d_lasR = self.rc * self.aS - (gamma + self.dS) * lasR
- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -