Hi Oleg,

Thanks. My Python function for the first derivative is (happily) identical in form. Such a function does seem like something of general utility that should perhaps be present in scitbx. Or maybe situations that require this are actually quite rare.

The code from David Thomas (e.g. RTMATS in Mosflm source) is nice because it does explicit calculation of up to the fourth derivative, after which they repeat. Although the need for that must surely be rather unusual.

Cheers

-- David


On 28 June 2012 11:34, Oleg Dolomanov <oleg@olexsys.org> wrote:
Hi David,

I did these some times ago, you can find them here:

For rotation around a vector:
at line 76
For 3 Euler angles:
at line 107

Also these derivatives present here:
at lines 61 and 130

Cheers,

Oleg.
From: dgwaterman@gmail.com
Date: Thu, 28 Jun 2012 11:13:45 +0100
To: cctbxbb@phenix-online.org
Subject: [cctbxbb] angle derivative of rotation matrix


Hi,

I'm using the axis_and_angle_as_r3_rotation_matrix method of scitbx.matrix.col in some Python code. Now I need to calculate the first derivative of the rotation matrix with respect to the rotation angle. Is there a function to do that somewhere in scitbx (or elsewhere in cctbx)? I had a brief look around but couldn't find anything that looked hopeful.

If there isn't something there already I have some Fortran code from David Thomas that I could adapt to do this. Of course it wouldn't be very efficient in pure Python, and calculating the matrix and its derivative separately, but that's okay for now as the code is just for testing out some ideas.

By the way, links to http://skal.planet-d.net/demo/matrixfaq.htm in the documentation (e.g. here) are broken. I found a document with the same name here: http://vamos.sourceforge.net/matrixfaq.htm, but I have no idea if this has the same content as the original.

Cheers

-- David

_______________________________________________ cctbxbb mailing list cctbxbb@phenix-online.org http://phenix-online.org/mailman/listinfo/cctbxbb

_______________________________________________
cctbxbb mailing list
cctbxbb@phenix-online.org
http://phenix-online.org/mailman/listinfo/cctbxbb