Skip to content

Instantly share code, notes, and snippets.

🎯
Focusing

João Faria j-faria

🎯
Focusing
Block or report user

Report or block j-faria

Hide content and notifications from this user.

Learn more about blocking users

Contact Support about this user’s behavior.

Learn more about reporting abuse

Report abuse
View GitHub Profile
View cl_arguments.py
import sys
def fun():
print 'fun'
def test():
print 'test'
if __name__ == "__main__":
if len(sys.argv) > 1 and sys.argv[1] == 'test':
@j-faria
j-faria / GLS
Last active Aug 29, 2015
Fortran implementation of the Generalized Lomb-Scargle (GLS) periodogram
View GLS
GLS periodogram
@j-faria
j-faria / pearsonr.f90
Created Apr 2, 2015
Calculate the Pearson correlation coefficient of two arrays
View pearsonr.f90
real(kind=8) function pearsonr(x, y) result(r)
! given two arrays x and y, this function computes their Pearson correlation coefficient r
implicit none
real(kind=8), dimension(:) :: x, y
real(kind=8), dimension(size(x)) :: xt, yt
real(kind=8) :: ax, ay, df, sxx, sxy, syy
integer :: n
if (size(x) /= size(y)) STOP 'Dimension mismatch in pearsonr'
@j-faria
j-faria / change_par_file.py
Last active Aug 29, 2015
Change parameter file, in place
View change_par_file.py
import fileinput
import os
# create a copy of the template parameter file
os.system('cp template_file new_file')
# for example, in the template you have
# par1 = 10
# par2 = 'some path'
# and you want to change to
View hello.py
def HelloWorld:
print 'Hello World, from Python'
@j-faria
j-faria / generate_from_cov_matrix.py
Created Oct 10, 2013
Generate data with a particular correlation matrix - resultados da primeira geek meeting
View generate_from_cov_matrix.py
#!/usr/bin/env python
# -*- coding: utf-8 -*-
# Generate data with a particular correlation matrix
#
# A useful fact is that if you have a random vector x with covariance matrix Σ,
# then the random vector Ax has mean AE(x) and covariance matrix Ω=AΣA^T.
# So, if you start with data that has mean zero, multiplying by A will not
# change that.
#
@j-faria
j-faria / fisher_yates.py
Last active Dec 25, 2015
Shuffle an array in place using the Knuth-Fisher-Yates algorithm.
View fisher_yates.py
#!/usr/bin/env python
# -*- coding: utf-8 -*-
#
# Shuffle an array in place using the Knuth-Fisher-Yates algorithm.
from random import randrange
def shuffle(x):
""" Shuffle x in place using the Knuth-Fisher-Yates algorithm """
for i in xrange(len(x)-1, 0, -1):
@j-faria
j-faria / jorge-noise-covmat.py
Created Nov 15, 2013
Calculate and show the covariance matrix for template noise.
View jorge-noise-covmat.py
from numpy import *
from pylab import *
from glob import glob
from scipy.stats import nanmean
NN = '088' # change this
fileglob = 'Residuals_NN' + NN + '_Index*'
@j-faria
j-faria / powernoise.py
Created Dec 14, 2013
Generate samples of power law noise.
View powernoise.py
from math import floor
from numpy import arange, exp, pi, flipud, append, conj, real
from numpy.fft import ifft
from numpy.random import randn
def powernoise(alpha, N, randpower=False, normalize=False):
"""
Generate samples of power law noise. The power spectrum
of the signal scales as f^(-alpha).
View configfile.py
# file configfile.py
# create a custom INI file reader and read a file called 'test.ini'
import ConfigParser
class iniReader(ConfigParser.ConfigParser):
def as_dict(self):
d = dict(self._sections)
for k in d:
d[k] = dict(self._defaults, **d[k])
You can’t perform that action at this time.