sergpolly

    header, body_stream = _headerops.get_header(instream)
    # # turn body_stream into file-like object aka stream:
    body_stream = _io.StringIteratorIO(body_stream)
    # # proceed ...
    header = _headerops.append_new_pg(header, ID=UTIL_NAME, PN=UTIL_NAME)
    outstream.writelines((l+'\n' for l in header))
    if outstream_rest:
        outstream_rest.writelines((l+'\n' for l in header))

    column_names = _headerops.extract_column_names(header)

sergpolly

def get_header(instream, comment_char='#'):
    '''Returns a header from the stream and a stream with the remaining
    lines.
    Parameters
    ----------
    instream : a file object
        An input stream.
    comment_char : str
        The character prepended to header lines (use '@' when parsing sams, 
        '#' when parsing pairsams).

sergpolly

package juicebox.tools.utils.juicer.hiccups;

import com.google.common.primitives.Doubles;
import com.google.common.primitives.Floats;
import jcuda.Pointer;
import jcuda.Sizeof;
import jcuda.driver.CUdeviceptr;
import jcuda.utils.KernelLauncher;
import juicebox.data.HiCFileTools;
import juicebox.data.MatrixZoomData;

sergpolly

sergpolly

sergpolly

sergpolly

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 66,
   "metadata": {},
   "outputs": [],
   "source": [
    "import numpy as np\n",
    "import matplotlib.pyplot as plt\n",

sergpolly

sergpolly

sergpolly