Andrew Cooke | Contents | Latest | RSS | 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

Choochoo Training Diary

Last 100 entries

Surprise Paradox; [Books] Good Author List; [Computing] Efficient queries with grouping in Postgres; [Computing] Automatic Wake (Linux); [Computing] AWS CDK Aspects in Go; [Bike] Adidas Gravel Shoes; [Computing, Horror] Biological Chips; [Books] Weird Lit Recs; [Covid] Extended SIR Models; [Art] York-based Printmaker; [Physics] Quantum Transitions are not Instantaneous; [Computing] AI and Drum Machines; [Computing] Probabilities, Stopping Times, Martingales; bpftrace Intro Article; [Computing] Starlab Systems - Linux Laptops; [Computing] Extended Berkeley Packet Filter; [Green] Mainspring Linear Generator; Better Approach; Rummikub Solver; Chilean Poetry; Felicitations - Empowerment Grant; [Bike] Fixing Spyre Brakes (That Need Constant Adjustment); [Computing, Music] Raspberry Pi Media (Audio) Streamer; [Computing] Amazing Hack To Embed DSL In Python; [Bike] Ruta Del Condor (El Alfalfal); [Bike] Estimating Power On Climbs; [Computing] Applying Azure B2C Authentication To Function Apps; [Bike] Gearing On The Back Of An Envelope; [Computing] Okular and Postscript in OpenSuse; There's a fix!; [Computing] Fail2Ban on OpenSuse Leap 15.3 (NFTables); [Cycling, Computing] Power Calculation and Brakes; [Hardware, Computing] Amazing Pockit Computer; Bullying; How I Am - 3 Years Post Accident, 8+ Years With MS; [USA Politics] In America's Uncivil War Republicans Are The Aggressors; [Programming] Selenium and Python; Better Walking Data; [Bike] How Fast Before Walking More Efficient Than Cycling?; [COVID] Coronavirus And Cycling; [Programming] Docker on OpenSuse; Cadence v Speed; [Bike] Gearing For Real Cyclists; [Programming] React plotting - visx; [Programming] React Leaflet; AliExpress Independent Sellers; Applebaum - Twilight of Democracy; [Politics] Back + US Elections; [Programming,Exercise] Simple Timer Script; [News] 2019: The year revolt went global; [Politics] The world's most-surveilled cities; [Bike] Hope Freehub; [Restaurant] Mama Chau's (Chinese, Providencia); [Politics] Brexit Podcast; [Diary] Pneumonia; [Politics] Britain's Reichstag Fire moment; install cairo; [Programming] GCC Sanitizer Flags; [GPU, Programming] Per-Thread Program Counters; My Bike Accident - Looking Back One Year; [Python] Geographic heights are incredibly easy!; [Cooking] Cookie Recipe; Efficient, Simple, Directed Maximisation of Noisy Function; And for argparse; Bash Completion in Python; [Computing] Configuring Github Jekyll Locally; [Maths, Link] The Napkin Project; You can Masquerade in Firewalld; [Bike] Servicing Budget (Spring) Forks; [Crypto] CIA Internet Comms Failure; [Python] Cute Rate Limiting API; [Causality] Judea Pearl Lecture; [Security, Computing] Chinese Hardware Hack Of Supermicro Boards; SQLAlchemy Joined Table Inheritance and Delete Cascade; [Translation] The Club; [Computing] Super Potato Bruh; [Computing] Extending Jupyter; Further HRM Details; [Computing, Bike] Activities in ch2; [Books, Link] Modern Japanese Lit; What ended up there; [Link, Book] Logic Book; Update - Garmin Express / Connect; Garmin Forerunner 35 v 230; [Link, Politics, Internet] Government Trolls; [Link, Politics] Why identity politics benefits the right more than the left; SSH Forwarding; A Specification For Repeating Events; A Fight for the Soul of Science; [Science, Book, Link] Lost In Math; OpenSuse Leap 15 Network Fixes; Update; [Book] Galileo's Middle Finger; [Bike] Chinese Carbon Rims; [Bike] Servicing Shimano XT Front Hub HB-M8010; [Bike] Aliexpress Cycling Tops; [Computing] Change to ssh handling of multiple identities?; [Bike] Endura Hummvee Lite II; [Computing] Marble Based Logic; [Link, Politics] Sanity Check For Nuclear Launch; [Link, Science] Entropy and Life

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

Permute Fucntion (Scheme)

From: "andrew cooke" <andrew@...>

Date: Tue, 14 Aug 2007 22:53:01 -0400 (CLT)

":";exec snow -- "$0" "$@"

(package* permute/v1.0.0
  (provide: (define (permute lst))))

