jmbr

	.text
	.align 4,0x90
	.globl __Z21matrix_vector_productPKdPd
__Z21matrix_vector_productPKdPd:
LFB18:
	vmovsd	8(%rdi), %xmm2
	vmovsd	(%rdi), %xmm1
	vmulsd	%xmm2, %xmm2, %xmm3
	vmovsd	16(%rdi), %xmm0
	vmovsd	32(%rdi), %xmm6

jmbr

#include <stdio.h>
#include <stddef.h>
#include <gsl/gsl_matrix.h>

int print_matrix(FILE *f, const gsl_matrix *m)
{
        int status, n = 0;

        for (size_t i = 0; i < m->size1; i++) {
                for (size_t j = 0; j < m->size2; j++) {

jmbr

#!/bin/sh

sbcl --eval "(ql:quickload '(:swank :cl-mpi))" \
     --eval "(defun toplevel () \
               (push #'mpi:mpi-finalize *exit-hooks*) \
               (mpi:mpi-init) \
               (let ((port (+ swank::default-server-port (mpi:mpi-comm-rank)))) \
                 (swank:create-server :port port :dont-close t)) \
               (sb-impl::toplevel-init))" \
     --eval "(save-lisp-and-die \"mpisbcl\" :executable t :toplevel #'toplevel)"

jmbr

This file goes in /usr/lib/gnome-flashback/

#! /bin/sh
exec gnome-session --session=gnome-flashback-stumpwm --disable-acceleration-check "$@"

jmbr

;;; StumpWM initialization file.

(ql:quickload "swank")

(in-package :stumpwm)

(setf *window-border-style* :none
      *resize-hides-windows* nil)

;;; Redirect output to ~/.stumpwm.d/stumpwm.log

jmbr

reset

load '~/gnuplot-colorbrewer/qualitative/Set2.plt' # Optional

set terminal push
set terminal epslatex size 3.22in, 2.255in
set output "fig-example.tex"

unset key
set border 3

jmbr

/*[ Radius of big circle ]*/
r0 = 150;

/*[ Radius of tube ]*/
r1 = 8; // r0/10.0;

/*[ Tolerance ]*/
rtol = 0.2;

module joint(r) {

jmbr

import numpy as np
from sklearn.decomposition import PCA
from scipy.spatial.distance import pdist


def get_points_from_square(num_points: int, ambient_dim: int) -> np.ndarray:
    """Get points from a grid of a linear manifold (a square in this case) endowed with the Euclidean metric.

    """
    aux = np.linspace(-1, 1, int(np.sqrt(num_points)))

jmbr

(defun switch-to-python-buffer ()
  "Switch to Python shell buffer quickly."
  (interactive)
  (let ((python-buffer (get-buffer "*Python*")))
    (if python-buffer
      (switch-to-buffer python-buffer)
      (run-python))))

(global-set-key (kbd "C-c z") 'switch-to-python-buffer)
(global-set-key (kbd "C-c C-z") 'switch-to-python-buffer)

jmbr

from typing import Tuple

import numpy as np
import scipy.spatial


DEFAULT_NUM_EIGENPAIRS: int = 10 + 1


def diffusion_maps(points: np.ndarray, epsilon2: float,