Skip to content
example scopeplot.png
#include "quadratic.h"
#include "utility/dspinst.h"
void AudioSynthWaveformRes::update(void)
{
audio_block_t *block;
uint32_t i, inc;
int32_t ph, modph;
int32_t y, tmp;
uint32_t quadrant;
#include "resonant.h"
#include "utility/dspinst.h"
void AudioSynthWaveformRes::update(void)
{
audio_block_t *block;
uint32_t i;
if (c or s)
{

Distortion is 95 dB THD or so, has fewer frequency components than linear-interpolated LUT, but overall THD level is similar.

Can only go up to fs/6, though, because signed a can only go to Q31 1.0, not 2.0.

Found in http://citeseerx.ist.psu.edu/viewdoc/summary?doi=10.1.1.84.1650 but didn't work until I combined with http://www.musicdsp.org/showArchiveComment.php?ArchiveID=10

This is even faster, but distorted at low levels: https://gist.github.com/endolith/14bbb3217f9f58248722

Crude timing results (Should have measured with many instances in parallel instead):

#!/usr/bin/env python
"""
Copyright 2012 endolith at gmail com
Copyright 2009 Ged Ridgway at gmail com
Translation and modification of
http://www.mathworks.com/matlabcentral/fileexchange/26026-bipolar-colormap
Based on Manja Lehmann's hand-crafted colormap for cortical visualisation
1   1.5708
2   1.1107
3   1.0472
4   1.0262
5   1.0166
6   1.0115
# -*- coding: utf-8 -*-
"""
Created on Thu Nov 13 21:36:42 2014
"""
from __future__ import division, print_function
import numpy as np
from numpy import zeros, abs, argmin, conj, argsort, concatenate, delete
MRI matlab hot.png

Install QPDF, and then edit this batch file with the correct path and create a shortcut to this batch file. Then you can Open With the shortcut file for quick access and it creates an unlocked version of a locked PDF in the same folder as the original.

# -*- coding: utf-8 -*-
"""
Created on Wed Jul 30 18:55:13 2014
Just intonation classes for music theory experiments in Python
"""
from __future__ import division, print_function
from fractions import gcd as _gcd
from fractions import Fraction
Something went wrong with that request. Please try again.