Craig Versek cversek

auth: Craig Wm. Versek ([email protected]) dates: gist created 2016/06/09 based directly from code written circa 2008 while I was a PhD candidate at UMass Amherst

These are simple pygame based pendulum physics simulations that started with a simple idea: perform the numerical integration (the "physics") within an infinite generator loop, supplying the parameters to the graphics rendering code on demand. The thought was to separate the number crunching from the application at large and simplify its organization and portability. The code can be configured to use one of a few numerical integration methods by swapping generators with the same interface. pygame is the only dependency that needs to be installed.

phys_gen.py illustrates the basic concept with no bells & whistles (nor graphics).

pendulum.py is a pygame driven simple pendulum interactive simulation -- try clicking and dragging the bob. Can use Heun, Steormer_Verlet, or RK4 methods.

doubled_pedulum.py has two double pendula running in

	################################################################################
	# this is code to compute NDVI from Blue (Visible channel) and Red (Infrared)
	# such as pictures taken with "SuperBlue" camera
	from PIL import Image
	import numpy

	img = Image.open("IMG_0512.JPG")
	imgR, imgG, imgB = img.split() #get channels
	#convert to double precision floating point..is this overkill? probably, could try 'float32' or 'float16'
	arrR = numpy.asarray(imgR).astype('float64')

	from pylab import *
	cdict = {
	'red' : ((0.00, 0.00, 0.00),
	(0.20, 0.00, 0.00),
	(0.40, 0.00, 0.00),
	(0.60, 0.00, 0.00),
	(0.80, 1.00, 1.00),
	(1.00, 1.00, 1.00)),
	'green': ((0.00, 0.00, 0.00),
	(0.20, 0.00, 0.00),

	import numpy as np
	import matplotlib

	def make_cmap_guassianHSV( num_segs = 100, #number of segments
	bandwidth = 0.25,
	red_center = 1.00,
	green_center = 0.75,
	blue_center = 0.50,
	name = "gaussianHSV"
	):

	import time, socket



	doc = """
	<!DOCTYPE html>
	<html>
	<head> <title>ESP8266 Pins</title> </head>
	<body> <h1>ESP8266 Pins</h1>
	<table border="1"> <tr><th>Pin</th><th>Value</th></tr> s </table>

	#!/bin/bash -x
	#FIRST TIME SETUP
	#install esp-open-sdk, Espressif ESP-IDF, and micropython dependencies
	echo "Installing esp-open-sdk, Espressif ESP-IDF, and micropython dependencies..."
	sudo apt-get update
	sudo apt-get install -y build-essential git make unrar-free unzip \
	autoconf automake libtool libtool-bin gcc g++ gperf \
	flex bison texinfo gawk ncurses-dev libexpat-dev \
	python sed libreadline-dev libffi-dev pkg-config \
	help2man python-dev python-serial wget linux-image-extra-$(uname -r)

	import numpy as np

	def speedx(sound_array, factor):
	""" Multiplies the sound's speed by some `factor` """
	indices = np.round( np.arange(0, len(sound_array), factor) )
	indices = indices[indices < len(sound_array)].astype(int)
	return sound_array[ indices.astype(int) ]

	def stretch(sound_array, f, window_size, h):
	""" Stretches the sound by a factor `f` """

	from digitalio import DigitalInOut, Direction, Pull
	import audioio
	import board
	import array
	import time
	import math

	FREQUENCY = 440 # 440 Hz middle 'A'
	SAMPLERATE = 8000 # 8000 samples/second, recommended!

	# The MIT License (MIT)
	#
	# Copyright (c) 2017 Dan Halbert for Adafruit Industries
	# Copyright (c) 2017 Kattni Rembor, Tony DiCola for Adafruit Industries
	#
	# Permission is hereby granted, free of charge, to any person obtaining a copy
	# of this software and associated documentation files (the "Software"), to deal
	# in the Software without restriction, including without limitation the rights
	# to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
	# copies of the Software, and to permit persons to whom the Software is

	import numpy as np
	from scipy.misc import imsave

	PIXELS_X = 256
	PIXELS_Y = PIXELS_X

	COLOR_SOLID_WHITE = (255,255,255,255)
	COLOR_SOLID_BLACK = (0,0,0,255)
	COLOR_FULLY_TRANSPARENT_BLACK = (0,0,0,0)