svn log | ruby -e "puts STDIN.read.split(/\n/).select{|l| l =~ /^r\d+/}.map{|l| l.split('|')[1].strip}.uniq.sort"

Who has any guess why I did need this?

The ruby script basically automates the upgrade instructions from the DokuWiki main page. So the following actions are performed when executing the script:

1. Making a backup into /tmp/dokuwiki_backup_#{timestamp} of every installation.
2. Downloading the dokuwiki release (passed in as a parameter).
3. Extracting the files and copying everything to the installations (execept for the content of the /data directory).
4. Creating missing folders in the /data directory, making the owner www-data:www-data and chmodding them to 664.
5. Deleting files from a list of file from older revisions.

Within the script you need to specifiy this snipped for setting your DokuWiki installations:

# Definition of existing installation
INSTALLATIONS = [
  '/path/to/docu/wiki/installation1',
  '/path/to/docu/wiki/installation2'
].freeze

Then you can call for a new release as follows:

/path/to/script/upgrade_dokuwiki.rb http://www.splitbrain.org/_media/projects/dokuwiki/dokuwiki-2009-12-25c.tgz

If you want to improve the script, feel free to fork me on Gist.

The procedures for sample size calculation are well published, but some of those work in a manner like this:

def calculate_sample_size
  n = 4 # 4 is a good starting point
          # a we have 2 arms minimum and a trial with propability involved need 2 probands at least
  loop do
    return n if fullfills_a_some_mathematical_condition(n)
    n += 1
  end
end

So we need to return the smallest integer fulfilling a certain condition and we need to iterate up to this one, testing against a condition each time.

The above works, but has a big problem: I’m doing it for at least 10 million times, usually around 100 million times. Caching is not an option here, as we are dealing with stocastic we don’t have deterministic parameters (those parameters are hidden in the code above).

But, going back to my computer science books, I found ‘binary search’ which I applied for this – of course with a Ruby closure:

module IterationHelper
  # This is the entry point for the iteration
  def iterate_up start = 0, step = 1000, &block
    max = start
    # First we need to find the a maximum, because in theory we have infinite integers
    # I define the step parameters, as for some sample size calculation, there is a way to guess those to
    #   some extend.
    # I could of course start with Fixnum::MAX but I in my case this is faster
    until(yield(max))
      max += step
    end
 
    # Now that we know the range where to look for the smallest match
    find_smallest_satisfier(max - step, max, &block)
  end
 
 
  private
  # This works nearly exactly as textbooks binary sort
  def find_smallest_satisfier min, max, &block
 
    # If min and max differ less than 2, we have found the smallest n
    if max - min < 2
      return min if yield(min)
      return max
    end
 
    # Get the value in the middle
    n = min + (max-min)/2
    # If yield(n) ==> value might be smaller
    if yield(n)
      return find_smallest_satisfier(min, n, &block)
      # Unless yield(n) ==> value might be bigger
    else
      return find_smallest_satisfier(n, max, &block)
    end
  end
end

This leeds to a much cleaner call like this:

include IterationHelper
def sample_size
  iterate_up {|n| fullfills_a_some_mathematical_condition(n) }
 end

But – i did say – I didn’t do this for cleaner calls. So I did a litte benchmark:

require 'lib/helpers/iteration_helper'
include IterationHelper
 
require 'benchmark'
include Benchmark
 
N = 10**3
 
def complex_calculation n
  # some dummy calculations
  10.times {Math.sqrt Kernel.rand}
  n >= @n_0
end
@targets = [25, 105, 234, 500, 765]
 
 
bm(7) do |x|
  x.report("normal") do
  @targets.each do |t|
    @n_0 = t
      N.times do
        m = 4
        loop do
           break if complex_calculation(m)
           m += 1
        end
      end
    end
  end
 
  x.report("binary") do
    @targets.each do |t|
        @n_0 = t
        N.times do
        iterate_up(4) {|n| complex_calculation(n)}
      end
    end
  end
end

The result was very impressing:

user system total real
normal 12.210000 0.030000 12.240000 ( 15.631090)
binary 0.520000 0.010000 0.530000 ( 0.778788)

You should – of course – expect the iteration do do better, when the actual satisfaction test needs less mathematics, but perform worse if the mathematics are more time consuming (which mine are).

First of all we have a method in the DSL:

simulate do
# Some other definition ...
arms do
treatment N([0, 0.1, 0.3], 1)
placebo N(0,1)
end
# Some more definitions ...
end

The intersting thing happens when you call the N function. If you just call it with simple numeric arguments, you get back a normal distributed sampler. If you call it with an array as one or more arguments, you get back the multiple normal distributed samples – constructed from the cartesian product of all params. This means, the treatment will get the Array as if you called:

[N(0,1), N(0.1,1), N(0.3,1)]

The idea behind this is, that if you specify those multiple parameters for the simulation, the simulation is run multiple times – each time with a sligly different set of definitions.

