In my newest project I need to use Fortran routines and variables in C. From the source code point of view this requires only to declare the needed Fortran entities in your C code whilst taking the “right” types of the variables and the name mangling of the Fortran compiler into account.

Here is an example. Consider you have the following Fortran code:

module X
  real, dimension(0:1) :: v
 
  subroutine Foo(A, b, c, &
    & D)
    integer, intent(in) :: a, &
      & B
    real, intent(out) :: C, d
    ! <Foo's body>
  end subroutine
end module X

To use the array v and the subroutine Foo() in C, you need to add (something like) the following declarations to your C source code:

extern float x_mp_v_[2];
 
void x_mp_foo_(int* a, int* b, float* c, float* d);

If only a handful of Fortran routines are needed to be called from C, declaring them by hand may be feasible. This is however error-prone and becomes impracticable for more than a handful of routines. Especially if the routine signatures change frequently.

So what I’d like to have is a tool that generates a C header file containing all declarations of Fortran entities from a given Fortran source file. Or if such tool does not exist, a tool that extracts variable declarations and routine signatures from a Fortran source file into a format that is easier to parse than Fortran itself would also be helpful. Any suggestions?

This script is a very simple (and stupid) tool to get or set data in SLHA (SUSY Les Houches Accord) files. It can output complete blocks or individual data lines and replace existing data lines with user given input. It is stupid because it does not check after replacing data lines if the resultant file conforms to the SLHA. So care must be taken when using the “set” action.

But as so often, it is best described by examples. The basic things you can do with slha-io.sh are (MSSM.in is the SLHA file in these examples):

$ slha-io.sh MSSM.in get modsel
Block MODSEL            # Select model
 1    0                 # MSSM
 4    1                 # R-parity is violated
Block SMINPUTS          # Standard Model inputs

$ slha-io.sh MSSM.in get sminputs 6
 6   1.729000E+02       # m_top(pole)

$ slha-io.sh MSSM.in get sminputs 6 | awk '{ print $2 }'
1.729000E+02

$ slha-io.sh MSSM.in set sminputs 6 \
  " 6   1.730000E+02       # m_top(pole)"
$ slha-io.sh MSSM.in get sminputs 6
 6   1.730000E+02       # m_top(pole)

Here is another example of using slha-io.sh to create a batch of files where a parameter is constantly increased from 70 to 80 in steps of 0.5:

for NUM in `seq -f "%.8E" 70 0.5 80`; do
    FILE = "m0_$NUM.in"
    cp MSSM.in "$FILE"
    slha-io.sh "$FILE" set minpar 1 " 1  $NUM  # m0"
done

Maybe this script is useful for one or two!

By the way, more sophisticated tools for SLHA I/O are SLHALib and SLHAio.

To enable syntax highlighting in Vim for SLHA download the two files slha.vim and scripts.vim. Copy the slha.vim into ~/.vim/syntax/ and copy scripts.vim into ~/.vim/ or if you have already a scripts.vim file append the content of my file to it. Without further ado Vim should now highlight SLHA files. Here is the compulsory screenshot:

If you have further suggestions or found an error in the scripts, please don’t hesitate to contact me. I’ll try to keep the scripts up to date here.

Update 08.08.2010: The syntax file now includes the block names of the Flavour Les Houches Accord (FLHA, arXiv:1008.0762 [hep-ph]).

Geldautomat im abgesicherten Modus starten

Übrigens: Gratulation an das Debian-Projekt zur Veröffentlichung von Debian 5.0 alias Lenny.

Der Beweis dieser Behauptung ist einfach, wenn man als Kettenbruch entwickelt und dabei benutzt, dass und sind. Die erste Relation gilt offensichtlich, da das Argument der Gaußklammer immer kleiner gleich und größer als 1 ist. Für die zweite Relation betrachtet man die Ungleichungen und und sieht, dass beide zu Widersprüchen führen, die obige Relation also korrekt sein muss.

Behauptung: Für alle natürliche Zahlen gilt

und

Beispiele:

This is based on the following equation:

Which means that every natural number except powers of 2 can be written as a product of an odd natural number and a natural number . This can be easily understood. Just take the prime factors of an arbitrary natural number. If one of these prime factors is odd, the number can be written as a product of this odd prime factor and the product of the remaining prime factors (which can be even or odd). If there is no odd prime factor, all prime factors must be even and since 2 is the only even prime number the original number must be a power of 2.

I’m wondering if my solution is correct and how other people have written this statement in TNT.

