o o h?@sddlmZddlmZddlmZddlmZddlm Z ddl m Z ddl m Z ddlmZmZmZmZmZmZdd lmZdd lmZdd lmZmZmZdd lmZdd lmZddl m!Z!GdddeZ"ddZ#dS))Rational)S)symbols)sign)sqrt)gcd) Complement)BasicTuplediffexpandEqInteger)ordered)_symbol)solveset nonlinsolve diophantine total_degree)Point)corecsveZdZdZfddZeddZeddZedd Zd d Z d d Z ddZ ddZ ddZ dddZZS)ImplicitRegiona Represents an implicit region in space. Examples ======== >>> from sympy import Eq >>> from sympy.abc import x, y, z, t >>> from sympy.vector import ImplicitRegion >>> ImplicitRegion((x, y), x**2 + y**2 - 4) ImplicitRegion((x, y), x**2 + y**2 - 4) >>> ImplicitRegion((x, y), Eq(y*x, 1)) ImplicitRegion((x, y), x*y - 1) >>> parabola = ImplicitRegion((x, y), y**2 - 4*x) >>> parabola.degree 2 >>> parabola.equation -4*x + y**2 >>> parabola.rational_parametrization(t) (4/t**2, 4/t) >>> r = ImplicitRegion((x, y, z), Eq(z, x**2 + y**2)) >>> r.variables (x, y, z) >>> r.singular_points() EmptySet >>> r.regular_point() (-10, -10, 200) Parameters ========== variables : tuple to map variables in implicit equation to base scalars. equation : An expression or Eq denoting the implicit equation of the region. cs8t|ts t|}t|tr|j|j}t|||SN) isinstancer r lhsrhssuper__new__)cls variablesequation __class__o/var/www/html/construction_image-detection-poc/venv/lib/python3.10/site-packages/sympy/vector/implicitregion.pyr9s   zImplicitRegion.__new__cC |jdS)Nrargsselfr$r$r%r B zImplicitRegion.variablescCr&)Nr'r)r$r$r%r!Fr+zImplicitRegion.equationcCs t|jSr)rr!r)r$r$r%degreeJr+zImplicitRegion.degreec CsZ|j}t|jdkrtt||jdtjddfSt|jdkrT|jdkrTt|j|}\}}}}}}|dd||krI|j |\} } | | fS|j |\} } | | fSt|jdkr|j\} } } t ddD]3} t ddD]+} t| | | | | i|jdtjdj s| | tt| | | | | idfSqmqft|dkrt|dSt) a0 Returns a point on the implicit region. Examples ======== >>> from sympy.abc import x, y, z >>> from sympy.vector import ImplicitRegion >>> circle = ImplicitRegion((x, y), (x + 2)**2 + (y - 3)**2 - 16) >>> circle.regular_point() (-2, -1) >>> parabola = ImplicitRegion((x, y), x**2 - 4*y) >>> parabola.regular_point() (0, 0) >>> r = ImplicitRegion((x, y, z), (x + y + z)**4) >>> r.regular_point() (-10, -10, 20) References ========== - Erik Hillgarter, "Rational Points on Conics", Diploma Thesis, RISC-Linz, J. Kepler Universitat Linz, 1996. Available: https://www3.risc.jku.at/publications/download/risc_1355/Rational%20Points%20on%20Conics.pdf r,r)domaini )r!lenr listrrRealsr- conic_coeff_regular_point_parabola_regular_point_ellipserangesubsis_emptysingular_pointsNotImplementedError)r*r!coeffsabcdefx_regy_regxyzr$r$r% regular_pointNs,   &,zImplicitRegion.regular_pointc Cs"||fdko||fdko|dd||ko||fdk}|s"td|dkrUd||d||d|||d}} |dkrR| |} |||  d|} n5d}n2|dkrd||d||d|||d}} |dkr| |} |||  d|} nd}|r| | fStd)N)rrr/r0*Rational Point on the conic does not existrF) ValueError) r*r?r@rArBrCrDokd_dashf_dashrFrEr$r$r%r7s$8. . z&ImplicitRegion._regular_point_parabolac.s4d|||d}|}|std|dkr'|dkr'd} d||||} nY|dkr`|} d|d|dd||||d|||dd|d||d||d|} n |} d|d|dd||||d|d||} | dko| dko| dk }|stdt| d} t| d} | j| j} } | j| j} }t| |}|| |}| ||}| |  |}t|tt|d}t ||}t|tt|d}t ||}t|tt|d}t ||}tt|||}||}||}||}t||}||}||}||}t||}||}||}||}t||}||}||}||}t d\}}}||d||d||d}t |} t | dkritdd }!| D]}"t |"j}#t|#d |"d}$|$dkrd }!qmt|$ttfs|$j}%t |%d krtt|%}&ttjtt|$d|&tj}'tt|'|&<t |%dkrtt|%\}&}(tjD])})|$|&|)}*ttjtt|*d|(tj}+|+js|)|&<tt|+|(<nqt |#dkrtfd d|"D\}}}n|"\}}}d }!nqm|!rtd|||}|||}|||}||}||}|dkrY|dkrY||d|d|},||d|d|}-|,|-fS|dkr||d||||| },|||,|d|}-|,|-fS|d||||| }-|||-|d|},|,|-fS)Nr0r/rKrlJ)zx y zFr1Tr,c3|]}|VqdSrr:.0srepr$r% z8ImplicitRegion._regular_point_ellipse..) rLrlimit_denominatorpqrrrabsrrrr3r free_symbolsdictfromkeysrintrnextiterrrIntegersrr r4rr:r;tuple).r*r?r@rArBrCrDDrMKLk1k2l1l2ga1b1c1a2r1b2r2c2r3g1g2g3rGrHrIeq solutionsflagsolsymssol_zsyms_zr\p_valuesr]i subs_sol_zq_valuesrErFr$rWr%r8sf<         $           z%ImplicitRegion._regular_point_ellipsecCs6|jg}|jD] }|t|j|g7}qt|t|jS)a Returns a set of singular points of the region. The singular points are those points on the region where all partial derivatives vanish. Examples ======== >>> from sympy.abc import x, y >>> from sympy.vector import ImplicitRegion >>> I = ImplicitRegion((x, y), (y-1)**2 -x**3 + 2*x**2 -x) >>> I.singular_points() {(1, 1)} )r!r r rr4)r*eq_listvarr$r$r%r<s zImplicitRegion.singular_pointscCst|tr|j}|j}t|jD]\}}|||||}qt|}t|jdkr8|j}t dd|D}|S|}t |}|S)a Returns the multiplicity of a singular point on the region. A singular point (x,y) of region is said to be of multiplicity m if all the partial derivatives off to order m - 1 vanish there. Examples ======== >>> from sympy.abc import x, y, z >>> from sympy.vector import ImplicitRegion >>> I = ImplicitRegion((x, y, z), x**2 + y**3 - z**4) >>> I.singular_points() {(0, 0, 0)} >>> I.multiplicity((0, 0, 0)) 2 rcss|]}t|VqdSrr)rUtermr$r$r%rYPsz.ImplicitRegion.multiplicity..) rrr(r! enumerater r:r r3minr)r*point modified_eqrrtermsmr$r$r% multiplicity2s zImplicitRegion.multiplicitytrVNcsr|j}|j}|dkr0t|jdkr|fSt|jdkr-|j\}}tt||d}||fStd}|dkrR|dur=|}nt|dkrNt|d}n|}t|dkr|} | D]+} t | j } t | dt| dkr~t fdd| D} || |dkr| }nq`t|dkrt|} t|jD]\} }| |||| } qt| } d}}| jD]}t||kr||7}q||7}qd|}t|t s|f}t|jdkr5|d}|d krtd d d }ntd d d }t|d d }||jd||jd|i}||jd||jd|i}||||d|d}||||d|d}||fSt|jd kr|\}}d|vrMtdd d }ntdd d }t|d d }t|d d }||jd||jd||jd|i}||jd||jd||jd|i}||||d|d}||||d|d}||||d|d}|||fSt)a Returns the rational parametrization of implicit region. Examples ======== >>> from sympy import Eq >>> from sympy.abc import x, y, z, s, t >>> from sympy.vector import ImplicitRegion >>> parabola = ImplicitRegion((x, y), y**2 - 4*x) >>> parabola.rational_parametrization() (4/t**2, 4/t) >>> circle = ImplicitRegion((x, y), Eq(x**2 + y**2, 4)) >>> circle.rational_parametrization() (4*t/(t**2 + 1), 4*t**2/(t**2 + 1) - 2) >>> I = ImplicitRegion((x, y), x**3 + x**2 - y**2) >>> I.rational_parametrization() (t**2 - 1, t*(t**2 - 1)) >>> cubic_curve = ImplicitRegion((x, y), x**3 + x**2 - y**2) >>> cubic_curve.rational_parametrization(parameters=(t)) (t**2 - 1, t*(t**2 - 1)) >>> sphere = ImplicitRegion((x, y, z), x**2 + y**2 + z**2 - 4) >>> sphere.rational_parametrization(parameters=(t, s)) (-2 + 4/(s**2 + t**2 + 1), 4*s/(s**2 + t**2 + 1), 4*t/(s**2 + t**2 + 1)) For some conics, regular_points() is unable to find a point on curve. To calulcate the parametric representation in such cases, user need to determine a point on the region and pass it using reg_point. >>> c = ImplicitRegion((x, y), (x - 1/2)**2 + (y)**2 - (1/4)**2) >>> c.rational_parametrization(reg_point=(3/4, 0)) (0.75 - 0.5/(t**2 + 1), -0.5*t/(t**2 + 1)) References ========== - Christoph M. Hoffmann, "Conversion Methods between Parametric and Implicit Curves and Surfaces", Purdue e-Pubs, 1990. Available: https://docs.lib.purdue.edu/cgi/viewcontent.cgi?article=1827&context=cstech r,r/rr$Nc3rRrrSrTrWr$r%rYrZz:ImplicitRegion.rational_parametrization..rPrVs_T)realr1rr_)r!r-r3r r4rr=r<rJr r_r`rarfrrr:r r(rrr)r* parameters reg_pointr!r-rGrHy_parrr<spointrrrrhnhn_1r parameter1rVrx_par parameter2rz_parr$rWr%rational_parametrizationWs/                (( z'ImplicitRegion.rational_parametrization)rN)__name__ __module__ __qualname____doc__rpropertyr r!r-rJr7r8r<rr __classcell__r$r$r"r%rs '   7|%rc Cst|dkr t|d}|d}t|}||d}|||}||d}||d|d}||d|d}||d|d} |||||| fS)Nr/rr,)rrLr coeff) r r!rGrHr?r@rArBrCrDr$r$r%r6s r6N)$sympy.core.numbersrsympy.core.singletonrsympy.core.symbolr$sympy.functions.elementary.complexesr(sympy.functions.elementary.miscellaneousrsympy.polys.polytoolsrsympy.sets.setsr sympy.corer r r r r rsympy.core.sortingrr sympy.solversrrr sympy.polysrsympy.geometryrsympy.ntheory.factor_rrr6r$r$r$r%s&             _