ntessore

#ifndef EIGEN_SVDLIBCSUPPORT_H
#define EIGEN_SVDLIBCSUPPORT_H

#include <Eigen/Sparse>

namespace Eigen
{
	namespace svdlib_h
	{
		extern "C"

ntessore

#!/bin/bash

##############################################################################
# One command to edit files from the command line. Helps your muscle memory. #
#                                                                            #
# Do you instinctively write `mate` or `atom` to edit files, only to get     #
# told `command not found` because you forgot you are in a SSH session?      #
# The `edit` script is a simple wrapper that launches the command specified  #
# through the `$EDITOR`, `$VISUAL` or `$EDIT_EDITOR` environment variable.   #
# In this way, you will always use the same command to start editing files.  #

ntessore

ntessore

#!/usr/bin/env python

####
# small command line tool to get a BibTeX bibliography from ADS bibcodes
# 
#     $ adsbibq.py '2015A&A...580A..79T' '2016MNRAS.463.3115T'
#     
#     Query Results from the ADS Database
#     
#     Retrieved 2 abstracts, starting with number 1.  Total number selected: 2.

ntessore

#!/bin/bash

####
# list "unneccessary" installed packages (i.e. not needed by any other package)
# 
# save:  $ curl -Lo /usr/local/bin/brew-unnecessary https://git.io/brew-unnecessary
# chmod: $ chmod +x /usr/local/bin/brew-unnecessary
# run:   $ brew unnecessary
####

ntessore

"""
Module to apply window function matrix to Cls

Compatible with matrices generated by https://github.com/ntessore/clwin.

Configuration:

Add an entry of the form `m_<spectrum> = <matrix.fits>` for each spectrum that
needs to be modified.

ntessore

export FC='/opt/local/bin/gfortran-mp-4.8'

set -e

mkdir lensed
pushd lensed
    git clone --depth=1 https://github.com/farhanferoz/MultiNest
    pushd MultiNest/MultiNest_v3.11_cmake/multinest
        mkdir build
        pushd build

ntessore

# download file if not already downloaded
def dl(url, file=None):
    from urllib.request import urlretrieve
    from urllib.parse import urlparse
    import os.path
    if file is None:
        file = os.path.basename(urlparse(url).path)
    if not os.path.exists(file):
        urlretrieve(url, file)
    return file

ntessore

import numpy as np
from scipy.fft import rfft, irfft
from scipy.special import loggamma


# constants
LN_2 = np.log(2)
LN_10 = np.log(10)

ntessore

// compute the logarithm of Bessel functions of large argument and order
//
// author: Nicolas Tessore <n.tessore@ucl.ac.uk>
// date: 6 Dec 2020
//

#include <math.h>


// compute the Bessel function ratio J_{nu-1}(x)/J_nu(x) from a continued