Insgesamt muss ich aber sagen, dass ich froh bin, dass diese Wahl endlich zu Ende gegangen ist. Diese andauernden Nichtnachrichten über Sarah Palins Kopfsprünge in die großen Teiche voller Fett (die von belustigend bis beängstigend das volle emotionale Spektrum abgedeckt haben), oder die hundertste Umfrage, die Kandidaten X oder Y in Staat A oder B wieder vorne sah, waren zum Ende hin doch dann eher nervend. Trotzdem hat diese Wahl einen fantastischen Einblick in die Mediendemokratie der USA erlaubt. Teilweise sah es so aus, als ob Aussagen und Positionen der Kandidaten weniger wichtig waren, als die Summe der Werbespots mit denen Swing States überflutet wurden. Kann man sich in Amerika mit genügend Geld die Präsidentschaft erkaufen? Dieser Wahlkampf hat auch gezeigt, was ich vor einer Bundestagswahl in Deutschland nicht miterleben möchte. Kein Kanzlerkandidat sollte kurz vor de Wahl bei Schmidt & Pocher vorbeischauen, und es sollte auch egal sein in welchen Outfits unsere Politiker auftreten (ich hoffe ja auch mal, dass sowas nicht von Parteigeldern bezahlt werden darf).

Ich bin aber nicht nur froh darüber, dass diese Wahl nun vorbei ist, sondern auch über deren Ausgang. Zuersteinmal ist es hoch anzuerkennen, dass die USA einen Schwarzen zu ihrem Präsidenten gewählt haben. Zwar ist Obama kein Nachfahre von afrikanischstämmigen Sklaven, trotzdem wäre ein zukünftiger afroamerikanische Präsidentschaftskandidat kein Novum mehr, dass Ängste schüren würde. Ähnliches gilt nach Angela Merkel wahrscheinlich auch für zukünftige Kanzlerkandidatinnen in der Bundesrepublik. Aber wir sind noch meilenweit davon entfernt z.B. einen türkischstämmigen Politiker zum Kanzler zu wählen. Andererseits hat diese Wahl vielen jungen Menschen in Europa, die sich für Politik zu interessieren begannen als George W. zum Präsidenten (wieder-) gewählt wurde, ein Stück Hoffnung in die USA wiedergegeben und ihr Ansehen teilweise restauriert. Man kann sehr gespannt sein auf die nächsten Jahre, allerdings sollte man auch nicht zu viel erwarten. Nach Bush kann es aber nur besser werden, oder?

Bei Teil 4/11 wird es grausam. Auf die Frage: What is true and what is false about causes of climate change? antwortete Palin: [...] I don’t want to argue about the causes. What I want to argue about is how we are going to get there to positively affect the impacts [of climate change].

Mit so einer Einstellung bekommt man bestimmt auf die Idee die heiße Herdplatte mit Eiswürfeln zu kühlen, damit man sich nicht mehr die Finger verbrennt, anstatt den Herd einfach auszuschalten.

bts status file:- `bts select source:boinc archive:both` \
  | egrep "^originator" | cut -f2- | sort | uniq -ci | sort -n \
  | tail -n5

 2 Riccardo Stagni <unriccio@email.it>
 3 Martin Michlmayr <tbm@cyrius.com>
 6 Robert Millan <rmh@aybabtu.com>
 6 Thibaut VARENE <varenet@debian.org>
10 Sandro Tosi <matrixhasu@gmail.com>

Sweet – Thanks Debbugs and Devscripts maintainers!
And congratulations Sandro (or thank you – whatever you find appropriate)!

Update:
Lucas Nussbaum recently demonstrated how the Ultimate Debian Database (UDD) can be used to assemble cool stats about Debian bugs. This brought me to translate the above shell code into an SQL query for UDD. To retrieve boinc‘s top five bug reporter from UDD, one would use a query like this:

udd=> select submitter, count(*) from (select * from bugs
  where source = 'boinc' union select * from archived_bugs
  where source = 'boinc') as all_boinc_bugs group by submitter
  order by count desc limit 5;
              submitter               | count
--------------------------------------+-------
 Sandro Tosi <matrixhasu@gmail.com>   |     8
 Thibaut VARENE <varenet@debian.org>  |     6
 Robert Millan <rmh@aybabtu.com>      |     6
 Martin Michlmayr <tbm@cyrius.com>    |     3
 Michael Ablassmeier <abi@grinser.de> |     2
(5 rows)

The small difference between the two lists is because UDD does not handle bugs that are assigned to multiple packages correctly yet.

