st-bender / sciapy

					self.times, self.values, left=0., right=0.) * taufac

		return yp * self.dt

	def _lt_corr_grad(self, t, tau, tmax=60.):

		"""Lifetime corrected gradient

		Corrects for a finite lifetime by summing over the last `tmax`

		days with an exponential decay given of lifetime(s) `taus`.

		"""

		bs = np.arange(self.dt, tmax + self.dt, self.dt)

		ypg = np.zeros_like(t)

		tauexp = np.exp(-self.dt / tau)

		taufac = np.ones_like(tau)

		for b in bs:

			taufac *= tauexp

			ypg += np.interp(t - self.lag - b / self.days_per_time_unit,

					self.times, self.values, left=0., right=0.) * taufac * b

		return ypg * self.dt / tau**2

	def get_value(self, t):

		t = np.atleast_1d(t)

				self.t_adj = -44.25

		super(ProxyModel, self).__init__(*args, **kwargs)

	def _lt_corr(self, t, tau, tmax=60.):

		"""Lifetime corrected values

		Corrects for a finite lifetime by summing over the last `tmax`

		days with an exponential decay given of lifetime(s) `taus`.

		"""

		bs = np.arange(self.dt, tmax + self.dt, self.dt)

		yp = np.zeros_like(t)

		tauexp = np.exp(-self.dt / tau)

		taufac = np.ones_like(tau)

		for b in bs:

			taufac *= tauexp

			yp += np.interp(t - self.lag - b / self.days_per_time_unit,

					self.times, self.values, left=0., right=0.) * taufac

		return yp * self.dt

	def _lt_corr_grad(self, t, tau, tmax=60.):

		"""Lifetime corrected gradient