A direct implementation wouldn’t be to hard – at least we are doing ruby. But as this is not only related to the N method, but to a lot of other things as well, I wanted to extract this rather into a model so I could write the N method something like this:

module DistributionHelper
  extend MultipleParameters
  def N mean, variance
    Distribution::Gauss.new(mean,variance)
  end
  mutate_params :N
end

The MultipleParameters module is implemented following the path, Dave Thomas led in his 5th Episode of “The Ruby Object Model and Metaprogramming” for memorization. I’ll try to explain the module but first have the implementation:

module MultipleParameters
  def mutate_params(name)
    original_method = instance_method(name)
 
    define_method(name) do |*args|
      mutated_arguments = [[]]
      args.each do |arg|
        if arg.is_a?(Array)
          new_results = []
          mutated_arguments.each do |r|
            arg.each do |a|
              rm = r.clone
              rm < < a
              new_results << rm
            end
          end
          mutated_arguments = new_results
        else
          mutated_arguments.each { |r| r &lt;&lt; arg  }
        end
      end
 
      bound_method = original_method.bind(self)
      ret = []
      mutated_arguments.each do |mutated_args|
        ret << bound_method.call(*mutated_args)
      end
 
      (ret.size == 1) ? ret.first : ret
    end
  end
end

So this is a class-level module (i.e. you need to get access to class methods rather than object methods) so you extend your class/module with it rather than including it. After having done this, you can call the mutate_params method for any previously defined method. When called, your originally defined method will get boxed into an object (line 3) and it will be redefinded. The redefined method will create the cartesian product of the argumenes (i.e. ([0,1], [1,2]) yields to (0,1), (0,2), (1,1) and (1,2)) and than create an array of all the results from the calls. Before beeing able to call the originally defined method you have to bind the method to the current self (line 24), because the method will be called as an object method, and than passing all different parameter sets to it.

Finally you have an array of all return values which is returned. The last line is somewhat of convinience as it returns the object directly if you did not specify any multiple parameters.

Note: This approach does not work, if you original method needs an array to be passed. I will probably extends the approach in the future to work with ranges as well. ]]> http://jthoenes.bergischweb.de/2009/02/03/argument-fun-with-ruby/feed/ 0 http://jthoenes.bergischweb.de/2008/11/12/randi2-day-of-code/ http://jthoenes.bergischweb.de/2008/11/12/randi2-day-of-code/#comments Wed, 12 Nov 2008 13:36:49 +0000 Johannes Thönes http://jthoenes.bergischweb.de/?p=31 ich möchte euch an dieser Stelle kurz auf den RANDI2 Day of Code am 28.11 in der Hochschule Heilbronn hinweisen. Details dazu im Projektblog: http://randi2.hs-heilbronn.de/blog/index.php/2008/11/12/randi2-day-of-code

Nun trug es sich aber zu (;-)), dass heise.de sein Design geändert hat. Mich hat deren Design schon vorhher immer genervt, wenn ich aus meinem RSS-Reader auf die einzelnen Artikel bin. Aber jetzt war definitiv die Schmerzensgrenze überschritten. Um nicht alle Elemente nachstylen zu müssen (was ja auch sehr fragil ist, da die ids und classes ja jederzeit geändert werden können), habe ich mir überlegt, das ich erstmal die Styles der Printversion nehme. Das ging erstaunlich einfach:

// Erstmal JQuery ins Dokument einführen
// Der Code zum einbinden von JQuer ist ist von:
// http://www.joanpiedra.com/jquery/greasemonkey/ und ist unter MIT-Lizenz gestellt
// Kommentare sind angepasst.
var GM_JQ = document.createElement('script');
GM_JQ.src = 'http://jquery.com/src/jquery-latest.js';
GM_JQ.type = 'text/javascript';
document.getElementsByTagName('head')[0].appendChild(GM_JQ);
 
// Warten bis JQuery nachgeladen ist
function GM_wait() {
if(typeof unsafeWindow.jQuery == 'undefined') { window.setTimeout(GM_wait,100); }
else { $ = unsafeWindow.jQuery; letsJQuery(); }
}
GM_wait();
 
// Das normale Greasemonkey Skript, nur jetzt mit JQuery
function letsJQuery() {
// Lösche alle Stylesheets die nicht media=print sind
$('link[rel=stylesheet]').not('[media=print]').remove();
// Ändere das media von print auf screen
$('link[rel=stylesheet]').filter('[media=print]').attr('media', 'screen');
}

Ich habe das Skript speziell für heise.de noch etwas erweitert … wen das interessiert, der möge mich ansprechen.

Ich möchte an dieser Stelle – zumindest partiell – das Wahlgeheimnis aussetzen: Ich habe eine Stimme für den Senat Sandra Scheel gegeben. Und alle anderen empfehle ich, dass selbe zu tun (und ggf. auch darüber zu bloggen …).

