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
[email protected]
https://lists.sourceforge.net/lists/listinfo/jmol-users
