Skip to content

Instantly share code, notes, and snippets.

View mathurinm's full-sized avatar

mathurinm

View GitHub Profile
@mathurinm
mathurinm / sparse.py
Created July 13, 2023 11:17
scipy vs numpy for sparse
import os
os.environ["OMP_NUM_THREADS"] = "1"
from time import perf_counter
from sklearn.utils import check_random_state
import numpy as np
import matplotlib.pyplot as plt
from scipy import sparse
d = 1000
@mathurinm
mathurinm / norms_X_col_numba
Last active October 30, 2021 16:36
computing column norms of sparse matrix with numba
import time
import numpy as np
from numba import njit
from numpy.linalg import norm
from scipy.sparse.linalg import norm as snorm
from libsvmdata import fetch_libsvm
X, y = fetch_libsvm("finance")
@njit
from celer import Lasso
import matplotlib.pyplot as plt
from scipy.optimize import fmin_bfgs
from numpy.linalg import norm
import numpy as np
from celer.datasets import make_correlated_data
import seaborn as sns
c_list = sns.color_palette("colorblind")
@mathurinm
mathurinm / sgd_saga_svrg_gd.py
Last active November 23, 2020 16:41
Example svrg saga sgd gd
import numpy as np
from numba import njit
from numpy.linalg import norm
import matplotlib.pyplot as plt
from sklearn.linear_model import LinearRegression
from libsvmdata import fetch_libsvm
# from celer.datasets import make_correlated_data
import numpy as np
from numba import njit
import matplotlib.pyplot as plt
from scipy.linalg import toeplitz
from numpy.linalg import norm
from sklearn.utils import check_random_state
from sklearn.linear_model import LinearRegression
def data(n_samples, n_features, rho=0.5, seed=24):
@mathurinm
mathurinm / lbfgs_l1logistic.py
Created March 12, 2019 02:44 — forked from vene/lbfgs_l1logistic.py
Solving L1-regularized problems with l-bfgs-b
"""l-bfgs-b L1-Logistic Regression solver"""
# Author: Vlad Niculae <vlad@vene.ro>
# Suggested by Mathieu Blondel
from __future__ import division, print_function
import numpy as np
from scipy.optimize import fmin_l_bfgs_b
import matplotlib.pyplot as plt
from matplotlib import rc
import seaborn as sns
rc('font', **{'family': 'sans-serif',
'sans-serif': ['Computer Modern Roman']})
params = {'axes.labelsize': 12,
'font.size': 12,
'legend.fontsize': 12,
'xtick.labelsize': 10,
\begin{wrapfigure}{h!}{0.42\textwidth}
\centering
\includegraphics[width=0.42\textwidth]{prebuiltimages/KL_Ball.pdf}\\
\caption{\label{plot:q_set_example}Example of sets $\Qset_{\alpha}(\Prob_{Z})$ and $\Qset_{\alpha}(\Prob_{Z})$ with $\balpha < \alpha$.}
\end{wrapfigure}
@mathurinm
mathurinm / main.py
Created March 28, 2018 15:33
Running R function with rpy2
# obtained with priceless help from Olivier Grisel
# rpy2 is available via pip: pip install rpy2
import numpy as np
from rpy2 import robjects
import rpy2.robjects.packages as rpackages
from rpy2.robjects import numpy2ri
from rpy2.robjects import pandas2ri
# we use function ltsReg of package robustbase
if __name__ == "__main__":
@mathurinm
mathurinm / mne_mayavi.py
Created February 2, 2018 10:30
The two lines I have to execute in ipython before using mayavi
from pyface.qt import QtGui, QtCore
%matplotlib qt