Skip to content

Instantly share code, notes, and snippets.

@jstults
jstults / falknerskan.mac
Created January 29, 2012 20:50
Maxima definition of Falkner-Skan governing equations
/* Falkner-Skan ODE
Zhang, J., Chen, B., 'An iterative method for solving the
Falkner-Skan equation'
Nachtsheim, P., Swigert, P., 'Satisfaction of Asymptotic Boundary
Conditions in Numerical Solution of Systems of Nonlinear Equations
of Boundary-Layer Type,' NASA TN D-3004, Lewis Research Center,
Cleveland, Oct, 1965.
@jstults
jstults / varcovid2.py
Last active December 2, 2020 10:14
Fit a vector autoregression moving average model to COVID-19 data
#
# filename: varcovid.py
# author: josh stults
# date: 30 Nov 2020
#
# fit some vector auto regression models to deaths and cases
#
# data from:
# https://covid.ourworldindata.org/data/owid-covid-data.csv
#
@jstults
jstults / varcovid.py
Created November 20, 2020 19:20
Vector autoregression model fit to COVID-19 data
#
# filename: varcovid.py
# author: josh stults
# date: 20 Nov 2020
#
# fit some autoregressive models to deaths and cases
#
# data from:
# https://covid.ourworldindata.org/data/owid-covid-data.csv
#
//CUSTOMIZER VARIABLES
// Overall radius of the sphere
sphereRadius = 30;
// Radius of each hole at the surface of the sphere
holeTopRadius = 5.5;
// Radius of each hole at the bottom of the hole
holeBottomRadius = 5;
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
@jstults
jstults / solution_30.stl
Created January 24, 2019 11:30
inverse csg of openlsto example 3d cantilever
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
@jstults
jstults / ransac.conf
Created January 24, 2019 11:32
configuration file for InverseCSG RANSAC algorithm
target_sample_num: 50000
cluster_epsilon: 0.1
epsilon: 0.1
min_points: 500
normal_threshold: 0.8
probability: 0.5
max_iter: 100
max_diam: 0.1
gc_max_iter: 10
base: 5.0
@jstults
jstults / topopt2CSG.py
Created January 24, 2019 11:34
python script to run InverseCSG on an example part from OpenLSTO based on the run_tests script that ships with InverseCSG
import os
import sys
# Usage: python3 topopt2CSG.py <build_dir> <test_name>
# Example:
# python3 topopt2CSG.py ../build openlsto
#
if len(sys.argv) < 3:
print('Error: please specify the test case name.')
@jstults
jstults / pure-bending.tpd
Last active December 16, 2017 22:55
Pure bending ToPy problem definition
[ToPy Problem Definition File v2007]
# Author: Josh Stults
# pure bending topology optimization on a 640x720 domain
PROB_TYPE : comp
PROB_NAME : pure-bend
ETA : exp
DOF_PN : 2
VOL_FRAC : 0.3
@jstults
jstults / solidvoid.py
Created December 7, 2013 20:32
script to process 1-bit pngs that define solid and void regions for ToPy topology optimization problems
# read in black & white pngs and create PASV_ELEM and ACTV_ELEM lists
# for topy input decks
#
# void.png: black pixels indicate regions that should be void
# solid.png: black pixels indicate regions that should be solid
#
from PIL import Image
solidim = Image.open("solid.png")
voidim = Image.open("void.png")