; generates a list of all permutations
; assumes all values in original list are unique (under equal?)
(define (permute lst)
  (case (length lst)
    ((0) '())
    ((1) (list lst))
    (else (flat-map (lambda (head)
                      (map (lambda (perm) (cons head perm))
                           (permute (except head lst))))
                    lst))))

(define (flat-map f lst) (flat-map-acc f lst '()))
(define (flat-map-acc f lst acc)
  (if (null? lst)
      acc
      (flat-map-acc f (cdr lst) (append acc (f (car lst))))))

; this does not attempt to keep any kind of ordering
(define (except x lst) (except-acc x lst '()))
(define (except-acc x lst acc)
  (if (null? lst)
      acc
      (let ((head (car lst))
            (tail (cdr lst)))
        (if (equal? x head)
            (append acc tail)
            (except-acc x tail (cons head acc))))))


(test*
 (define (dup x) (list x x))
 (expect* (equal? '() (flat-map dup '())))
 (expect* (equal? '(1 1) (flat-map dup '(1))))
 (expect* (equal? '(1 1 2 2) (flat-map dup '(1 2)))))

(test*
 (expect* (equal? '() (except 0 '())))
 (expect* (equal? '() (except 1 '(1))))
 (expect* (equal? '(2) (except 1 '(1 2))))
 (expect* (equal? '(3 2 1) (except 0 '(1 2 3))))
 (expect* (equal? '(2 3) (except 1 '(1 2 3))))
 (expect* (equal? '(1 3) (except 2 '(1 2 3))))
 (expect* (equal? '(2 1) (except 3 '(1 2 3)))))

(test*
 (expect* (equal? '() (permute '())))
 (expect* (equal? '((1)) (permute '(1))))
 (expect* (equal? '((1 2)(2 1)) (permute '(1 2))))
 (expect* (equal? '((1 2 3)(1 3 2)(2 1 3)(2 3 1)(3 2 1)(3 1 2))
                  (permute '(1 2 3)))))

Improved Permutation Function (Start of List Library)

From: "andrew cooke" <andrew@...>

Date: Wed, 15 Aug 2007 12:07:38 -0400 (CLT)

":";exec snow -- "$0" "$@"

; library of list-related functions

(package* ac-lists/v1.0.0
  (provide:
      ; abstract the typical patter used to recurse over lists
      (define-syntax acl-list-process
        (syntax-rules (null?)
          ((_ ((null? list) null-body)
              ((head tail) body))
           (if (null? list)
               null-body
               (let ((head (car list))
                     (tail (cdr list)))
                 body)))))
    (define (acl-filter pred lst))
    (define (acl-filter-acc pred lst acc))
    (define (acl-except x lst))
    (define (acl-flat-map fun lst))
    (define (acl-flat-map-acc fun lst acc))
    (define (acl-permute lst))
    ))

; remove all values that match predicate from a list
(define (acl-filter pred lst) (reverse (acl-filter-acc pred lst '())))
(define (acl-filter-acc pred lst acc)
  (acl-list-process
   ((null? lst) acc)
   ((head tail) (if (pred head)
                    (acl-filter-acc pred tail acc)
                    (acl-filter-acc pred tail (cons head acc))))))

(test*
 (define (one? x) (= 1 x))
 (expect* (equal? '() (acl-filter one? '())))
 (expect* (equal? '() (acl-filter one? '(1))))
 (expect* (equal? '(2) (acl-filter one? '(2))))
 (expect* (equal? '(2 3) (acl-filter one? '(2 1 3)))))

; exclude instances of a particular value
(define (acl-except x lst) (acl-filter (lambda (y) (equal? x y)) lst))

(test*
 (expect* (equal? '() (acl-except 0 '())))
 (expect* (equal? '() (acl-except 1 '(1))))
 (expect* (equal? '(2) (acl-except 1 '(1 2))))
 (expect* (equal? '(1 2 3) (acl-except 0 '(1 2 3)))))

; map with append instead of cons
(define (acl-flat-map fun lst) (reverse (acl-flat-map-acc fun lst '())))
(define (acl-flat-map-acc fun lst acc)
  (acl-list-process
   ((null? lst) acc)
   ((head tail) (acl-flat-map-acc fun tail (append (fun head) acc)))))

(test*
 (define (dup x) (list x x))
 (expect* (equal? '() (acl-flat-map dup '())))
 (expect* (equal? '(1 1) (acl-flat-map dup '(1))))
 (expect* (equal? '(1 1 2 2) (acl-flat-map dup '(1 2)))))

; a list of all permutations
; assumes all values in original list are unique (under equal?)
(define (acl-permute lst)
  (define (perm-rest x)
    (define (cons-x rest) (cons x rest))
    (let ((rest (acl-except x lst)))
      (map cons-x (acl-permute rest))))
  (acl-list-process
   ((null? lst) '())
   ((head tail) (if (null? tail)
                    (list lst)
                    (acl-flat-map perm-rest lst)))))

(test*
 (expect* (equal? '() (acl-permute '())))
 (expect* (equal? '((1)) (acl-permute '(1))))
 (expect* (equal? '((1 2)(2 1)) (acl-permute '(1 2))))
 (expect* (equal? '((1 3 2)(1 2 3)(2 3 1)(2 1 3)(3 2 1)(3 1 2))
                  (acl-permute '(1 2 3)))))

Comment on this post