rawlik

rawlik

import pylab as pl
from scipy.misc import derivative
import inspect

def cramer_rao(model, p0, X, noise, show_plot=False):
    """Calulate inverse of the Fisher information matrix for model
    sampled on grid X with parameters p0. Assumes samples are not
    correlated and have equal variance noise^2.
    
    Parameters

rawlik

rawlik

rawlik

CC=clang++

FLAGS= `root-config --cflags` -std=c++11
LINK_FLAGS= `root-config --glibs`

DEPS = breedSFC.h
PROGRAMS = playground tests
PROGRAMS_OBJS = $(foreach X,$(PROGRAMS),$(X).o)
OBJS = breedSFC.o

rawlik

class careful_autoscale:
    def __init__(self, ax=None):
        self.ax = ax if ax else gca()
    
    def __enter__(self):
        self.xl = self.ax.get_xlim()
        self.yl = self.ax.get_ylim()
        self.lines = self.ax.get_lines()
        self.lines_visibility = [ l.get_visible() for l in self.lines ]
        [ l.set_visible(False) for l in self.lines ]

rawlik

def brighten(color, v):
    return np.array(matplotlib.colors.to_rgb(color)) ** (1 - v)

rawlik

# use color ANSI escape sequences on windows
if os.name == "nt":
    import colorama
    colorama.init()

def welcome():
    BLACK   = "\033[0;30m"
    RED     = "\033[0;31m"
    GREEN   = "\033[0;32m"
    YELLOW  = "\033[0;33m"

rawlik

import coloredlogs
import logging

if __name__=="__main__":
    # only configure logging in the main program, everywhere else just import and use it
    
    # set up logging
    logfile_folder = "N:/data/log/"
    logfile_date = datetime.datetime.now().strftime("%Y%m%d_%H%M%S")
    logfile_filename = f"bunker5_{logfile_date}.log"

rawlik