jtilly

% Mean utilities
mu1 = 2;
mu2 = 3;
mu3 = 1;

% Scale parameter
sigma = 0.4;

% Draw Type I extreme value errors 
R = -evrnd(0,1,100000,3);

jtilly

real(dp), PARAMETER :: NAN = TRANSFER((/ Z'00000000', Z'7FF80000' /), 1.0_8)

jtilly

#!/bin/bash
sudo apt-get install docker-engine=1.8.3-0~trusty
sudo gpasswd -a ${USER} docker
sudo service docker restart
newgrp docker

jtilly

#!/bin/bash
curl -O http://r.research.att.com/libs/gfortran-4.8.2-darwin13.tar.bz2
sudo tar fvxz gfortran-4.8.2-darwin13.tar.bz2 -C /

jtilly

#!/bin/bash
cat -e -t -v $1

jtilly

#!/bin/bash

# get the latest version
VER=$(wget --no-check-certificate -qO- https://s3.amazonaws.com/rstudio-desktop/current.ver)
# download it
wget https://download1.rstudio.org/rstudio-${VER}-amd64.deb
# install it
sudo dpkg -i rstudio-${VER}-amd64.deb
# clean up
rm -f rstudio-${VER}-amd64.deb

jtilly

---
title: 'Random Number Generation in Parallel: Reproducibility'
author: "Jan Tilly"
date: "January 15, 2016"
output: html_document
---

```{r global_options, include=FALSE}
knitr::opts_chunk$set(results = 'hold')
```

jtilly

function ols(y, x, n, k) result (beta) 

    implicit none

    external DGELS

    integer, intent(in) :: n, k
    real(dp), allocatable, intent(in) :: y(:), x(:, :)
    integer :: info, lwork
    real(dp) :: beta(k)

jtilly

print_graph = function(g, type = NULL, filename = NULL, height = 6, width = 6, path = ifelse(type == "tikz", "tex", "pdf")) {
    
    print(g)
    
    if(!is.null(type)) {
        
        if(type == "tikz") {
            tikz(paste0(path, "/", filename), height = height, width = width)
            print(g)
            dev.off()

jtilly

alias time='/usr/bin/time -f "\nCPU: %Us\tReal: %es\tRAM: %MKB"'