Thanks, Bob
Aha, the mysteries of matrix multiplication  :-D

I circumvented my problem creating dummy atoms and applying the rotation 
to them. But anyway I'd like to do it this way.
Still, I cannot make it work

 rotaProt = compare( {*.CA/2.1 & polymer=8} {*.CA/1.1 & polymer=8} )

 show rotaProt

// rotaProt = 
[[-0.7424616,0.2323488,-0.6283032,132.39798],[0.14965674,0.9717479,0.1
8250754,65.91373],[0.6529577,0.04147502,-0.756258,41.888985],[0.0,0.0,0
.0,1.0]]

 p1 = {112.7159 -1.6228517 -76.52575}
 p2 = {67.10559 -0.5669861 -24.276031}
 draw li @p1 @p2  // the line is displayed

  p11 = rotaProt * p1

 show p1
// p1 = {112.7159 -1.6228517 -76.52575}

 show p11
// p11 = "NaN"

 p11 = @{rotaProt * p1}

 show p11
// p11 = ""

 p11 = p1 * rotaProt 

 show p11
// p11 = {17326.207 -249.45786 -11763.213}

What am I missing?


El 11 Nov 2015 a las 17:27, Robert Hanson escribió:
> That should be
> 
>   myRot * myPoint
> 
> not
> 
>  myPoint * myRot


------------------------------------------------------------------------------
_______________________________________________
Jmol-users mailing list
Jmol-users@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/jmol-users

Reply via email to