Skip to content

Instantly share code, notes, and snippets.

View larsoner's full-sized avatar

Eric Larson larsoner

View GitHub Profile
======================================================================
ERROR: test_decomp_update.TestQRdelete_d.test_delete_last_p_row
----------------------------------------------------------------------
Traceback (most recent call last):
File "/usr/lib/python2.7/dist-packages/nose/case.py", line 197, in runTest
self.test(*self.arg)
File "/home/larsoner/.local/lib/python2.7/site-packages/scipy/linalg/tests/test_decomp_update.py", line 308, in test_delete_last_p_row
q1, r1 = qr_delete(q, r, 0, a.shape[0], 'row')
File "scipy/linalg/_decomp_update.pyx", line 1652, in scipy.linalg._decomp_update.qr_delete (scipy/linalg/_decomp_update.c:27659)
ValueError: Reorthogonalization Failed, unable to perform row deletion.
# -*- coding: utf-8 -*-
"""
An example where effecting jittering of triggers occurs when
downsampling before epoching.
"""
import numpy as np
import matplotlib.pyplot as plt
# 1 sec of data @ 1000 Hz
fs = 1000. # Hz
@larsoner
larsoner / test_sph.py
Last active October 22, 2015 11:04
Samu validation
from __future__ import print_function
import sys
from os import path as op
import numpy as np
from scipy.io import loadmat
from mne.io import read_info
from mne import pick_types, pick_info
from mne.preprocessing.maxwell import (_concatenate_sph_coils, _sh_negate,
_bases_real_to_complex, _deg_order_idx,
import numpy as np
import time
from mne.time_frequency import cwt_morlet
dur = 3.
sfreq = 1000.
x = np.empty((306, int(dur * sfreq)))
def test():
return cwt_morlet(x, sfreq, np.arange(5, 30))
# -*- coding: utf-8 -*-
# see:
# http://eeweb.poly.edu/iselesni/lecture_notes/least_squares/
# LeastSquares_SPdemos/deconvolution/html/deconv_demo.html
import numpy as np
from scipy import linalg, sparse
import matplotlib.pyplot as plt
from pyeparse.utils import pupil_kernel
# -*- coding: utf-8 -*-
import numpy as np
import time
n_choices = 18
def switch_fun(x):
if x == 0:
# -*- coding: utf-8 -*-
# Authors: Eric Larson <[email protected]>
#
# License: BSD (3-clause)
from __future__ import print_function
from copy import deepcopy
import time
@larsoner
larsoner / test_ave.py
Created April 5, 2016 15:20
Time averaged evoked
# -*- coding: utf-8 -*-
"""
Example for creating a time-averaged Evoked file
"""
import numpy as np
import mne
data_path = mne.datasets.sample.data_path()
fname = data_path + '/MEG/sample/sample_audvis-ave.fif'
evoked = mne.read_evokeds(fname, condition='Right Auditory')
@larsoner
larsoner / failures.txt
Created May 3, 2016 20:32
MKL failure log
larsoner@bunk:~$ python
Python 2.7.11+ (default, Apr 17 2016, 14:00:29)
[GCC 5.3.1 20160413] on linux2
Type "help", "copyright", "credits" or "license" for more information.
>>> import scipy
>>> scipy.test('full')
Running unit tests for scipy
NumPy version 1.12.0.dev0+1fc180b
NumPy relaxed strides checking option: True
NumPy is installed in /home/larsoner/.local/lib/python2.7/site-packages/numpy-1.12.0.dev0+1fc180b-py2.7-linux-x86_64.egg/numpy
# -*- coding: utf-8 -*-
import os.path as op
import numpy as np
import mne
def select_center_vertex_in_label(src, label):
""" Select vertex closest (in Euclidean distance) to the center of mass of
the given label.