Andrew Cooke | Contents | Latest | RSS | Twitter | Previous | Next

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.

Personal Projects

Lepl parser for Python.

Colorless Green.

Photography around Santiago.

SVG experiment.

Professional Portfolio

Calibration of seismometers.

Data access via web services.

Cache rewrite.

Extending OpenSSH.

C-ORM: docs, API.

Last 100 entries

[Link, Programming] MySQL UTF8 Broken; NCOA for $20; [Link, Books] Latin American Authors; [Link, Computing] Optimizatin Puzzle; [Link, Books, Politics] Orwell Prize; [Link] What the Hell Is Happening With Qatar?; [Link] Deep Learning + Virtual Tensor Machines; [Link] Scaled Composites: Largest Wingspan Ever; [Link] SCP Foundation; [Bike] Lessons From 2 Leading 2 Trailing; [Link] Veg Restaurants in Santiago; [Link] List of Contemporary Latin American Authors; [Bike] FTHR; [Link] Whoa - NSA Reduces Collection (of US Residents); [Link] Red Bull's Breitbart; [Link] Linux Threads; [Link] Punycode; [Link] Bull / Girl Statues on Wall Street; [Link] Beautiful Chair Video; Update: Lower Pressures; [Link] Neat Python Exceptions; [Link] Fix for Windows 10 to Avoid Ads; [Link] Attacks on ZRTP; [Link] UK Jazz Invasion; [Review] Cuba; [Link] Aricle on Gender Reversal of US Presidential Debate; {OpenSuse] Fix for Network Offline in Updater Applet; [Link] Parkinson's Related to Gut Flora; Farellones Bike Park; [Meta] Tags; Update: Second Ride; Schwalbe Thunder Burt 2.1 v Continental X-King 2.4; Mountain Biking in Santiago; Books on Ethics; Security Fail from Command Driven Interface; Everything Old is New Again; Interesting Take on Trump's Lies; Chutney v6; References on Entropy; Amusing "Alexa.." broadcast; The Shame of Chile's Education System; Playing mp4 gifs in Firefox on Opensuses Leap 42.2; Concurrency at Microsoft; Globalisation: Uk -> Chile; OpenSuse 42.2 and Synaptics Touch-Pads; Even; Cherry Jam; Lebanese Writer Amin Maalouf; C++ - it's the language of the future; Learning From Trump; Chinese Writer Hu Fayun; And; Apricot Jam; Also; Excellent Article on USA Politics; Oh Metafilter; Prejudice Against The Rurals; Also, Zizek; Trump; Why Trump Won; Doxygen + Latex on CentOS 6; SMASH - Solve 5 Biggest Problems in Physics; Good article on racism, brexit, and social divides; Grandaddy are back!; Consciousness From Max Entropy; Democrats; Harvard Will Fix Black Poverty; Modelling Bicycle Wheels; Amusing Polling Outlier; If Labour keeps telling working class people...; Populism and Choice; Books on Defeat; Enrique Ferrari - Argentine Author; Transcript of German Scientists on Learning of Hiroshima; Calvert Journal; Owen Jones on Twitter; Possible Japanese Authors; Complex American Literature; Chutney v5; Weird Componentized Virus; Interesting Argentinian Author - Antonio Di Benedetto; Useful Thread on MetaPhysics; RAND on fighting online anarchy (2001); NSA Hacked; Very Good LRB Article on Brexit; Nussbaum on Anger; Tasting; Apple + Kiwi Jam; Hit Me; Sudoku - CSP + Chaos; Recycling Electronics In Santiago; Vector Displays in OpenGL; And Anti-Aliased; OpenGL - Render via Intermediate Texture; And Garmin Connect; Using Garmin Forerunner 230 With Linux; (Beating Dead Horse) StackOverflow; Current State of Justice in China; Axiom of Determinacy; Ewww; Fee Chaos Book

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

Uniformly Random, Correlated Numbers in Matlab/Octave

From: andrew cooke <andrew@...>

Date: Sat, 23 Jul 2011 11:27:19 -0400

Someone just asked me how to implement the work at
http://acooke.org/random.pdf in Matlab.  Here's the code for Octave, which is
very similar:


n = 1000;  # the number of points
t = 0.1; # the 't' parameter from the paper

sq2 = sqrt(2);

# first, generate points in the u,v coordinates as in figure 4.
# u goes from 0 to sq2, v from -t*sq2 to t*sq2
uv = rand(n,2) * [sq2,0;0,2*t*sq2] - repmat([0,t*sq2],n,1);
# now rotate that by 45 degrees to xy
xy = uv * [1/sq2,1/sq2;-1/sq2,1/sq2];
# finally, "fold in" points outside the unit square
for i = 1:n
  #printf("%d x %f y %f\n", i, xy(i,1), xy(i,2));
  if (xy(i,1) < 0) xy(i,1) = -xy(i,1); endif
  if (xy(i,2) < 0) xy(i,2) = -xy(i,2); endif
  if (xy(i,1) > 1) xy(i,1) = 2-xy(i,1); endif
  if (xy(i,2) > 1) xy(i,2) = 2-xy(i,2); endif
endfor;

#xy
plot(xy(:,1),xy(:,2),'*')


If anyone with more experience than me know of a better way of folding the
points (ie without the explicit "for" loop), please let me know!

Andrew

Avoiding For Loops

From: Jonathan Dursi <ljdursi@...>

Date: Sat, 23 Jul 2011 12:26:57 -0400

I think this is pretty easily done with the vector-friendly min & max 
operations:

function xy = randminmax(n)
     t = 0.1; # the 't' parameter from the paper

     sq2 = sqrt(2);

     # first, generate points in the u,v coordinates as in figure 4.
     # u goes from 0 to sq2, v from -t*sq2 to t*sq2
     uv = rand(n,2) * [sq2,0;0,2*t*sq2] - repmat([0,t*sq2],n,1);
     # now rotate that by 45 degrees to xy
     xy = uv * [1/sq2,1/sq2;-1/sq2,1/sq2];

     # finally, "fold in" points outside the unit square
     xy = max(xy,-xy);
     xy = min(xy,2-xy);
endfunction


t1=now(); xy=randfor(100000); t2=now(); t2-t1
ans =  1.8261e-05
t1=now(); xy2=randminmax(100000); t2=now(); t2-t1
ans =  1.4284e-07

-- 
Jonathan Dursi | SciNet, Compute/Calcul Canada

Re: Avoiding For Loops

From: andrew cooke <andrew@...>

Date: Sat, 23 Jul 2011 13:50:57 -0400

Thanks!  I was looking for something like a map - never ocurred to me to look
for specific functions.  I'll forward that to the person who asked.

(And even when I read your suggestion, before I looked at the code, I couldn't
see how the upper bound would work.  I'm clearly not getting the matrix-based
approach at all... guess I better practice!)

Andrew

Comment on this post