The actual conversion to Git was done with git-svn. I cloned the Subversion repositories with git-svn‘s --prefix=svn-import/ and --stdlayout options so that the trunk and all branches and tags were imported as remote branches with svn-import/ prepended. Then I created local branches of the previous Subversion branches and removed the remote remnants, for example:

git checkout -b debian/backports/etch svn-import/branches/etch-backports
git branch -d -r svn-import/branches/etch-backports

The Subversion tags are imported as svn-import/tags/<version> where <version> is the Debian version number. I don’t need branches for every version of my packages but I wanted to convert these tags to real Git tags but without losing the actual commit dates and messages. To achieve this I wrote this small script:

#!/bin/bash

for branch in `git branch -r`; do
    if [ `echo $branch | egrep "svn-import/tags/.+$"` ]; then
        version=`basename $branch`
        subject=`git log -1 --pretty=format:"%s" $branch`
        export GIT_COMMITTER_DATE=`git log -1 --pretty=format:"%ci" $branch`

        echo "Tag $version [Y/n]?"
        read yesno
        if [ -z $yesno ] || [ $yesno = "Y" ]; then
            git tag -s -f -m "$subject" "debian/$version" "$branch^"
            git branch -d -r $branch
        fi
    fi
done

For each remote branch that contains svn-import/tags/ it gets the version number, the commit message of the tag and exports GIT_COMMITTER_DATE (man git-tag) with the value of the date the version was tagged in Subversion and then asks you if it should tag the parent of the tag branch with the original commit date and message and remove the then unnecessary tag branch. The “safety measure” is there because people sometimes commit directly to Subversion tags (ugh!) and then "$branch^" would not be the commit you want to tag. To decide which Subversion tags are safe to tag with this script, gitk --all can be of great help.

P.S.: This is my first post on Planet Debian, so hello Planet!

Today the voting period for the Debian Project Leader Election 2008 has started, and here is the ballot I've just submitted:

-----BEGIN PGP SIGNED MESSAGE-----
Hash: SHA1

- - - -=-=-=-=-=- Don't Delete Anything Between These Lines =-=-=-=-=-=-=-=-
d81e16a2-03b6-4340-84f2-51de89b8185e
[ 2 ] Choice 1: Steve McIntyre
[ 1 ] Choice 2: Raphael Hertzog
[ 3 ] Choice 3: Marc Brockschmidt
[ 4 ] Choice 4: None Of The Above
- - - -=-=-=-=-=- Don't Delete Anything Between These Lines =-=-=-=-=-=-=-=-
-----BEGIN PGP SIGNATURE-----
Version: GnuPG v1.4.6 (GNU/Linux)

iD8DBQFH73cBft6HNdxCZCkRAmIgAJ49VGjuI6g+tIcgXhqMHXqE/kTNjQCcDqy5
oQSqQHJ5YCYtYE7y05mNdms=
=BwB2
-----END PGP SIGNATURE-----

This is the first election where I can actively participate since I became a Debian Developer on December 4th, 2007, and therefore I'm very interested in this year's DPL candidates and the election in general. I'm however a little bit disappointed because in my perception the candidates and their goals are less controversial than in the years before. It is also interesting to note that since 2004 this is the first election that has fewer candidates than the election in the preceding year (2008: #3, 2007: #8, 2006: #7, 2005: #6, 2004: #3).

Artikel 27
[ Handelsflotte ]

Alle deutschen Kauffahrteischiffe bilden eine einheitliche Handelsflotte.

Ich kann mir nicht helfen, aber ich stelle mir hier große Segelschiffe à la Gorch Fock vor, mit deutscher Flagge und Kanonen bewaffnet, und jederzeit bereit, um auf Kaperfahrt zu gehen... :-)

Anscheinend ist meine Vorstellung teilweise richtig. Ich frage mich, welche Relevanz dieser Artikel wohl hat.

Sehr bedauenswert, dass die Abgeordnete Dagmar Metzger von der eigenen Partei aufgrund ihrer Aufrichtigkeit aus dem hessischen Landtag gemobbt wird. Wer solche Parteifreunde hat, dem wünscht man nicht seine schlimmsten Feinde an den Hals; aber genau das wird sich Andrea Ypsilanti auch gedacht haben, als Dagmar Metzger am Freitag vor die Presse trat, um zu verkünden, dass sie ihre Landesvorsitzende nicht zur Ministerpräsidentin wählen werde.

Ich bin mir nicht sicher wie ich gestern darauf gekommen bin, aber ich wollte mir noch einmal die Rede von Werner Schulz (Bündnis 90/Die Grünen) vor dem 15. Deutschen Bundestag bei der Aussprache zur Vertrauensfrage anschauen, und tatsächlich habe ich Auszüge daraus in der PHOENIX Bibliothek gefunden. Eine komplette Aufzeichnung der Rede konnte ich bisher noch nicht finden, aber in Textform liegt sie z.B. hier vor.

