# C[omp]ute

Welcome to my blog, which was once a mailing list of the same name and is still generated by mail. Please reply via the "comment" links.

Always interested in offers/projects/new ideas. Eclectic experience in fields like: numerical computing; Python web; Java enterprise; functional languages; GPGPU; SQL databases; etc. Based in Santiago, Chile; telecommute worldwide. CV; email.

© 2006-2015 Andrew Cooke (site) / post authors (content).

## Chinese Writer Hu Fayun

From: andrew cooke <andrew@...>

Date: Fri, 2 Dec 2016 16:36:08 -0300

http://www.nybooks.com/daily/2016/11/28/china-inside-and-outside-system-hu-fayun/
http://www.raggedbanner.com/

Andrew

## Previous Entries

### Apricot Jam

From: andrew cooke <andrew@...>

Date: Mon, 21 Nov 2016 14:37:12 -0300

Fruit 2.5kg (after removing stones - bought 3kg)
Sugar (rubia) 1.25kg
Juice 1 lemon
5tsp (old) vanilla extract (at end)

Halved fruit and included skins.  Cooked until "superheated" and
catching n bottom pf pan.  Still appeared quite runny.  Initial taste
is quite tart (reduced sugar - previously using 2:3 instead of 1:2).

Early in cooking, foam overlfowed pan.  After cleaning up, adding some
(olive) oil seemed to stop that.

Andrew

### Excellent Article on USA Politics

From: andrew cooke <andrew@...>

Date: Mon, 21 Nov 2016 07:51:44 -0300

Missed this earlier.  Has a comment on police that is obvious, but
something I had bever seen or thought of before.

Andrew

(From https://ask.metafilter.com/302848/Help-me-build-a-post-Trump-reading-list)

### Why Trump Won

From: andrew cooke <andrew@...>

Date: Thu, 10 Nov 2016 12:49:30 -0300

Some good links here -

Andrew

### Doxygen + Latex on CentOS 6

From: andrew cooke <andrew@...>

Date: Tue, 1 Nov 2016 10:55:37 -0300

The following magic incantations were necessary to get Doxygen latex
output to convert to pdf using pdflatex on CentOS 6.8 (which does not
have luatex):

# increase by a factor of 10 the values of
# main_memory, pool_size and save_size
sudo emacs -nw /usr/share/texmf/web2c/texmf.cnf

# don't know if this is really needed
sudo texhash

# this ends with an error, but appears to work anyway(?!)
sudo fmtutil-sys --all

Together they address the "TeX capacity exceeded, sorry" errors.

Andrew

### SMASH - Solve 5 Biggest Problems in Physics

From: andrew cooke <andrew@...>

Date: Mon, 31 Oct 2016 15:07:29 -0300

http://www.sciencealert.com/this-new-hypothesis-claims-to-have-solved-5-of-the-biggest-problems-in-physics

Andrew

### Good article on racism, brexit, and social divides

From: andrew cooke <andrew@...>

Date: Thu, 27 Oct 2016 19:08:41 -0300

https://www.theguardian.com/commentisfree/2016/oct/25/white-working-class-brexiters-politicians-bigotry

In particular, the way this is framed seems to get round the "these
are just racists" response that you see on places like MeFi (which,
curiously, welcomed this -
https://www.metafilter.com/163116/The-Weaponisation-of-the-Working-Class)

Andrew

From: andrew cooke <andrew@...>

Date: Thu, 27 Oct 2016 16:13:06 -0300

Fuck.  Mark Riley just played "Way We Won't".

Andrew

### Consciousness From Max Entropy

From: andrew cooke <andrew@...>

Date: Thu, 20 Oct 2016 21:55:26 -0300

http://physicsworld.com/cws/article/news/2016/oct/18/consciousness-is-tied-to-entropy-say-researchers

https://arxiv.org/abs/1606.00821

Hmmm.

Andrew

### Democrats

From: andrew cooke <andrew@...>

Date: Sat, 15 Oct 2016 21:43:55 -0300

http://www.lse.ac.uk/newsAndMedia/videoAndAudio/channels/publicLecturesAndEvents/player.aspx?id=3606

Not a happy listen.

Andrew

### Harvard Will Fix Black Poverty

From: andrew cooke <andrew@...>

Date: Sat, 15 Oct 2016 21:33:29 -0300

https://www.instagram.com/p/BLhgSQZgZhx/

Andrew

### Modelling Bicycle Wheels

From: andrew cooke <andrew@...>

Date: Thu, 13 Oct 2016 10:50:46 -0300

I'm trying to do some simple (2D, simplified physics) modelling of
bicycle wheels.  Unfortunately, it's turning out to be harder than I
thought.

The immediate problem is that I am not able to find (efficiently and
reliably) the stable solution.  It seems that my equations for forces
and energy are correct, because occasionally a solution converges.
The difficulty is "numerical" in some sense.

I can think of three problems:

- Local forces

I calculate the force at each "spoke hole" as the sum of three
components - the spoke (in tension, connected to the fixed hub), and
the two rim segments (in compression, spanning to adjacent holes).
The simplified physics assumes a "hinge" at each hole.

This sounds fairly reasonable (in particular, I have a Jobst quote
somewhere supporting the hinge model), except that it assumes that
the rim segment endpoints are in equilibrium.  In fact there may be
net forces there (when not at equilibrium) which are "passed along"
and also affect the hole under examination.

I don't see any way to avoid this.  It's the ugly compromise one
gets when finding a static solution to a system that would go
through some dynamic evolution.  Ideally you would model the entire
dynamics.  Instead, you iterate the static minimisation.  But some
of the (GSL) minimisation routines I am using seem to handle this
poorly.

- Local movement

I describe the system as a collection of (x,y) coordinates of the
holes (spoke end points).  So a simple physical action like
contracting the rim inwards in response to spoke tension involves
altering all the variables.

I am not sure this is a problem for methods that use a matrix-based
solution (since the solution can move multiple points), but it seems
like it could seriously affect simplex minimisation.

Instead, I could use some kind of Fourier scheme where zeroth order
is a constant radial shift, etc etc.  The problem here is that it
seems to make calculating derivatives much more expensive (although
derivatives are not needed for simplex).

- Unstable collapse

The hinge model assumes that the rim segments always have an
internal (hub facing) angle less than 180 degrees.  If this is
exceeded then the wheel becomes unstable and collapses.

This may be a source of instabilities in the solution when an
additional force is applied.  It would be useful to have some kind
of minimisation that could easily detect this condition.

Given that I am currently stalled it seems like I have two options.  I
can implement the Fourier based approach (second above) or try to do
some kind of customized solution "by hand".

A hand-rolled approach could deal with item 1 above by using an
iterative alorithm (quasi-dynamic).  It could also handle the third
point if it used calculations tailored exactly to the configuration
used.

So I could calculate the approximate solution for each hole in
isolation as a pair of 1D roots (zero force perpendicular and parallel
to the chord from neighbouring holes).  This would be quick and allow
detection of the instability.  By moving the hole only a fraction of
the distance to the solution and then iterating the entire wheel I
would move towards a consistent solution.

Andrew