*h)SrSSKJrJr SSKJr SSKJr "SS\5r SS /S S S S /S 4S jr S r /SQr Sr SS /SS /SS SS S S 4SjrSS /SS /SSSS S S 4Sjr\ \ l \ \ l \ \ l \\ l\\ lg )z Plotting (requires matplotlib) ) hsv_to_rgb hls_to_rgb) NoConvergence)xrangec \rSrSr\\\\4rSr g)VisualizationMethods N) __name__ __module__ __qualname____firstlineno__ ValueErrorArithmeticErrorZeroDivisionErrorr plot_ignore__static_attributes__r Q/home/mande/repo/quber/.venv/lib/python3.13/site-packages/mpmath/visualization.pyr r s0A=QKrr Nc U(aSnSn U(d%SSKn U R5n U RS5n[U[[ 45(dU/nUup/SQn [ U5GHWupURXX- [U5- 5n/n/nSn[[U55GHInUS:wa&UH nUUS- U::dMUUU:dM[e U"UU5nURU5(d[U5S:a[e[US5(ayUR(ah[UR 5n[UR5nU(dS nUR#U5 /nUR#[UU5UU45 MU(aSnUR#U5 /n[US 5(a UR nUR#[UU5U45 GML U(aUR#U5 UHnUVs/sHnUSPM nnUVs/sHnUSPM nnU(dM6X[U 5-n[US5S :XaAUVs/sHnUS PM nnUR'UUS U-S S9 UR'UUSU-S S9 MUR'UUUS S9 M GMZ UR)UVs/sHn[U5PM sn5 U(a+UR+UVs/sHn[U5PM sn5 UR-S5 UR/S5 UR1S 5 U (a(U(aW R3XVS9 gW R55 gg!UR$a U(aUR#U5 /nGM f=fs snfs snfs snfs snfs snf)a Shows a simple 2D plot of a function `f(x)` or list of functions `[f_0(x), f_1(x), \ldots, f_n(x)]` over a given interval specified by *xlim*. Some examples:: plot(lambda x: exp(x)*li(x), [1, 4]) plot([cos, sin], [-4, 4]) plot([fresnels, fresnelc], [-4, 4]) plot([sqrt, cbrt], [-4, 4]) plot(lambda t: zeta(0.5+t*j), [-20, 20]) plot([floor, ceil, abs, sign], [-5, 5]) Points where the function raises a numerical exception or returns an infinite value are removed from the graph. Singularities can also be excluded explicitly as follows (useful for removing erroneous vertical lines):: plot(cot, ylim=[-5, 5]) # bad plot(cot, ylim=[-5, 5], singularities=[-pi, 0, pi]) # good For parts where the function assumes complex values, the real part is plotted with dashes and the imaginary part is plotted with dots. .. note :: This function requires matplotlib (pylab). Nro)brgmkFrgu<7~imagTrealz--) linewidth:xzf(x)dpi)pylabfigure add_subplot isinstancetuplelist enumeratearangefloatrlenrisnanabshasattrr!r"appendrplotset_xlimset_ylim set_xlabel set_ylabelgridsavefigshow)ctxfxlimylimpointsfiler) singularitiesaxesfigr*arcolorsnfuncr'segmentssegment in_complexisingvreimsycz_s rr8r8 s>8  C llns# a% ' ' C DA &FQ< JJqac5=0 1 AA 6 -QqS6T>addl",,!.1J99Q<<3q6E>$$1f%%!&&qvvBqvvB%%)  0"$NNE!A$KR#89!%*  0"$q&))FFNNE!A$K#341 :  OOG $G&'w!1wA'&'w!1wA'3v;'A71:!##*+7aQqT7+ !QQ! 4 !QA 3 !QQ / I ` MMT*T58T*+  .AuQx./OOCOOFIIdO  MM$M ( JJL  1?? OOG, (' , +.sD.NNCNANN7-N<-OOO *N43N4c8URU5(agURU5(agSn[URU55UR-SUR-- nUS-S-nS[SS[ U5S--- 5- n[ X4S 5$) N?r\r\?r^r^(-DT! @r$r^r\r333333?皙?)isinfr4r2argpir5r)r@rXrdrIrs rdefault_color_functionreos yy|| yy|| B swwqz SVV ##&&1A SCA eAs3q63;'((A aC  r) )rhrh)gffffff)皙?g?r^)g)rhr^r\)g)g?rara)rhr[)g?)r\?r`)r^)rjr^rh)gffffff?)gffffff?rirh)r\rg)@rgcURU5(agURU5(agSn[URU55U- n[ [ US5S5n[ S[[55H\n[USU:dM[US- unupgn[Uun upn X5- X- - n XjU- U --X{U- U --XU- U --4s $ g)Nr[r]r_r\rfrr) rbr4r2rcmaxminranger3blue_orange_colors)r@rXrdwrPrIragabarrbgbbbrUs rphase_color_functionrxs yy|| yy|| B cggajBA C3KA 1S+, - a  #a '015OA|03OA| A"uai<UAIrb5!)|; ; .riFc UbUS:Xa URnUS:Xa URnSSKn U(aSn Sn U (d!U R5n U R S5n UupUupX- nX- n[ UR UU-U- 5S-5n[ UR UU-U- 5S-5nU RXU5nU RXU5nU RUUS45n[U5Hon[U5H1nURUUUU5nU"U"U55nUUUU4'M3 U(dML[[U5S -[U5-5 Mq XX4Vs/sHn[U5PM snuppU RUXX4S S 9 U R!S 5 U R#S 5 U (a(U(aU R%XxS9 gU R'5 gg!URa SnNf=fs snf)ak Plots the given complex-valued function *f* over a rectangular part of the complex plane specified by the pairs of intervals *re* and *im*. For example:: cplot(lambda z: z, [-2, 2], [-10, 10]) cplot(exp) cplot(zeta, [0, 1], [0, 50]) By default, the complex argument (phase) is shown as color (hue) and the magnitude is show as brightness. You can also supply a custom color function (*color*). This function should take a complex number as input and return an RGB 3-tuple containing floats in the range 0.0-1.0. Alternatively, you can select a builtin color function by passing a string as *color*: * "default" - default color scheme * "phase" - a color scheme that only renders the phase of the function, with white for positive reals, black for negative reals, gold in the upper half plane, and blue in the lower half plane. To obtain a sharp image, the number of points may need to be increased to 100,000 or thereabout. Since evaluating the function that many times is likely to be slow, the 'verbose' option is useful to display progress. .. note :: This function requires matplotlib (pylab). Ndefaultphaserrrr#r]z of lower)extentoriginzRe(z)zIm(z)r()rerxr*r+r,intsqrtlinspacezerosrmpcrprintstrr2imshowr;r<r>r?)r@rArSrTrDcolorverboserEr)rGr*rHrearebimaimbdredimMNr'rVrqrKrrXrRrYs rcplotrs@ }*** ((  C llns#HCHC )C )C CHHVCZ^ $Q &'A CHHVCZ^ $Q &'A s#A s#A  Q1IA AYA!ad#A $!A$KAacF  7 #a&6/CF* +.1s,@A,@q%(,@ACcKK3S.wK?OOGOOG  MM$M ( JJL  ?? $# $ Bs"G11H 1HHdTc ^SSKn SSKJn U(aSn Sn U (d+U R5n U RR U 5n UupUunnX- nUU- n[ U[[45(dXD/nUunnU RXU5nU RUUU5n[S5Vs/sHnU RUU45PM snunnn[S5Vs/sHnSS/PM snunnn[U5Hn[U5HnU"URUU5URUU55nUuUUU4'UUU4'UUU4'UUU4U4UUU4U4UUU4U44H$unn UU S:aUU S'UU S:dMUU S'M& M M U(aU RUUUSSS9 OU RUUUSSS9 U R!S5 U R#S5 U R%S 5 U(aUUU4V s/sHn U SU S- PM sn un!n"n#['U!U"U#5n$U!U$:a)U$U!- n%U R)USU%S - - USU%S - -5 U"U$:a)U$U"- n%U R+USU%S - - USU%S - -5 U#U$:a)U$U#- n%U R-USU%S - - USU%S - -5 U (a(U(aU R/XxS 9 gU R15 ggs snfs snf![a UUUUUsUUU4'UUU4'UUU4'GNf=fs sn f) a Plots the surface defined by `f`. If `f` returns a single component, then this plots the surface defined by `z = f(x,y)` over the rectangular domain with `x = u` and `y = v`. If `f` returns three components, then this plots the parametric surface `x, y, z = f(u,v)` over the pairs of intervals `u` and `v`. For example, to plot a simple function:: >>> from mpmath import * >>> f = lambda x, y: sin(x+y)*cos(y) >>> splot(f, [-pi,pi], [-pi,pi]) # doctest: +SKIP Plotting a donut:: >>> r, R = 1, 2.5 >>> f = lambda u, v: [r*cos(u), (R+r*sin(u))*cos(v), (R+r*sin(u))*sin(v)] >>> splot(f, [0, 2*pi], [0, 2*pi]) # doctest: +SKIP .. note :: This function requires matplotlib (pylab) 0.98.5.3 or higher. rNr#r)rstridecstrider'rVrXrkr()r*mpl_toolkits.mplot3dmplot3dr+axes3dAxes3Dr-r/r.rrrconvert TypeErrorplot_wireframe plot_surfacer;r< set_zlabelrm set_xlim3d set_ylim3d set_zlim3dr>r?)&r@rAurRrD keep_aspect wireframerEr)rGr*rrHuaubvavbdudvrrrPr'rVrXxabyabzabrKrfdatarWcabdxdydzmaxddeltas& rsplotrsi4*  C lln~~$$S) FB FB B bB ftUm , ,! DAq rq!A r2q!A,21I6Iqu{{Aq6"I6GAq!%+AY/YaVY/MCc AYAckk!A$'QqT):;E ;).&!A#!A#!A#acFC=1QqS63-!AaC&#G3s1v:CFs1v:CF H  Aq!Q: !Q1a8OOCOOCOOC25sCA#c!fs1voA B2r2 92IE OOCFUS[0#a&53;2F G 92IE OOCFUS[0#a&53;2F G 92IE OOCFUS[0#a&53;2F G  MM$M ( JJL  C7/  ;)*1qtU&!A#!A#!A# ;Bs$/K3 K8?K= L*=&L'&L')__doc__colorsysrrlibmpr libmp.backendrobjectr r8rerprxrrr rrrs, !R6R!f4$D4aF !  <VAt4 $TJXFr!fSd$TPf!.D+,@)""r