mò ¼v'Bc@s/dZdZdZyeWnd\ZZnXy|dkZdkZdkZdkZdk Z ei i ei i e ieiƒƒdddgƒƒd klZlZWn1yeiid eƒWn eGHnX‚nXeidiei ƒd Zeid ƒdZeiƒdZd „Zd„Zd„Zd„Zd„Zd„Ze djoÈei!dƒordk"Z"dk#Z#e"i$ƒe#i%ƒZ&e&i'dƒoee(e)ee(e)ƒƒq'eddeddƒƒq+e*eiƒdjoee+eidƒƒGHq+eiidƒndS(s€calculate ipv6 address using Tony Hain's draft provider-independent method refer to draft-hain-ipv6-PI-addr-01.txt Oct. 2001s> hain_address.py -- calculate provider-independent ipv6 address Copyright (C) 2005 John Comeau This program is free software; you can redistribute it and/or modify it under the terms of the GNU General Public License as published by the Free Software Foundation; either version 2 of the License, or (at your option) any later version. This program is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details. You should have received a copy of the GNU General Public License along with this program; if not, write to the Free Software Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA. s5Not all needed libraries found, upgrade or check pathiiNitlibtpython(sgpls jclicenses%s iÿÿÿÿt.cGstS(N(tFalse(twhatever((t./hain_address.pyt DebugPrint1scGs0x)t|dƒtijo|d}qW|S(Ni(ttypetargsttypest TupleType(R((Rtscalars4sc Gstd|ƒyt|ƒ\}}}Wnt|ƒ\}}d}nXtdt|ƒƒ}tdt|ƒƒ}td||ƒ|d}ddd >}t ||d ƒ}t t |ƒ|d ƒ}|djo|dd >O}nt||ƒtd ||fƒd |t||ƒfS(smain routine, calculates address from input lat and long quote from draft: Creating the Provider-Independent address prefix is accomplished by dividing the WGS-84 latitude and longitude readings in degrees (for latitudes between +/- 15 degrees measured to 5 places right of the decimal ie: deg.xxxxx) by the incremental angle. For 22-bit values the incremental angle is 0.0000858306884765625 degrees (360/(2^22)). This will result in uniquely identifiable areas with 2096832 along the latitude axis and 4193664 along the longitude axis. The specific sequence for address formation is: 1. for east/west : convert all values to east (ie: 90w = 270e) 2. divide demarcation degree values by 0.0000858306884765625 3. convert each of the integers to 22-digit binary 4. for north/south (set MSB for south) 5. bit interleave latitude, longitude into 44-bit result 6. prepend FP xxxx (IANA assigned) to form 48-bit prefix (0001 is recommended to avoid breaking up any of the unassigned 3-bit spaces)sargs:t1tlatitudet longitudes lat, lon:if360.0iif0.5is%x %xs%s%011xN(RRR R Rtprefixt canonicalizetstrtlattlont convertertintteastwesttabst northsoutht interleave( RRRRRRRR R((Rthainaddr:s&     cGs$t|ƒ}|\}}t}t|||ƒy…|ddgjotd|‚n^|djo*ti dƒi |ƒp td‚q²n'ti dƒi |ƒp td‚nWn-tj o!}td|id ƒd SnX|djo‹|i ƒi d ƒo|d d !t}}n/|i ƒi d ƒo|d d !t}}nt|ƒ}|o|d jo | }q n¦|i ƒi dƒo|d d !t}}n/|i ƒi dƒo|d d !t}}nt|ƒ}|o|d jod|}n|d jod|}n|S(NR Rscannot grok direction %ss[+-]?\d+(?:\.\d{0,})?[NSns]?$simproper latitude formats[+-]?\d+(?:\.\d{0,})?[EWew]?$simproper longitude formats ValueError: if0.0tsiÿÿÿÿtntwteih(R Rt directiontvalueRtadjustRt ValueErrortretcompiletmatchtinstancetlowertendswithtTruetfloat(RRR&R R!((RRbsB        cGsãt|ƒ}y+t|ddƒt|ddƒ}}Wntj o|\}}nXd}xetddƒD]T}d|}||d|>@djO}|dK}||d|>@djO}|dK}qoW|dL}t d|ƒ|S(sFinterleave binary digits of two numbers according to Tony Hain's draftiiiiis result: %xN( R RRtatbt TypeErrortresulttrangetpositiontshiftR(RR+R,R1R.R0((RR„s" +   cGsFt|ƒ\}}}d|||dd!|dd!|dd!fGHdS(s9use HTML form to make this script available as CGI scripts0Content-type: text/html
Latitude: Example: 30.73N
Longitude: Example: 117.18W
IPv6 Address: /48
iiii N(R RRRtipv6addr(RRR2R((Rtcgiform–st__main__t SERVER_NAMERisUsage: %s LATITUDE LONGITUDE (ii(,t__doc__t Copyrightt errormessageR)RtsystosR R#tpwdtpathtappendtseptjointgetpwuidtgeteuidt com.jcomeautgplt jclicensetstderrtwritetargvtsplittselftcommandt originalselfRR RRRR3t__name__tgetenvtcgitbtcgitenabletparsetcgivarsthas_keyRRtlenttuple(R7RR3RRRCRIR#R)RORRRKR RNR;R R9RR8RRDRJR:((Rt?sH-8    ( "