Ich kenne Sandra vor allem von Theater – sie hat sich damals extrem schnell in die Gruppe eingefunden und auch Verantwortung übernommen. Sie war immer zuverlässig (was man beim besten Willen nicht von allen bei uns behaupten kann!) und hat das, was sie uns zugesagt hat, stets auch getan. Ich glaube, das sie die richtige ist, uns Studenten im Senat zu vertreten und bin mir sicher, dass sie für Anliegen ein offenes Ohr hat.

Insofern: VOTE FOR SANDRA!

P.S.: Bevor es jemand anmerkt. Sandra tritt in einem Team von vier Studenten auf – kann aber einzeln gewählt werden.

der eine oder andere von euch hat in den letzten Tagen vielleicht schon das eine oder andere meinem neusten Projekt gehört. Ab jetzt geht es also los:

Der MIdizini Podcast

In diesem Podcast möchte ich – zusammen mit Hauke Hund – ca. alle 4 Wochen Projekte innerhalb des MI Studiums beleuchten. Dazu werden wir mit Teilnehmern der verschiedenen Projekte Interviews führen. Den Anfang hat Herr Haag gemacht – natürlich mit dem CAMPUS Projekt. Ich hoffe das diese Angebot viele Studenten interessiert. Hauke und ich sind schon gespannt auf eure Rückmeldungen … also lasst und wissen, was euch gefällt und was nicht.

Wer mehr über die Idde und Hintergründe zum Podcast erfahren will, kann dies erfahren unter: http://midizini.bergischweb.de/blog/derpodcast/

Die Länge des Podcasts ist übrigens genau auf eine Zugfahrt zwischen HD und HN bzw. zurück zugeschnitten. Vorbild für den Podcast war – zumindest für mich – der exzellente Podcast SE-Radio (http://www.se-radio.net) den ich übrigens auch jedem nur empfehlen kann …

Aber, ich habe jetzt was gefunden was ich echt hilfreich finde: Empfohle Blogpost im Google Reader. Hier kann ich meine Empfehlungen zu Posts – die ich persönlich gut gefunden habe – als RSS-Feed weitergeben. Natürlich funktioniert das nur, wenn euch auch die Dinge gefallen die mir gefallen. Ist das der Fall, so könnt ihr den Feed (ab jetzt Rechts in meinem Blog zu finden) abonnieren. Wenn nicht … naja, gut stört mich auch nicht.

Was ich habe gerne hätte wären andere Menschen, die dieses  Feature (oder ähnliche Features, denn ich glaube nicht dass es das nur bei Google Reader gibt) nutzen und mir die Feedadresse mitteilen. Man kann ja nicht alle Feeds lesen und möchte nichts verpassen. Also: Her mit euren Empfehlungen.

Ruby 1.9

Matz und Koichi Sasada haben das neue Ruby 1.9 vorgestellt. Und einen Überblick darüber gegeben, was war, was ist und was noch kommen wird. Neben den technischen Details (die man übrigens sehr gut auch aus dem Google TechTalk Vortrag von Matz entnehmen kann) hat Matz vorallem herausgehoben, dass für ihn die Innovation der Sprache Ruby sehr wichtig ist. Er möchte nicht mit dem zufrieden sein was ist, sondern möchte weiter gehen und Ruby weiter entwickeln.
Weiterhin ist natürlich zu erwähnen: Ruby programming is fun …

JRuby

Für mich war JRuby die absolute Überraschung. Bisher hat man sich immer gefragt, was eigentlich JRuby soll – nun im Prinzip ist JRuby dazu da, Ruby Code zu schreiben und gleichzeitig Java-Klassen zu verwenden. Whant so see some?
(Nach Installation via apt-get install jruby1.0 und interaktiv, also jirb1.0)

require 'java'
include_class 'javax.swing.JFrame'
include_class 'javax.swing.JLabel'
include_class 'java.awt.event.ActionListener'
include_class 'javax.swing.JOptionPane'
 
frame = JFrame.new "Das JRuby Fenster"
frame.setSize 400,200
# Überlagerter Setter von
# frame.set_default_close_operation(JFrame::EXIT_ON_CLOSE)
frame.default_close_operation=JFrame::EXIT_ON_CLOSE
label = JLabel.new "Hello World"
frame.add label
 
frame.visible= true

Dinge wie obejct.method etc gehen natürlich auch, sodass man super neuer Java-Bibliotheken über jruby ausprobieren kann, selbst wenn man sie ‘nur’ für Java selbst verwenden will.
Übrigens wurde die 1.1 Release während der Konferenz bekannt gegeben.

Weitere Vorträge

Wurden auch gehalten. Alle waren sehr interessant und haben vorallem eines gezeigt: Mit Ruby kann man verdammt viel mehr machen, als bloß Rails. Und dass heißt ja nun wirklich etwas … Ich fands jedenfalls sehr cool (schon alleine wegen der T-Shirts )

Die Folien sollten in Kürze http://www.euruko2008.org zu finden sein.