Während ich mir den Zusammenschnitt der Aussprache nochmal anschaute ist mir ein lustiger Versprecher mit prophezeihendem Charakter bei der Rede von Angela Merkel aufgefallen. Sie sagte: »Rot-Grün kann unser Land nicht mehr regieren, die PDS darf unser Land nicht regieren, CDU und CSU gemeinsam mit der SPD...«. Natürlich erntete sie mit diesem Versprecher schallendes Gelächter aus fast allen Fraktionen im Hohen Haus, doch hätten die Genossen geahnt wie recht Angela Merkel damit hatte, hätten sie sich wahrscheinlich das Lachen verkniffen.

Cause: “DAM Approval: Approved on 2007-10-21 by Joerg Jaspert (joerg)”. And this happened only 23 days after Front Desk checked the completeness of my report. I'm somehow impressed, thought this would take longer. Thanks to Joerg Jaspert!

Now I'm anxiously awaiting my account creation, let's see how long this takes...

While being a Debian New Maintainer applicant, I decided to help the QA team by working on orphaned packages once I'm a Debian Developer. Currently there are 321 packages which are orphaned (O) or which a maintainer intents to adopt (ITA), see the packages overview for the Debian QA Group. Okay, so I've the choice to work on 321 different packages now, but which one should I choose? A tool which can help with this decision is wnpp-alert, it checks for installed packages that are up for adoption, orphaned or for which help is requested (RFH). It narrows down my selection to four orphaned packages, great! Now I want to sort these packages by their popularity (What's this?). For this task I've written a little script which uses the excellent grep-dctrl and ept-cache tools:

#!/bin/sh

MAINT=`(test -z "$1" && echo packages@qa.debian.org) || echo "$1"`

for x in `grep-status -F Maintainer "$MAINT" -n -s Package`; do
    ept-cache show $x | grep-dctrl '' -n -s Popcon,Package | paste -sd '\t';
done | sort -nr

It searches for installed packages that are maintained by the Debian QA Group and sorts them by their popcon score. If an argument is given to this script it will use it as search pattern for the Maintainer field instead of packages@qa.debian.org.

Note: You need to enable the popcon data provider for ept-cache for this script to work. See /usr/share/doc/ept-cache/README for details.

In the winter term 2006/07 I attended the lecture "Relativistic Quantum Field Theory" hold by Prof. Rückl at the University of Würzburg. Since I was (and still am) interested in this topic and wanted to learn more about the mathematical methods used in QFTs I decided to do a mini research project about the decay of the muon in the electroweak SM. Recently I finished the paper of my project which is unfortunately only available in German. You can download it here: Myon-Zerfall.pdf, and for your convenience I cite the abstract here:

Das Myon ist ein Elementarteilchen der zweiten Generation im Sinne des Standardmodells, welches identisch zum Elektron die Ladung −e trägt, einen Spin von 1/2 hat und bis auf die starke Wechselwirkung allen anderen Wechselwirkungen unterliegt. Es unterscheidet sich vom Elektron durch eine etwa 200 mal größere Masse und die Tatsache, dass seine Lebensdauer endlich ist. Der dominierende Zerfallskanal des Myons ist der Zerfall in ein Elektron, ein Myon-Neutrino und ein Elektron-Antineutrino, der im elektroschwachen Standardmodell durch ein W-Boson vermittelt wird. Diese Arbeit ist in drei Kapitel unterteilt. In dem ersten Kapitel wird nach Einführung der Feynman-Regeln der schwachen Wechselwirkung aus dem Feynman-Diagramm des Myon-Zerfalls die Feynman-Amplitude bestimmt und deren spin-gemitteltes Betragsquadrat im Ruhesystem des Myons berechnet. Um damit zur Zerfallsbreite zu gelangen wird im zweiten Kapitel der Drei-Teilchen-Phasenraum betrachtet und das entsprechende Phasenraumelement ausgewertet. Im dritten Kapitel werden dann die Ergebnisse aus den vorangegangenen Kapiteln kombiniert und das Energiespektrum der entstehenden Elektronen sowie die Lebensdauer des Myons berechnet und diese mit dem experimentellen Wert verglichen. Zum Schluß dieser Arbeit werden noch Korrekturen zu den Ergebnissen diskutiert und es wird auf die Bedeutung des Myon-Zerfalls eingegangen.