jzrake

from OpenGL.GLUT import *
from OpenGL.GLU import *
from OpenGL.GL import *


rotationX = 0.0
rotationY = 0.0
last_x = 0
last_y = 0

jzrake

local function plot(series, tpause)
   local gp = io.popen("gnuplot", 'w')

   local lines = { }
   for k,v in pairs(series) do
      table.insert(lines, string.format(" '-' u 1:2 title '%s'", k))
   end
   
   gp:write("plot" .. table.concat(lines, ",") .. "\n")
   for k,v in pairs(series) do

jzrake

#!/usr/bin/env python

import numpy as np


def Quantile(data, q, precision=1.0):
    """
    Returns the q'th percentile of the distribution given in the argument
    'data'. Uses the 'precision' parameter to control the noise level.
    """

jzrake

#!/usr/bin/env python


def load_from_many(pattern):
    """
    Take a pattern and return a numpy array containing the concatenation of all
    the files matching that pattern. Only the first column and rows after the
    first are returned, and the data values are offset by the difference between
    the first and second columns of the first row.
    """

jzrake

#!/usr/bin/env python

import numpy as np
from matplotlib import pyplot as plt
from scipy.optimize import leastsq
from scipy.ndimage import binary_closing, grey_closing


def LineModel(pt1, pt2, Nx=100, Ny=100, h=1.0, sig=0.1):
    """

jzrake

void draw_vbo()
{
    GLuint vbo,ibo;
    GLfloat verts[8][3] = {{0.0, 0.0, 0.0},
			   {0.0, 0.0, 0.1},
			   {0.0, 0.1, 0.0},
			   {0.0, 0.1, 0.1},
			   {0.1, 0.0, 0.0},
			   {0.1, 0.0, 0.1},
			   {0.1, 0.1, 0.0},

jzrake

#!/usr/bin/env python

# ------------------------------------------------------------------------------
#
# Authors: Jonathan Zrake, Bez Laderman: NYU CCPP
# Date: May 7th, 2012
#
# This piece of code implements the left and right eigenvectors of the ideal
# special relativistic hydrodynamics equations. The formulas are an exact
# translation of those given in the literature:

jzrake

#include <iostream>
#include <stdexcept>


class ZrakeFailed : public std::runtime_error
{
public:
  ZrakeFailed(const std::string &msg) : std::runtime_error("zraker failed. " + msg) { }

jzrake

/*------------------------------------------------------------------------------
 * FILE: euler.c
 *
 * Copyright (C) 2011 Jonathan Zrake, NYU CCPP: zrake <at> nyu <dot> edu
 *
 * This program is free software: you can redistribute it and/or modify it under
 * the terms of the GNU General Public License as published by the Free Software
 * Foundation, either version 3 of the License, or (at your option) any later
 * version.

jzrake

int H5Lexists_safe(hid_t base, const char *path)
// -----------------------------------------------------------------------------
// The HDF5 specification only allows H5Lexists to be called on an immediate
// child of the current object. However, you may wish to see whether a whole
// relative path exists, returning false if any of the intermediate links are
// not present. This function does that.
// http://www.hdfgroup.org/HDF5/doc/RM/RM_H5L.html#Link-Exists
// -----------------------------------------------------------------------------
{