Re: Faddeeva

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

 






On Wednesday, December 14, 2016 6:46 AM, Patrick Dupre <pdupre@xxxxxxx> wrote:



Thank very good

However, there is something that I do not understand

I define
W1=0.2
z(x) = (x + W1*{0,1})
w(x) = faddeeva (z (x))
and
Voi (x) = voigt (x, W1)

Voi (x) and real (w(x)) are exactly the same
Now,
d_w (x) = - 2 * z (x) * w (z (x)) + 2*{0,1}/sqrt(pi) " which is the derivative
D_voi(x,del) = (Voi (x+del) - Voi (x-del)) / 2 / del

now 
plot [-4:4] D_voi(x,0.001)
and
plot [-4:4] real (d_w (x))

are different except for W1=0

What am I doing wrong?
It seems that the issue is with the faddeeva, or when I my calls!


Thank.

===========================================================================
Patrick DUPRÉ                                 | | email: pdupre@xxxxxxx
Laboratoire de Physico-Chimie de l'Atmosphère | |
Université du Littoral-Côte d'Opale           | |
Tel.  (33)-(0)3 28 23 76 12                   | | Fax: 03 28 65 82 44
189A, avenue Maurice Schumann                 | | 59140 Dunkerque, France
===========================================================================
_______________________________________________
users mailing list -- users@xxxxxxxxxxxxxxxxxxxxxxx
To unsubscribe send an email to users-leave@xxxxxxxxxxxxxxxxxxxxxxx



I am now a bit confused.  {0,1} represents the number 0 + 1i correct?
Then I would believe that Mr. Howell is correct in that the magnitude should be called upon since to get a real number function.  His comment
<quote>
I am not a mathematician, but I have programmed extensively
with complex numbers.  W(x)=faddeeva(z(x)) will return a complex,
right?
        The magnitude of a complex is its full value, so you would need
the sqrt(real^2+imag^2) would you not?

</quote>
z(x) is the one defined with the z(x) = (x + W1*{0,1}) which would be z(x) = x + 0*i which is just x because the W is 0.  If the function is defined correctly and I have a strong hunch that it is, then there is a problem with gnuplot and it has to be a bug somewhere.  We can try plotting the function using plotshare.com which is a frontend to gnuplot. 


There is no output from the plots check in 


http://www.plotshare.com/index.ws/plot/432515199
It actually plots one function only and the other one does not come up.  


Best Regards,



Antonio A. Olivares
_______________________________________________
users mailing list -- users@xxxxxxxxxxxxxxxxxxxxxxx
To unsubscribe send an email to users-leave@xxxxxxxxxxxxxxxxxxxxxxx



[Index of Archives]     [Older Fedora Users]     [Fedora Announce]     [Fedora Package Announce]     [EPEL Announce]     [EPEL Devel]     [Fedora Magazine]     [Fedora Summer Coding]     [Fedora Laptop]     [Fedora Cloud]     [Fedora Advisory Board]     [Fedora Education]     [Fedora Security]     [Fedora Scitech]     [Fedora Robotics]     [Fedora Infrastructure]     [Fedora Websites]     [Anaconda Devel]     [Fedora Devel Java]     [Fedora Desktop]     [Fedora Fonts]     [Fedora Marketing]     [Fedora Management Tools]     [Fedora Mentors]     [Fedora Package Review]     [Fedora R Devel]     [Fedora PHP Devel]     [Kickstart]     [Fedora Music]     [Fedora Packaging]     [Fedora SELinux]     [Fedora Legal]     [Fedora Kernel]     [Fedora OCaml]     [Coolkey]     [Virtualization Tools]     [ET Management Tools]     [Yum Users]     [Yosemite News]     [Gnome Users]     [KDE Users]     [Fedora Art]     [Fedora Docs]     [Fedora Sparc]     [Libvirt Users]     [Fedora ARM]
  Powered by Linux