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