.fdZddlZddlZddlmZGddej ZGddej ZGdd ej Zy) zGeodesic testsN)GeodesicceZdZdZgdgdgdgdgdgdgdgd gd gd gd gd gdgdgdgdgdgdgdgdgZdZdZdZy) GeodesicTestzGeodesic test suite) gHhA@gz6>WnagQ[@g|͍Y&gVc}Qg}_A)`@gw aAgpow T@gh'$WAgB1A?g轅?g݀X[B) gиp $K@g>>tAZ@gd"6@g BS@gեh@gz+G[@g0QOAgSErB@gΪB6MAgEoڟ?gߡ?gϻb B) g|,}5g=,Ԛa@gv8@gVD@g?X@gDg `Ag FR@gWAgk?g %?g&) gx]`Pg \@g$ge@gvi)g_^}tq@g8@g_fDeAg&u$Y@guJWAg:D0ƿgmTƿgAؓ) gLqUwm1g B) g{l@@gE*-%H@g(b@gۅ:$LgwR~RIi@g^H@g3oAg5b@gGsHAgK:yg"I~gAyx) gVF@g4h^J@g;8%W@gʡEgdzi@g V40V@ghOpAg0g9Dc@g!&EAg?/g2gs) g}RHIg/L N0g1#\Zgn/i@gJiWWgٚ)7GgBm ]XAgTR M@gY!_TAgp$?g?g-XB) gP6 wMgT㥛!g[$a@gs!gS?o*R`@gcA3@g@DlfAg>]pZ@g1TNvWAg,пg{Rѿgl5) gIK5Qgy@ٔ+RgɓqW@g"lxzPIg@aӻ gZ+ SA@g@vd0NAgæ?A@g )zKLAgf?g!0 ?ggo) gq@%g$~ @gk>EOUgx&1@gwN#`glmRTg]eAg3Y@gАbvWAg5ˆtͿg_vͿgbqጎB) g%1 5gd@gpQ=@g!J]H@gB=j@g)E@g)lZaAg>}T@g=3:XAg!!>\?g-&?gJ  `@gJneYg';U&g˾+'@g] ٔUglƯhAgm)]@gqM5UAg?v\ܿg,d ܿgN,B) g}iƺ9gP6 bgC4cgJiWLg-Rpg,Hg haAg ︇*U@g ',4?XAghؒL?gnLG?gbYB) gE>'Dgףp= ^@gDb,@g4FJg|DL$L`@g~P%@gZ~MAg +A@gv/aKAgUo?gŖʪOl?gێr) gm&@g"Ha@goS@g[rP@g)?n@geYQY@gf5fAgZ@gK#B.(WAg ӿgeuӿgd$B) g*D/x=gs.UeW@gxdgmnLw;g46JQgx&{/gtiAgnUR^@g1?TAgҷ=}gg-Bc rtjD]#}|\ }}}}}}}} } } } } tjj ||||tj tj z}|j||dd|j||dd|j||dd|j||dd|j| |dd|j| |d d|j| |d d |j| |d d |j| |d d&y)z/Helper function for testing inverse calculationlon2vIh%<=deltaazi1azi2s12:0yE>a12m12M12V瞯<M21S12皙?N)r testcasesrWGS84InverseALL LONG_UNROLLassertAlmostEqual)selfllat1lon1r lat2rr r rrrrrinvs _C:\Users\user\Documents\project_loop\venv\Lib\site-packages\geographiclib/test/test_geodesic.py test_inversezGeodesicTest.test_inverse]sG  # # ;'(%tT4tT Cc3 NN " "4tT\\H$8$88 :c T3v;> T3v;> T3v;> S#e*d; S#e*e< S#e*d; S#e*e< S#e*e< S#e*c: ;c DtjD] }|\ }}}}}}}} } } } } tjj ||||tj tj z}|j||dd|j||dd|j||dd|j| |dd|j| |dd|j| |d d |j| |d d |j| |d d y)z.Helper function for testing direct calculationr rr rr rrrrrrrrN)rrrrDirectrrrrrrrr r rr r rrrrrdirects r" test_directzGeodesicTest.test_directns/  # # >'(%tT4tT Cc3~~$$T4s%-\\H4H4H%HJf T6&>5A T6&>5A T6&>5A S&-? S&-> S&-? S&-? S&-= >r$c DtjD] }|\ }}}}}}}} } } } } tjj |||| tj tj z}|j||dd|j||dd|j||dd|j||dd|j| |dd|j| |d d |j| |d d |j| |d d y)z>Helper function for testing direct calculation with arc lengthr rr rr r rrrrrrrN)rrrr ArcDirectrrrr's r"test_arcdirectzGeodesicTest.test_arcdirect~s/  # # >'(%tT4tT Cc3~~''dD#(0 x7K7K(KMf T6&>5A T6&>5A T6&>5A S&-> S&-> S&-? S&-? S&-= >r$N)__name__ __module__ __qualname____doc__rr#r)r,r$r"rrsCCCGEFECFDDCCEEECCEF[PG)d;"> >r$rceZdZdZdZdZdZdZdZdZ dZ d Z d Z d Z d Zd ZdZdZdZdZdZdZdZdZdZdZdZdZdZdZdZdZdZ dZ!d Z"d!Z#y")# GeodSolveTestzGeodSolve testsctjjdddd}|j|ddd|j|d d d|j|d d d y) GeodSolve0gLD@g33333sRgJ)""H@gffffff@r gj@+0J@h㈵>r r grt[@r i*PY?Nrrrrrr!s r"test_GeodSolve0zGeodSolveTest.test_GeodSolve0sg .. uk4 @C3v;&A3v; 6B3u:w8VD6&>7FC6&>9fEr$ctdd}|jdddd}|j|ddd |j|d dd |j|d d d |jdddd}|j|ddd |j|d dd |j|d d d y)z4Check fix for antipodal prolate bug found 2010-09-04jXAO贁N{g x#?rg x#r gk V@r6r r r 2r7rggK4V@Nrrrrgeodr!s r"test_GeodSolve2zGeodSolveTest.test_GeodSolve2s E8 $D ,,w8S 1C3v;&A3v;&A3u:x= ,,sAtS )C3v;&A3v;&A3u:x=r$cvtjjdddd}|j|dddy ) z-Check fix for short line bug found 2010-05-21g&?B@rg&?B@gHת>r g;On?Mb@?r Nr8r9s r"test_GeodSolve4zGeodSolveTest.test_GeodSolve4s8 .. !!2I ?C3u:uf=r$cRtjjdddd}|j|ddd|d dkr8|j|d d d|jt |d d dy |j|d dd|j|d ddy )z.Check fix for point2=pole bug found 2010-05-03g+@4?rcAr Zr6r rijr rBN)rrr&rabsr=s r"test_GeodSolve5zGeodSolveTest.test_GeodSolve5s ^^ " "#3RD AF6&>2v> f~ VF^T6B S0#vF VF^R@ VF^Q?r$cZtjjdddd}|j|dddtjjd dd d }|j|dd dtjjd ddd}|j|dddy)zCheck fix for volatile sbet12a bug found 2011-06-25 (gcc 4.4.4 x86 -O3). Found again on 2012-03-27 with tdm-mingw32 (g++ 4.6.1).gB V@rgB VgQdf@r gDlsArIr gP'PV@gP'PVg6*f@g]sAg*qQUV@g*qQUVgOePf@gunsANr8r9s r"test_GeodSolve6zGeodSolveTest.test_GeodSolve6s .. !!13I KC3u:|VD .. !!13I KC3u:|VD .. !!79N PC3u:|VDr$cvtjjdddd}|j|dddy ) zACheck fix for volatile x bug found 2011-06-25 (gcc 4.4.4 x86 -O3)g s)L@rg s)LgjEuf@r gPd=sArIr Nr8r9s r"test_GeodSolve9zGeodSolveTest.test_GeodSolve9s: .. !!13I KC3u:|VDr$cvtjjdddd}|j|dddy ) zRCheck fix for adjust tol1_ bug found 2011-06-25 (Visual Studio 10 rel + debug)gοW+idJ@rgͿW+idJgHMtf@r gsArIr Nr8r9s r"test_GeodSolve10zGeodSolveTest.test_GeodSolve10s< .. !!79O QC3u:|VDr$cvtjjdddd}|j|dddy ) zRCheck fix for bet2 = -bet1 bug found 2011-06-25 (Visual Studio 10 rel + debug)g]QBH@rg[QBHg G0sf@r gMb)sArIr Nr8r9s r"test_GeodSolve11zGeodSolveTest.test_GeodSolve11s< .. !!68N PC3u:|VDr$ctdd}|jdddd}|j|ddd |j|d d d |j|d d d y)zCheck fix for inverse geodesics on extreme prolate/oblate ellipsoids Reported 2012-08-29 Stefan Guenther ; fixed 2012-10-07g33333sV@gHzGrr gzG^@{Gz?r r gIZ@r g33333p@rNrDrEs r"test_GeodSolve12zGeodSolveTest.test_GeodSolve12so D% D ,,q!S# &C3v;=3v;=3u:ud;r$cNtjjdddtj}|j tj |d|j tj |d|j tj |dy)z*Check fix for inverse ignoring lon12 = nanrr r r Nrrrmathnan assertTrueisnanr9s r"test_GeodSolve14zGeodSolveTest.test_GeodSolve14sk .. Aq$(( 3COODJJs6{+,OODJJs6{+,OODJJs5z*+r$ctdd}|jddddtj}|j|ddd y ) zcInitial implementation of Math::eatanhe was wrong for e^2 < 0. This checks that this is fixed.r@rAr_ri\r7r N)rr&AREArrrFr(s r"test_GeodSolve15zGeodSolveTest.test_GeodSolve15sC E8 $D [[Aq!X]] 3F6%=%=r$ctjjddddtjtjz}|j |ddd|j |d d d|j |d d dtjj ddd}|jdtjtjz}|j |ddd|j |d d d|j |d d dtjjdddd}|j |ddd|j |d d d|j |d d d|jd}|j |ddd|j |d d d|j |d d dy)z1Check fix for LONG_UNROLL bug found on 2015-05-07(irZsAr ir_r rir ViN)rrr&STANDARDrrLinePositionrr(lines r"test_GeodSolve17zGeodSolveTest.test_GeodSolve17s ^^ " "2sC#+#4#4x7K7K#KMF6&>3:6&>4;6&>4; >>  r3 ,D ]]3 1 1H4H4H H IF6&>3:6&>4;6&>4; ^^ " "2sC 5F6&>3:6&>3:6&>4; ]]3 F6&>3:6&>3:6&>4;r$ctdd}|jddddtj}|j|ddd y ) z3>6&>4?6%=&#> ^^ # #AsAs$,$5$58L8L$LNF6&>3>6&>3>6%=&#>r$ctjjdddd}|j|ddd|j|ddd|j|dd d tjjdddd }|j|dd d|j|dd d|j|ddd tjjdddd}|j|ddd|jt |ddd|j|ddd tjjdddd}|j|ddd|jt |ddd|j|ddd tdd}|jdddd}|j|ddd|j|ddd|j|ddd |jdddd}|j|ddd|jt |ddd|j|ddd |jdddd}|j|ddd|jt |ddd|j|ddd tdd}|jdddd}|j|ddd|j|ddd|j|ddd |jdddd}|j|ddd|j|ddd|j|ddd |jddd d}|j|ddd|j|ddd|j|ddd |jdddd}|j|ddd|jt |ddd|j|ddd y)zCheck max(-0.0,+0.0) issues 2015-08-22 (triggered by bugs in Octave -- sind(-0.0) = +0.0 -- and in some version of Visual Studio -- fmod(-0.0, 360.0) = +0.0.rrr gV@r6r r r i 0r7pf@gʡEK@g/$_@i>0rBf@i[<1r_im/r@i|1rCgO贁Nkg:X0@@gz(_b@io2i1N)rrrrrO)rr!rFs r"test_GeodSolve33zGeodSolveTest.test_GeodSolve33#sb .. Aq# .C3v;&A3v;&A3u:x= .. Aq% 0C3v;&A3v; 6B3u:x= .. Aq# .C3v;@3s6{+YG3u:x= .. Aq# .C3v;@3s6{+YG3u:x= E1 D ,,q!Q $C3v;&A3v;&A3u:x= ,,q!Q $C3v;@3s6{+YG3u:x= ,,q!Q $C3v;@3s6{+YG3u:x= E8 $D ,,q!Q $C3v;&A3v;&A3u:x= ,,q!Q $C3v;&A3v;&A3u:x= ,,q!S# &C3v;&A3v; 6B3u:x= ,,q!Q $C3v;@3s6{+YG3u:x=r$ctjjtjddd}|j tj |d|j tj |d|j tj |dtjjtjddd}|j tj |d|j tj |d|j tj |dy)zoCheck fix for nan + point on equator or pole not returning all nans in Geodesic::Inverse, found 2015-09-23.rrNr r r Nr`r9s r"test_GeodSolve55zGeodSolveTest.test_GeodSolve55Vs .. 1a 4COODJJs6{+,OODJJs6{+,OODJJs5z*+ .. 1b! 4COODJJs6{+,OODJJs6{+,OODJJs5z*+r$ctjjdddd}|j|ddd|j|d d d|j|d d d y)z>Check for points close with longitudes close to 180 deg apart.+=r|rBr g(#=gt=r r gf@r gr~qAr~Nr8r9s r"test_GeodSolve59zGeodSolveTest.test_GeodSolve59bsk .. $4b# >C3v;(97K3v;(:GL3u:'94Hr$cftjjddddtjtjz}|j |ddd|j |d d d|j t |d d dtjjddd d}|jdtjtjz}|j |ddd|j |d d d|j t |d d dy)z0Make sure small negative azimuths are west-going-rgeKrMr g,g~5F@r6r rLr rBPN) rrr&rrrrrO InverseLinertrus r"test_GeodSolve61zGeodSolveTest.test_GeodSolve61is ^^ " "2q*?#+#4#4x7K7K#KMF6&>8VD6&>4@3vf~.VD >> % %b!R1F GD ]]3 1 1H4H4H H IF6&>8VD6&>4@3vf~.VDr$ctjjddddtj}|j dtjtj z}|j |ddd |j |d d d |j t|d d d |j |ddd |j |dd d |j t|dd d |j |ddd |j |ddd |j |ddd |j |ddd |j |ddd |j |dd d |j d!tjtj z}|j |ddd |j |d d d |j t|d d d |j |dd"d |j |dd#d |j |dd d |j |dd$d |j |dd%d |j |dd&d |j |dd'd |j |dd(d |j |dd)d y*)+zCheck for bug in east-going check in GeodesicLine (needed to check for sign of 0) and sign error in area calculation due to a bogus override of the code for alp12. Found/fixed on 2015-12-19.rLgC]r2irBrMrg>@r6r rr rr ggNrr r ir7rgj@+0V@riarg.AUr}rgAV?gHj>rrrog7ُ >fi-1g~f@iFgUg/1瀻nBN)rrrrrtrrrOrrvr(s r"test_GeodSolve65zGeodSolveTest.test_GeodSolve65vs >> % %b*?c&.ll 4D ]]3 x/C/C C DF6&>8vF6&>8vF3vf~. 6J6&>9vF6&>8vF3vf~. 6J6%=8sC6%=8vF6%=7sC6%=:uE6%=9vF6%=1sC ]]3 x/C/C C DF6&>8vF6&>8vF3vf~. 6J6&>9vF6&>:vF6&>8vF6%=8sC6%=9vF6%=5sC6%=:uE6%=:vF6%=+> % %b*7vF6&>:vF6&>8vF ]]3>#,,x/C/CCEF6&>9vF6&>8vF6&>:vFr$cNtjjdddd}|jd|jztj tj z}|j|ddd |j|d d d |j|d d d y)zCheck that DirectLine sets s13.r_rgrrMr7r gQ>@r6r rggDioB@r gIQ,K@N)rr DirectLinertrrrrrrs r"test_GeodSolve71zGeodSolveTest.test_GeodSolve71s >> $ $Q2s 3D ]]3>#,,x/C/CCEF6&>8VD6&>8VD6&>8VDr$c(tjjdddd}|j|ddd|j|d d d|j|d d d|j t j d |d d kDy)a0Check for backwards from the pole bug reported by Anon on 2016-02-13. This only affected the Java implementation. It was introduced in Java version 1.44 and fixed in 1.46-SNAPSHOT on 2016-01-17. Also the + sign on azi2 is a check on the normalizing of azimuths (converting -0.0 to +0.0).rNr|rBg.r gnBT@r6r rrpr rr_N)rrr&rrcracopysignr=s r"test_GeodSolve73zGeodSolveTest.test_GeodSolve73s ^^ " "2r3 5F6&>8VD6&>4@6&>1f=OODMM!VF^4q89r$ctjjddddtj}|j |ddd|j |d d d|j |d d d|j |d dd|j |dd d|j |ddd|j |ddd|j |dddy)zNCheck fix for inaccurate areas, bug introduced in v1.46, fixed 2015-10-16.gΪVK@g?.@gec]K@gQI.@r gOK@r~r r gZ"a)K@r gK;:C@rg1f7L7?rrg$P?rrgMMZArIN)rrrrrr9s r"test_GeodSolve74zGeodSolveTest.test_GeodSolve74s .. '7G!) /C3v; dC3v; dC3u: dC3u:dC3u: dC3u:dC3u:dC3u:Er$ctjjdddd}|j|ddd|j|d d d|j|d d d y)zNThe distance from Wellington and Salamanca (a classic failure of Vincenty)gDg"""""e@g{D@gr g( d@dy=r r g]E!3@r gT- sAƠ>Nr8r9s r"test_GeodSolve76zGeodSolveTest.test_GeodSolve76sl ..  !+\ ;C3v;I3v;I3u:Hr$ctjjdddd}|j|ddd|j|d d d|j|d d d y)z5An example where the NGS calculator fails to convergeg333333;@rg;rr g_YF@rr r g۪I`@r g@, sArNr8r9s r"test_GeodSolve78zGeodSolveTest.test_GeodSolve78sg .. sE5 9C3v;I3v;I3u:Hr$ctjjddddtj}|j |ddd|j |dddtjjddddtj}|j |dd d|j |dd dtjjd dd dtj }|j |d dd |j |d dd|j |ddd |j |ddd |j |ddd|j |dd d|j |dd d|j |ddd|j tjd |d dkD|j tjd |d dkD|j tjd |ddkDtjjddddtj }|j |d dd |j |d dd|j |ddd |j |ddd |j |ddd|j |dd d|j |dd d|j |dddtjjd ddtj}|jdtj}|j tj|d y)zSome tests to add code coverage: computing scale in special cases + zero length geodesic (includes GeodSolve80 - GeodSolve83) + using an incapable line.rrNrg3Kug|=r rgư>r_gK7A4@rrr rr rBr rrrrrr7rgN)rrr GEODESICSCALErrrcrarrsLATITUDErtEMPTYrd)rr!rvr(s r"test_GeodSolve80zGeodSolveTest.test_GeodSolve80s, .. Aq"h.D.D EC3u:~wG3u:~wG .. AtT83I3I JC3u:q':3u:q': .. FAx|| DC3u:q%83u:q$73v;U;3v;U;3u:q483u:q%83u:q%83u:q%8OODMM!SZ0145OODMM!SZ0145OODMM!SZ0145 .. QC >C3u:q%83u:q$73v;593v;U;3u:q$73u:q%83u:q%83u:'8#F >>  q!R):): ;D ]]4 0FOODJJve}-.r$c tjjdddtj}|j tj |d|j tj |d|j tj |dtjjdddtj}|j tj |d|j tj |d|j tj |dtjjddtjd}|j tj |d|j tj |d|j tj |dtjjddtjd}|j tj |d|j tj |d|j tj |dtjjdtjdd}|j |ddk(|j tj |d|j |ddk(tjjdtjdd}|j |ddk(|j tj |d|j |ddk(tjjtjddd}|j tj |d|j tj |d|j tj |dtjjtjddd}|j tj |d|j tj |d|j tj |dy) zTests for python implementation to check fix for range errors with {fmod,sin,cos}(inf) (includes GeodSolve84 - GeodSolve91).rrNr rr rrN)rrr&rainfrcrdrbr=s r"test_GeodSolve84zGeodSolveTest.test_GeodSolve84s<^^ " "1aTXX 6FOODJJvf~./OODJJvf~./OODJJvf~./ ^^ " "1aTXX 6FOODJJvf~./OODJJvf~./OODJJvf~./ ^^ " "1a4 8FOODJJvf~./OODJJvf~./OODJJvf~./ ^^ " "1a4 8FOODJJvf~./OODJJvf~./OODJJvf~./ ^^ " "1dhhD 9FOOF6Na'(OODJJvf~./OOF6Nb() ^^ " "1dhhD 9FOOF6Na'(OODJJvf~./OOF6Nb() ^^ " "488QD 9FOODJJvf~./OODJJvf~./OODJJvf~./ ^^ " "488QD 9FOODJJvf~./OODJJvf~./OODJJvf~./r$ctjjdddd}|j|ddd|j|d d d|j|d d d y)zCheck fix for inaccurate hypot with python 3.[89]. Problem reported by agdhruv https://github.com/geopy/geopy/issues/466 ; see https://bugs.python.org/issue43088gB@gm^gB@g6Ea^r g8V@gHz>r r g*rV@r gL7A`?rINr8r9s r"test_GeodSolve92zGeodSolveTest.test_GeodSolve92$sn .. !3Z!)[ BC3v; TD3v; TD3u:VDr$cNtjjddtjd}|j tj |d|j tj |d|j tj |dy)zPCheck fix for lat2 = nan being treated as lat2 = 0 (bug found 2021-07-26)rrNr r r Nr`r9s r"test_GeodSolve94zGeodSolveTest.test_GeodSolve94.sm .. Atxx 4COODJJs6{+,OODJJs6{+,OODJJs5z*+r$ctdd}|jddddtj}|j|ddd y ) zFailure with long doubles found with test case from Nowak + Nowak Da Costa (2022). Problem was using somg12 > 1 as a test that it needed to be set when roundoff could result in somg12 slightly bigger that 1. Found + fixed 2022-03-30.iRag!6Zwk?rgr\ N@g-YvV@rlWrYr7r NryrEs r"test_GeodSolve96zGeodSolveTest.test_GeodSolve966sE G_ -D ,,q!-/? OC3u:~sCr$N)$r-r.r/r0r:r>rGrJrPrRrTrVrXr]rerlrwrzrrrrrrrrrrrrrrrrrrr1r$r"r3r3s=F >> @ EE EE<,><*C G ?1>f ,I EKB GE : FII&/P"0HE,Dr$r3ceZdZdZej j dZej j dZe dZ e dZ dZ dZ dZd Zd Zd Zd Zd ZdZdZy)PlanimeterTestzPlanimeter testsFTctjj|D](}tjj|d|d*tjj ddS)zHelper function for polygonsrr_FT)rpolygonClearAddPointComputepointsps r" PlanimeterzPlanimeterTest.PlanimeterEs\  " 2%%adAaD12  ! ! ) )% 66r$ctjj|D](}tjj|d|d*tjj ddS)zHelper function for polylinesrr_FT)rpolylinerrrrs r" PolyLengthzPlanimeterTest.PolyLengthMs\!!# 3&&qtQqT23  " " * *5$ 77r$cddgddgddgddgg}tj|\}}}|j|dd|j|d d d dgd dgd dgd dgg}tj|\}}}|j|dd|j|d d dd gd dgdd gd dgg}tj|\}}}|j|dd|j|dd ddgddgddgg}tj|\}}}|j|dd |j|dd tj|\}}}|j|dd |j t j |y)z2Check fix for pole-encircling bug found 2011-03-16YrrNrBigvH#Ag-C6?r g8=Br_ig8=gaӋ'#Ag(BignBi}1N)rrrrrcrardrr_ perimeterareas r"test_Planimeter0zPlanimeterTest.test_Planimeter0Us1gBx"cRI 6F'226:Ay$9k4@4:Ahb C:Sz :F'226:Ay$9k4@4;"gAwAA /F'226:Ay$9k4@4:1g1v2w 'F'226:Ay$9h:4!11='226:Ay$9h:OODJJt$%r$cddgddgddgg}tj|\}}}|j|dd|j|ddy ) z;Check fix for Planimeter pole crossing bug found 2011-06-24rrgfffffV@g|fi:r_r g4=BNrrrrs r"test_Planimeter5zPlanimeterTest.test_Planimeter5ms[3i"db&\ 2F'226:Ay$9fa84:r$cdddgddgddgg}tj|\}}}|j|dd|j|ddddgddgddgg}tj|\}}}|j|dd|j|ddddgddgddgg}tj|\}}}|j|dd|j|ddddgddgddgg}tj|\}}}|j|dd|j|ddy ) zjaD TAaC5!))11%?jaD TAaC5!))11$=jaD TA2a45!))11$>jaD TA2a4"9c:7;r$ctjjtjjddtjj ddtjj ddtjj ddtjj dd\}}}|j |dd y ) zACheck fix to transitdirect vs transit zero handling inconsistencyrrNriFTg.Ar\r N)rrrrAddEdgerr)rrrs r"test_Planimeter29z PlanimeterTest.test_Planimeter29s  "##Aq)""B-""Q-""3-''//trsQ +D>8$$D>LoDH%%oDb [:X&&[:r$