JohannesBuchner

#!/usr/bin/env python3
"""
Manage a group of processes, keeping only N running at a time.

Usage:
  procmanage.py N STRING

Arguments:
  N            An integer representing the number of processes to send the CONT signal to.
  STRING       A string representing the full name to match against.

JohannesBuchner

#!/usr/bin/env python3
# 
# SYNOPSIS: stripgif.py input.odp output.odp
# 
# 

import os
import sys
import zipfile
import lxml.etree as ET

JohannesBuchner

import matplotlib as mpl
import matplotlib.pyplot as plt

# make a colormap that can be used to look up colors
norm = mpl.colors.Normalize(vmin=0, vmax=3)
zcmap = plt.cm.ScalarMappable(norm=norm, cmap=plt.cm.viridis)

# use it to plot stuff
plt.errorbar(
    x=1, xerr=0.1, y=2, yerr=0.2,

JohannesBuchner

def sbpl(x, norm, lam1, lam2, x0, xbrk, Lambda):
    """Smoothly bending powerlaw

    Parameterization from Ryde99
    https://ui.adsabs.harvard.edu/abs/1999ApL%26C..39..281R/abstract

    Parameters
    ----------
    x: array of independent variable
    xmax: only consider x values up to this value

JohannesBuchner

from getdist.mcsamples import MCSamples
import getdist.chains

def highest_density_interval_from_samples(xsamples, xlo=None, xhi=None, probability_level=0.68):
    """
    Compute the highest density interval (HDI) from posterior samples.

    Parameters
    ----------
    xsamples : array_like

JohannesBuchner

import warnings
import sys
import astropy.io.fits as pyfits
from astropy.table import Table

print("opening FITS file ...")
f = pyfits.open(sys.argv[1])

print("extensions: ", [e.name for e in f])