Hi, thanks for examples to make scaled units (nm, km, GJ, etc..). I'll get to this soon. BTW, sorry for sending previous email twice. I pressed "send" keyboard shortcut too early. Today I wanted to finish my quantum wavefunction harmonic oscillator. I was working on it, while in parallel compiling & trying to install boost 1.48, which currently was unsuccessful (debian squeeze). BTW, if anybody has a hint about how to install boost 1.48 on debian squeeze using preferably using the "legal" `dpkg-buildpackage -rfakeroot -b` way, it would be great. But still I am quite happy today, because I managed to get a correctly working plots of this harmonic oscillator, but without units, unfortunately. One small problem is that I couldn't get to higher energies because double precision has limited factorials. I would need a working factorial of 100! which currently seems unlikely. And my currently maximum is only 25! I didn't check yet, but maybe boost has some arbitrary precision library, which would work with Boost.Units and boost::math::factorial ? So, for your pleasure I am presenting you with a harmonic oscillator code :) Unfortunately without units :) But if you will try to convert this to units, I might learn even more about using this library. Especially because dimension of wavefunction in 1 dimensional space is 1/sqrt(meter), and I have no idea how to declare this derived_dimension, this didn't work: typedef derived_dimension<length_base_dimension,root<2> >::type quantum_wavefunction_1D; In the attachment you will find a plot and code for frequency=1Hz, hbar=1 J*s, mass=1kg (plotted against a classical harmonic oscillator). Which you can compare with wikipedia and admire how similar they are :) (note that I take Planck's constant to be 1, for start). So things to do: - install boost 1.48 on my PC - make promised scaled SI units - convert this harmonic oscillator code to use Boost.Units - find some arbitrary precision library, that will work with Units and VERY big factorials, like 100! or even better 1000! I hope you enjoy this code. I like it :) best regards -- Janek Kozicki http://janek.kozicki.pl/ |
