topic Co-ordinate transformations in IntelĀ® Fortran Compiler
https://community.intel.com/t5/Intel-Fortran-Compiler/Co-ordinate-transformations/m-p/1106892#M128293
<P>I hate co-ordinate transformations.</P>
<P>Interestingly Harrison's old code and his book included an interesting coding example for co-ordinate transformations. In the book he published the 3 by 3 matrices for his sample problem. I had checked several of the matrix samples, but not the last one , turns out the code produces a slight variation on the matrix, the published answers agree with mine, but is that because the code is wrong or the book is wrong.</P>
<P>In moving from a 4 by 4 degree of freedom method to a 6 by 6 so I can match the shell elements to the beam elements I took Harrison's code, and then started looking at the final numbers -- I noticed the third beam has some numbers negative when the book says positive.</P>
<P>Turns out that in the book Harrison draws the beam going down the negative z axis - gamma = -90 , but the sample has gamma, but the coded example has gamma = -90</P>
<P>Ah the joys of coding.</P>
<P>John</P>Thu, 23 Jun 2016 19:08:23 GMTJohnNichols2016-06-23T19:08:23ZCo-ordinate transformations
https://community.intel.com/t5/Intel-Fortran-Compiler/Co-ordinate-transformations/m-p/1106892#M128293
<P>I hate co-ordinate transformations.</P>
<P>Interestingly Harrison's old code and his book included an interesting coding example for co-ordinate transformations. In the book he published the 3 by 3 matrices for his sample problem. I had checked several of the matrix samples, but not the last one , turns out the code produces a slight variation on the matrix, the published answers agree with mine, but is that because the code is wrong or the book is wrong.</P>
<P>In moving from a 4 by 4 degree of freedom method to a 6 by 6 so I can match the shell elements to the beam elements I took Harrison's code, and then started looking at the final numbers -- I noticed the third beam has some numbers negative when the book says positive.</P>
<P>Turns out that in the book Harrison draws the beam going down the negative z axis - gamma = -90 , but the sample has gamma, but the coded example has gamma = -90</P>
<P>Ah the joys of coding.</P>
<P>John</P>Thu, 23 Jun 2016 19:08:23 GMThttps://community.intel.com/t5/Intel-Fortran-Compiler/Co-ordinate-transformations/m-p/1106892#M128293JohnNichols2016-06-23T19:08:23Z