<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40"><head><meta http-equiv=Content-Type content="text/html; charset=utf-8"><meta name=Generator content="Microsoft Word 15 (filtered medium)"><!--[if !mso]><style>v\:* {behavior:url(#default#VML);}
o\:* {behavior:url(#default#VML);}
w\:* {behavior:url(#default#VML);}
.shape {behavior:url(#default#VML);}
</style><![endif]--><style><!--
/* Font Definitions */
@font-face
        {font-family:"Cambria Math";
        panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
        {font-family:Calibri;
        panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
        {font-family:Verdana;
        panose-1:2 11 6 4 3 5 4 4 2 4;}
@font-face
        {font-family:Tahoma;
        panose-1:2 11 6 4 3 5 4 4 2 4;}
@font-face
        {font-family:"DejaVu Sans";
        panose-1:0 0 0 0 0 0 0 0 0 0;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
        {margin:0in;
        margin-bottom:.0001pt;
        font-size:12.0pt;
        font-family:"Times New Roman",serif;}
a:link, span.MsoHyperlink
        {mso-style-priority:99;
        color:blue;
        text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
        {mso-style-priority:99;
        color:purple;
        text-decoration:underline;}
p
        {mso-style-priority:99;
        mso-margin-top-alt:auto;
        margin-right:0in;
        mso-margin-bottom-alt:auto;
        margin-left:0in;
        font-size:12.0pt;
        font-family:"Times New Roman",serif;}
span.EmailStyle18
        {mso-style-type:personal-reply;
        color:black;}
.MsoChpDefault
        {mso-style-type:export-only;
        font-size:10.0pt;}
@page WordSection1
        {size:8.5in 11.0in;
        margin:1.0in 1.0in 1.0in 1.0in;}
div.WordSection1
        {page:WordSection1;}
/* List Definitions */
@list l0
        {mso-list-id:1261987007;
        mso-list-template-ids:-1977056432;}
ol
        {margin-bottom:0in;}
ul
        {margin-bottom:0in;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]--></head><body lang=EN-US link=blue vlink=purple><div class=WordSection1><p class=MsoNormal><span style='color:black'>Hi Al<o:p></o:p></span></p><p class=MsoNormal><span style='color:black'><o:p> </o:p></span></p><p class=MsoNormal><span style='color:black'>I couldn’t agree more.<o:p></o:p></span></p><p class=MsoNormal><span style='color:black'><o:p> </o:p></span></p><p class=MsoNormal><span style='color:black'>One thing that observations of the principles of physics has illustrated for me is that the problems all boil down to very basic and easy to understand starting principles.  <o:p></o:p></span></p><p class=MsoNormal><span style='color:black'>In this way I agree very much with John Archibald Wheeler when he said…</span><i><span style='font-size:10.0pt;font-family:"Arial",sans-serif;color:#333333'>“It is my opinion that everything must be based on a simple idea. And it is my opinion that this idea, once we have finally discovered it, will be so compelling, so beautiful, that we will say to one another, “yes, how could it have been any different?””<o:p></o:p></span></i></p><p class=MsoNormal><i><span style='font-size:10.0pt;font-family:"Arial",sans-serif;color:#333333'><o:p> </o:p></span></i></p><p class=MsoNormal><span style='font-size:10.0pt;font-family:"Arial",sans-serif;color:#333333'><o:p> </o:p></span></p><p class=MsoNormal><span style='font-size:10.0pt;font-family:"Arial",sans-serif;color:#333333'>There are some questions which might help us in this quest…<o:p></o:p></span></p><p class=MsoNormal><span style='font-size:10.0pt;font-family:"Arial",sans-serif;color:#333333'><o:p> </o:p></span></p><p class=MsoNormal><span style='font-size:10.0pt;font-family:"Arial",sans-serif;color:#333333'>When it is difficult and complex to really get the basic concepts to start with, is it likely we are on the right path?  <o:p></o:p></span></p><p class=MsoNormal><span style='font-size:10.0pt;font-family:"Arial",sans-serif;color:#333333'>Are we really explaining the concepts well? Or is it rather that we MUST manufacture artificially complex solutions in order to keep the basis we trust to be correct? <o:p></o:p></span></p><p class=MsoNormal><span style='font-size:10.0pt;font-family:"Arial",sans-serif;color:#333333'>Have we done our job and tested and questioned the foundations upon which we are building?<o:p></o:p></span></p><p class=MsoNormal><span style='font-size:10.0pt;font-family:"Arial",sans-serif;color:#333333'>Can we expect to build the correct theory from the existing foundations? Or is it the foundations which are preventing our further understanding? <o:p></o:p></span></p><p class=MsoNormal><span style='font-size:10.0pt;font-family:"Arial",sans-serif;color:#333333'>Is it errors in the foundations which has caused the stagnation in physics? <o:p></o:p></span></p><p class=MsoNormal><span style='font-size:10.0pt;font-family:"Arial",sans-serif;color:#333333'>Is it possible that the proposed solutions are so complex simply because they are wrong? And again, are they wrong because we have misunderstood some of the foundational principles of nature?<o:p></o:p></span></p><p class=MsoNormal><span style='font-size:10.0pt;font-family:"Arial",sans-serif;color:#333333'><o:p> </o:p></span></p><p class=MsoNormal><span style='font-size:10.0pt;font-family:"Arial",sans-serif;color:#333333'>Is it possible that when we do understand correctly the basic principles of nature, that the solutions will be much easier, paradox free, coherent, and quite unified?  I think the answer to this question is a resounding yes.<o:p></o:p></span></p><p class=MsoNormal><span style='font-size:10.0pt;font-family:"Arial",sans-serif;color:#333333'><o:p> </o:p></span></p><p class=MsoNormal><span style='font-size:10.0pt;font-family:"Arial",sans-serif;color:#333333'>This is not to say that any specific theory is wrong.  But it is at least a healthy approach to attempt to keep us well founded in our quest. But it takes work.<o:p></o:p></span></p><p class=MsoNormal><span style='font-size:10.0pt;font-family:"Arial",sans-serif;color:#333333'><o:p> </o:p></span></p><p class=MsoNormal><span style='font-size:10.0pt;font-family:"Arial",sans-serif;color:#333333'>Chip</span><span style='color:black'><o:p></o:p></span></p><p class=MsoNormal><span style='color:black'><o:p> </o:p></span></p><div><div style='border:none;border-top:solid #E1E1E1 1.0pt;padding:3.0pt 0in 0in 0in'><p class=MsoNormal><b><span style='font-size:11.0pt;font-family:"Calibri",sans-serif'>From:</span></b><span style='font-size:11.0pt;font-family:"Calibri",sans-serif'> General [mailto:general-bounces+chipakins=gmail.com@lists.natureoflightandparticles.org] <b>On Behalf Of </b>af.kracklauer@web.de<br><b>Sent:</b> Friday, November 13, 2015 8:41 AM<br><b>To:</b> general@lists.natureoflightandparticles.org<br><b>Cc:</b> Stephen Leary <sleary@vavi.co.uk>; Mark, Martin van der <martin.van.der.mark@philips.com>; Nature of Light and Particles - General Discussion <general@lists.natureoflightandparticles.org>; pete@leathergoth.com; David Williamson <david.williamson@ed.ac.uk>; Nicholas Bailey <Nicholas.Bailey@glasgow.ac.uk><br><b>Subject:</b> Re: [General] SU(2) equation set<o:p></o:p></span></p></div></div><p class=MsoNormal><o:p> </o:p></p><div><div><div><p class=MsoNormal><span style='font-size:9.0pt;font-family:"Verdana",sans-serif'>Hi:<o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:9.0pt;font-family:"Verdana",sans-serif'> <o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:9.0pt;font-family:"Verdana",sans-serif'>If I may, I'd like to stick a word here in.  For the record (as is said): I'm not up to speed, although I've given it a once-over.  <o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:9.0pt;font-family:"Verdana",sans-serif'> <o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:9.0pt;font-family:"Verdana",sans-serif'>In the past, goofy, incomplete, mystical Physics theories have resulted from faulty or contradictory or inapproporate (and usually also covert-implicit) input.  Then the theory get developed, embelished, expanded, partially verified and what not until its "too big to fail"!  An antidote for this syndrom is FIRST to produce a proposal or schimatic "for dummies."  One that "your mother could understand."  It's hard to produce these things, because one has to return to a stage of understanding, so as to speak to a 'dummy' on his level. that the typical proponent has long left behind.  But the great advantage is that, it is hard to fudge the basic inputs on the basis of known, but sophisticated knowledge, information that in principle should not be referenced or input at the initial level, so as to keep the logic clean.<o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:9.0pt;font-family:"Verdana",sans-serif'> <o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:9.0pt;font-family:"Verdana",sans-serif'>If this step is ommited, and you are very lucky (at first!), you may find yourself way, way, out on the limb with a lot of *** on your face!<o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:9.0pt;font-family:"Verdana",sans-serif'> <o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:10.0pt;font-family:"Verdana",sans-serif'>For what it is worth, Al   </span><span style='font-size:9.0pt;font-family:"Verdana",sans-serif'><o:p></o:p></span></p></div><div><p class=MsoNormal><span style='font-size:9.0pt;font-family:"Verdana",sans-serif'>  <o:p></o:p></span></p><div style='border:none;border-left:solid #C3D9E5 1.5pt;padding:0in 0in 0in 8.0pt;margin-left:7.5pt;margin-top:7.5pt;margin-right:3.75pt;margin-bottom:3.75pt;word-wrap: break-word;-webkit-nbsp-mode: space;-webkit-line-break: after-white-space' name=quote><div style='margin-bottom:7.5pt'><p class=MsoNormal><b><span style='font-size:9.0pt;font-family:"Verdana",sans-serif'>Gesendet:</span></b><span style='font-size:9.0pt;font-family:"Verdana",sans-serif'> Freitag, 13. November 2015 um 14:57 Uhr<br><b>Von:</b> "John Williamson" <<a href="mailto:John.Williamson@glasgow.ac.uk">John.Williamson@glasgow.ac.uk</a>><br><b>An:</b> "Nicholas Bailey" <<a href="mailto:Nicholas.Bailey@glasgow.ac.uk">Nicholas.Bailey@glasgow.ac.uk</a>><br><b>Cc:</b> "Stephen Leary" <<a href="mailto:sleary@vavi.co.uk">sleary@vavi.co.uk</a>>, "David Williamson" <<a href="mailto:david.williamson@ed.ac.uk">david.williamson@ed.ac.uk</a>>, "Nature of Light and Particles - General Discussion" <<a href="mailto:general@lists.natureoflightandparticles.org">general@lists.natureoflightandparticles.org</a>>, "<a href="mailto:pete@leathergoth.com">pete@leathergoth.com</a>" <<a href="mailto:pete@leathergoth.com">pete@leathergoth.com</a>>, "Mark, Martin van der" <<a href="mailto:martin.van.der.mark@philips.com">martin.van.der.mark@philips.com</a>><br><b>Betreff:</b> Re: [General] SU(2) equation set<o:p></o:p></span></p></div><div name=quoted-content><div><div><p class=MsoNormal><span style='font-size:10.0pt;font-family:"Tahoma",sans-serif;color:black'>Dear everyone,<br><br>Chip - your background is ideal. I always tell my aero students (I teach them a maths course) - that they will have a better skill set than physicists or mathematicians for solving differential equations because fluid mechanics is (relatively) hard, compared to mere Maxwell or QM (even relativistic QM!). Also the algebra (The space-time algebra (STA)- Cl(1,3)) is a kind of 4D extension of the quaternion algebra (which it contains as a sub-algebra). This is not to say it will be easy for you - there are concepts in there which will still make your head explode.<br><br>First of these is: the quaternion algebra is a division algebra (technically a division ring). The STA is not. This is both a problem and a solution. It is a problem in that there are limits where differentiation explodes. It is a solution in that division has implicit scalings (for which we need Martin and my paper on division) which allow (sorry -force!) extra constraints which are required - as the Maxwell equations alone) are under-constrained. Same is true for ordinary QM by the way - as one needs (at least) an extra external normalisation condition. Note, for example, that my eq 21 does not contain one - it is self-quantised.<br><br>Coming back to the code. I'm delighted that you professionals don't really seem to mind on which platform it runs. My preference would be for Python, or failing that something I could, at least read, like C. Object-oriented stuff makes my poor old head hurt!<br><br>Hmm on speed ... fast is good of course, but I think i would rather it was clever. What I am looking for is generating solutions to the particle spectrum. I have an analytic solutions for the photon, and a candidate one for the electron and the positron.  This would allow for testing of the numerical program in those cases. In the first instance a numerical system would look to parallel these.<br><br>What is needed further is a model for the quarks (I think I know how to do this, as I talked about first at CYBCOM 2008) and the interactions between them. Also one needs a model for the neutrinos - and how they propagate. Martin and I have some ideas for this - but there is room for (a lot of) development here. The best thing would be to get all of us into a room and not come out till we had a specification. I think we need a week in the first instance. Starting with me and Martin teaching you about how things work! Despite the existing papers, email back and forwards is just going to be inadequate<br><br>To brass tacks: a good starting point (even better than Stephen's code - as it is more extensive) for the algebra input is Martin's 16 by 16 matrix.<br><br>The second thing is a sharing of the division algebra paper (are we ready for this Martin?). This is a second-level thing though - and not needed for the first steps in the photon solution.<br><br>Anyway Stephen and Martin - this is down to you guys to decide whether or not (and how much) of this stuff you want to make public.<br><br>The next is the development of a proper "grid". A merely spatial (xyz) grid is not going to hack this. Indeed, it will be counter productive since it imposes a form which is just too simple (and one is already lost). What the 4-vector differential actually DOES (in terms of process - per differential) is transform a set of eight even terms (the fields, root-mass and quadrivector) to a set of odd terms (the currents and the angular momenta) and vice -versa. The constraint here is that the process must be UNITARY.  This means going to transformations (in terms of process) from an eight degree of freedom system to another eight degree of freedom system, using unitary differential transformations. 8D to 8D then. Not for the faint-hearted! Any "motion" or "oscillation" is really a transformation between these sets. I have invented a mathematics to try to think about this, based on fundamental process. I cannot email it to you as it is symbolic. It could be implemented (I think) as a process in a program. I anyway do not want to just "give it away", publically. Best method to teach it in the first instance is on a whiteboard!<br><br>The final step, however, is that everything must project back onto a grid on which one can measure it. That is 3D space and time - the basis set.<br><br>That is it. The process should parallel the process observed in reality - and project onto the basis of what may be measured. Hopefully, then, what is observed parallels what is simulated.<br><br>That brings me to a second question. Are we going to make the code and the process public or private? If private, under the "ownership" of what body? There are a few considerations here. Firstly, I noted recently that there is a prize for the solutions of a set of problems in mathematics - the Clay prizes. Solve one, win a million (dollars - i think). I got fairly excited, because I thought I may be able to solve one of them with (aspects of) the new maths. I phoned Nick Bailey for help - and he and I had a session last Saturday. Upshot of that was - no - probably cannot solve it in short order - partly because the way the question is formulated is on an un-realistic basis. Best that can be hoped for is to suggest "oracle" like solutions to particular practical problems.  Another problem has to do with the "momentum gap" in Yang-Mills theories. Again - cannot solve this directly (as Yang-Mills - though close to the truth - are themselves unphysical - in my view) - but should be able to solve the ACTUAL underlying problem - so (morally) should win that one (if we can do it)). Also, if we have code modelling both light and matter, one can begin to think about simulating molecules and crystals - and hence inventing new materials, devices and systems. That can be very lucrative!<br><br>My usual feeling is that all code should be open - but if there is serious money to be made down the line, which we could put into a foundation for the support and training of young international scientists, for example, then I would like to do that.<br><br>Anyway - the good news is that Hilbert once said that physics was getting to hard for the physicists. If we can solve one or two of the Clay prizes that will tend to show that maths just got too hard for the mathematicians!<br><br>Regards, John.<br><br>  <o:p></o:p></span></p><div><div class=MsoNormal align=center style='text-align:center'><span style='color:black'><hr size=2 width="100%" align=center></span></div><div id=divRpF118720><p class=MsoNormal><b><span style='font-size:10.0pt;font-family:"Tahoma",sans-serif;color:black'>From:</span></b><span style='font-size:10.0pt;font-family:"Tahoma",sans-serif;color:black'> Nicholas Bailey<br><b>Sent:</b> Friday, November 13, 2015 11:21 AM<br><b>To:</b> John Williamson<br><b>Cc:</b> Nature of Light and Particles - General Discussion; Mark, Martin van der; Stephen Leary; Joakim Pettersson<br><b>Subject:</b> Re: [General] SU(2) equation set</span><span style='color:black'><br> <o:p></o:p></span></p></div><div><p class=MsoNormal><span style='color:black'> <o:p></o:p></span></p></div><div><p style='margin:0in;margin-bottom:.0001pt'><span style='color:black'>Do you want it to be <i>fast</i>?<o:p></o:p></span></p><p style='margin:0in;margin-bottom:.0001pt'><span style='color:black'> <o:p></o:p></span></p><p style='margin:0in;margin-bottom:.0001pt'><span style='color:black'>Python is fast (and fun!) to write but slow to run. Although actually, not that slow. <a href="http://percival-music.ca/artifastring/" target="_blank"><span style='color:#0057AE'>Artifastring</span></a>, a modal physical model written by my ex-postgrad in C++ as the engine behind <a href="http://percival-music.ca/vivi.html" target="_blank"><span style='color:#0057AE'>Vivi</span></a>, the virtual violinist, actually ran faster when bound to Python than when linked to a larger C++ program! I didn't look into why. Curious.<o:p></o:p></span></p><p style='margin:0in;margin-bottom:.0001pt'><span style='color:black'> <o:p></o:p></span></p><p style='margin:0in;margin-bottom:.0001pt'><span style='color:black'>Anyway, if you want fast, I'd write in C++ with as much built-in operator overloading as you can so that you can read your own code. You can test that to hell, then pull it in to Python using <a href="http://www.swig.org/" target="_blank"><span style='color:#0057AE'>swig</span></a> or <a href="https://riverbankcomputing.com/software/sip/intro" target="_blank"><span style='color:#0057AE'>sip</span></a>.<o:p></o:p></span></p><p style='margin:0in;margin-bottom:.0001pt'><span style='color:black'> <o:p></o:p></span></p><p style='margin:0in;margin-bottom:.0001pt'><span style='color:black'>Alternatively, I suppose one could think of starting with Python and extending the functionality of a numeric module such as <a href="http://www.numpy.org/" target="_blank"><span style='color:#0057AE'>numpy</span></a>.<o:p></o:p></span></p><p style='margin:0in;margin-bottom:.0001pt'><span style='color:black'> <o:p></o:p></span></p><p style='margin:0in;margin-bottom:.0001pt'><span style='color:black'>The other advantage of having at least the possibility of running C++ instead of Python is that Python doesn't support threads (because of the <a href="https://www.jeffknupp.com/blog/2012/03/31/pythons-hardest-problem/" target="_blank"><span style='color:#0057AE'>Global Intepreter Lock</span></a> problem) and if you want to do real heavy duty number crunching on multiple cores, you'd best use C++ with <a href="http://www.boost.org/" target="_blank"><span style='color:#0057AE'>boost</span></a> or something like that. Boost is proposed to be included in the standard library in C++17 anyway.<o:p></o:p></span></p><p style='margin:0in;margin-bottom:.0001pt'><span style='color:black'> <o:p></o:p></span></p><p style='margin:0in;margin-bottom:.0001pt'><span style='color:black'>Nick/.<o:p></o:p></span></p><p style='margin:0in;margin-bottom:.0001pt'><span style='color:black'> <o:p></o:p></span></p><p style='margin:0in;margin-bottom:.0001pt'><span style='color:black'>On Thursday 12 November 2015 13:37:03 John Williamson wrote:<o:p></o:p></span></p><p style='mso-margin-top-alt:9.0pt;margin-right:30.0pt;margin-bottom:9.0pt;margin-left:30.0pt'><span style='color:black'>Firstly need to implement a system with primitive operations which parallel the workings of the space-time algebra. Multiplications, divisions and differentiations in particular.  It would be nice if such a framework existed already, but I do not think so. Stephen has written a suite in PYTHON to deal with the multiplications, which may be a start. Not sure what is best here: C++?, Python?, Something else?<o:p></o:p></span></p><p style='margin:0in;margin-bottom:.0001pt'><span style='color:black;background:white'>Hi John and Martin,<br><br>I am thrilled with the opportunities now when John has raised the funding question and Chip finally took a dive into the math and calculus!<br>Hope it stays at that level or continues to the next!<br><br>Some things I believe I could be most helpful with in your work are listed below. Tell me what to do and I will work free on it a few hours a week, odd week numbers. I could also offer my services more intensively as a consultant if that is an option. My normal tariff is 850 SEK/h (pretty standard in Sweden) and I only charge for productive time (not setting up tools, training, coffe-breaks etc). Likewise goes for my collegues at Join Business and Technology AB where there are few people that already work with FEM modeling and simulations regularly.</span><span style='color:black'><br> <o:p></o:p></span></p><ol start=1 type=1><li class=MsoNormal style='color:black;mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;mso-list:l0 level1 lfo1'>simpy/numpy Python framework for interactive calculus:<o:p></o:p></li><li class=MsoNormal style='color:black;mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;mso-list:l0 level1 lfo1'><a href="https://github.com/joakimbits/Quflow-and-Perfeco-tools/blob/master/lib/algebra.py" target="_blank">https://github.com/joakimbits/Quflow-and-Perfeco-tools/blob/master/lib/algebra.py</a> (outdated now but I can fix it to work in Python 3)<o:p></o:p></li><li class=MsoNormal style='color:black;mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;mso-list:l0 level1 lfo1'>Ejs Java framework for interactive 2D-3D simulations (and a clever web interface): <a href="http://www.um.es/fem/EjsWiki/Main/ExamplesLorentzForce" target="_blank">http://www.um.es/fem/EjsWiki/Main/ExamplesLorentzForce</a><o:p></o:p></li><li class=MsoNormal style='color:black;mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;mso-list:l0 level1 lfo1'>Usage of DX12 or its up-coming open-source derivates to make rough FEM/numeric calculations on arbitary gpu/cpu platforms.<o:p></o:p></li><li class=MsoNormal style='color:black;mso-margin-top-alt:auto;mso-margin-bottom-alt:auto;mso-list:l0 level1 lfo1'>I could also bring my dual quad-core 32GB RAM server online if it helps.<o:p></o:p></li></ol><p class=MsoNormal><span style='color:black'>Which calculus framework are you experimenting with today? I used to work in Mathematica, Matlab and Mathcad but nowadays I use Python or Jython for pretty much everything.<br><br>Best regards,<br>/Joakim <o:p></o:p></span></p><p><span style='color:black'> <o:p></o:p></span></p><div id=divBdy><div><div><div><div><p class=MsoNormal><span style='color:black'>Hi John W.<o:p></o:p></span></p></div><div><p class=MsoNormal><span style='color:black'> <o:p></o:p></span></p></div><div><p class=MsoNormal><span style='color:black'>Some of my background is in aerodynamic simulation, where I have written more than 500K lines of ‘C’ code and created one of the world’s most accurate real-time simulations of flight. Lots of use of quaternions and fluid mechanics of course.<o:p></o:p></span></p></div><div><p class=MsoNormal><span style='color:black'> <o:p></o:p></span></p></div><div><p class=MsoNormal><span style='color:black'>I have written code in many different languages, so practically any we choose will work for me.<o:p></o:p></span></p></div><div><p class=MsoNormal><span style='color:black'> <o:p></o:p></span></p></div><div><p class=MsoNormal><span style='color:black'>In order to start we will need to do as you suggest and look at the fundamentals to create a library of functions we can call as required.<o:p></o:p></span></p></div><div><p class=MsoNormal><span style='color:black'> <o:p></o:p></span></p></div><div><p class=MsoNormal><span style='color:black'>So if I can look at Stephen’s Python code for multiplications that would help me get a handle on how to write the division and differentiation low level code.<o:p></o:p></span></p></div><div><p class=MsoNormal><span style='color:black'> <o:p></o:p></span></p></div><div><p class=MsoNormal><span style='color:black'>Some personal experience and observations…<o:p></o:p></span></p></div><div><p class=MsoNormal><span style='color:black'>Object oriented code generally takes more time to design and runs a little slower, but it protects and manages the data and functions operating on the data as individual objects. Object oriented code generally uses more memory as well. Computation time, and memory capacity, are really critical for us because of the tremendous computational load it will take to do the simulation at any useful scale. So when we were writing flight simulation code, which had to compute all the forces and moments in real time, at about 120 times a second, we wrote it all in ‘C’.  I am bringing this up because even now, when I am doing physics modeling, which is much simpler than the task at hand, I often have to wait an hour for the results from our fairly large and very fast workstations here in the office. It would be nice to create the simulation in the fastest running code we can, so that we have some hope of getting results in reasonable time, without having to rent expensive time on, or build a supercomputer.<o:p></o:p></span></p></div><div><p class=MsoNormal><span style='color:black'> <o:p></o:p></span></p></div><div><p class=MsoNormal><span style='color:black'>But those opinions are not as important as getting the foundation in place so we can start doing some real work in simulating particles, and interactions.  Exciting stuff.<o:p></o:p></span></p></div><div><p class=MsoNormal><span style='color:black'> <o:p></o:p></span></p></div><div><p class=MsoNormal><span style='color:black'>Chip<o:p></o:p></span></p></div></div></div></div></div><p style='margin:0in;margin-bottom:.0001pt'><span style='color:black'> <o:p></o:p></span></p></div></div></div><p class=MsoNormal><span style='font-size:9.0pt;font-family:"DejaVu Sans",serif'>_______________________________________________ If you no longer wish to receive communication from the Nature of Light and Particles General Discussion List at <a href="mailto:af.kracklauer@web.de">af.kracklauer@web.de</a> <a href="http://lists.natureoflightandparticles.org/options.cgi/general-natureoflightandparticles.org/af.kracklauer%40web.de?unsub=1&unsubconfirm=1" target="_blank">Click here to unsubscribe </a><o:p></o:p></span></p></div></div></div></div></div></div></div></body></html>