libs5-python 1.0.5__py3-none-any.whl

libs5_python-1.0.5.dist-info/METADATA

@@ -0,0 +1,49 @@
+Metadata-Version: 2.4
+Name: libs5-python
+Version: 1.0.5
+Summary: libs5 python scripts, and helper classes (useg e.g. in BEANS software)
+Project-URL: Homepage, https://gitlab.com/ahypki/libs5-python
+Project-URL: Issues, https://gitlab.com/ahypki/libs5-python/-/issues/
+Author: Arkadiusz Hypki
+License-Expression: MIT
+License-File: LICENSE
+Classifier: Operating System :: OS Independent
+Classifier: Programming Language :: Python :: 3
+Requires-Python: >=3.9
+Description-Content-Type: text/markdown
+
+# libs5-python
+
+`libs5-python` is a package with `Python` scripts, and helper classes which are used e.g. in `BEANS` software
+
+
+## Installation
+
+To install `libs5-python` type in the console:
+
+```
+pip install libs5-python
+```
+
+To upgrade `libs5-python` type in the console:
+
+```
+pip install libs5-python --upgrade
+```
+
+## Roadmap
+
+There is no specific roadmap for `libs5-python`. New features are added if they are needed.
+
+## Contributing
+
+Feel free to contribute to this project by sending me your opinion, patches, or requesting some features through gitlab issue system.
+
+## Links
+
+Home page: https://gitlab.com/ahypki/libs5-python/
+Issues   : https://gitlab.com/ahypki/libs5-python/-/issues/
+
+## License
+
+`libs5-python` is released under the MIT license.

libs5_python-1.0.5.dist-info/RECORD

@@ -0,0 +1,22 @@
+net/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
+net/hypki/File.py,sha256=52wnCUBcDlBGH5swwR-PHJ9ofaW6BYloNgCLR605o9g,1007
+net/hypki/Ini.py,sha256=938r36Pd1hRnwjX5ZjXqZIbRyIlklh1hjZfA1DN8ZPU,3044
+net/hypki/Logger.py,sha256=vvJXYtb5qnLWmSoEnrtM937BWT2DAQUquBm-DN2umo8,1650
+net/hypki/StringUtils.py,sha256=UuUnPIlkIe0WzUo-Q0o8HmE0kQP6MJwNZzK2fm39J-Y,438
+net/hypki/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
+net/hypki/apiutils.py,sha256=E_0EKQ46BbdTHqVABky8DqIe1WetpglmsPWgnjUGeJQ,1036
+net/hypki/args.py,sha256=pXCJs2-LP1w8-3V-6KtOZuw4V9Q3VP7typ54qJ62GqM,516
+net/hypki/asserts.py,sha256=awT2uysmbsPFpOBf6vaS0r39QBvsI5QGrwreXk0bwhE,198
+net/hypki/io.py,sha256=H1ImSJojFbVa8Ww1iQ1Nr-sIjz_hEM0HkhMmSkFZEzo,1878
+net/hypki/log.py,sha256=Qlbm1nZVBgs7EXF3Y12upGOU5Mpg0tL5YhJ3e6HGthQ,1730
+net/hypki/math.py,sha256=-ZcIvUcSD5VTtPer6chvrARjL4BKEGqwsgnmAyfld4Q,197
+net/hypki/random.py,sha256=svaPZV2B5DQj_GToIX49pOKxbgdKYY8hmu743Gyvw_c,231
+net/hypki/regex.py,sha256=t80iy6fXiiQOYB1FIjpgpfrLdG2QOwKtRgnrCAn-9Ws,934
+net/hypki/system.py,sha256=9MOAIJ221BH4GR4KKpPXexJGmFaWGZ1_wgi49kaFpmE,558
+net/hypki/tests/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
+net/hypki/tests/kepler.py,sha256=3mEdHecMvM_ztUO17ZVvhxIXTze12-xx71A6Y3iaf_M,2307
+net/hypki/tests/numpytest.py,sha256=FAyGjE2u_rqcoXUtj8HP4VbWP2f_LhOJ3hYVp37tDsA,378
+libs5_python-1.0.5.dist-info/METADATA,sha256=mU-_CkJ3BqnoPIumvMAz762Y0w91WJG-464PLI73lQ4,1257
+libs5_python-1.0.5.dist-info/WHEEL,sha256=qtCwoSJWgHk21S1Kb4ihdzI2rlJ1ZKaIurTj_ngOhyQ,87
+libs5_python-1.0.5.dist-info/licenses/LICENSE,sha256=jUGenC6Eo6P1BkABG3AQ2ThWGSDHgxa58dRIB8iR8Pc,1072
+libs5_python-1.0.5.dist-info/RECORD,,

libs5_python-1.0.5.dist-info/WHEEL

@@ -0,0 +1,4 @@
+Wheel-Version: 1.0
+Generator: hatchling 1.27.0
+Root-Is-Purelib: true
+Tag: py3-none-any

libs5_python-1.0.5.dist-info/licenses/LICENSE

@@ -0,0 +1,21 @@
+MIT License
+
+Copyright (c) 2025 Arkadiusz Hypki
+
+Permission is hereby granted, free of charge, to any person obtaining a copy
+of this software and associated documentation files (the "Software"), to deal
+in the Software without restriction, including without limitation the rights
+to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
+copies of the Software, and to permit persons to whom the Software is
+furnished to do so, subject to the following conditions:
+
+The above copyright notice and this permission notice shall be included in all
+copies or substantial portions of the Software.
+
+THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
+IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
+FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
+AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
+LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
+OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
+SOFTWARE.

net/hypki/File.py

@@ -0,0 +1,52 @@
+'''
+Created on Oct 14, 2025
+
+@author: ahypki
+'''
+import os
+import shutil
+from pathlib import Path
+
+class File(object):
+    '''
+    classdocs
+    '''
+    __absPath = None
+
+    def __init__(self, path):
+        '''
+        Constructor
+        '''
+        if isinstance(path, str):
+            self.setAbsPath(path)
+        elif len(path) == 1:
+            self.setAbsPath(path[0])
+        else:
+            self.setAbsPath("/".join(path))
+            # p2 = "/"
+            # for p1 in path:
+            #     p1 += 
+        
+    def getAbsPath(self):
+        return self.__absPath
+
+
+    def setAbsPath(self, path):
+        self.__absPath = Path(path).expanduser()# os.path.abspath(path)
+        
+    def exists(self):
+        return os.path.exists(self.__absPath)
+    
+    def mkdirs(self):
+        os.makedirs(self.__absPath, exist_ok=True)
+
+    
+    def copy(self, destination):
+        shutil.copy2(self.__absPath, destination)
+
+    
+    def remove(self):
+        os.remove(self.__absPath)
+    
+    
+    

net/hypki/Ini.py

@@ -0,0 +1,102 @@
+'''
+Created on Oct 15, 2025
+
+@author: ahypki
+'''
+import re
+
+class Ini(object):
+    '''
+    classdocs
+    '''
+    __path = None
+    __iniLines = None
+
+    def __init__(self, path):
+        '''
+        Constructor
+        '''
+        self.__path = path
+        
+    def __getIni(self):
+        if self.__iniLines is None:
+            self.__iniLines = open(self.__path, "r").readlines()
+                
+        return self.__iniLines
+    
+    def __isKey(self, line):
+        k = self.__getKey(line)
+        return False if k is None else True
+    
+    def __getKey(self, line):
+        m = re.match("^[\s]*(\w[\w\d]*)[\s]*=[\s]*.*$", line)
+        if m:
+            return m.group(1)
+        return None
+    
+    def __isSection(self, line):
+        sec = self.__getSection(line)
+        return False if sec is None else True
+    
+    def __getSection(self, line):
+        m = re.match("^[\s]*\[(\w[\w\d]+)\][\s]*$", line)
+        if m:
+            return m.group(1)
+        return None
+        
+    def getSections(self):
+        sections = []
+        for line in self.__getIni():
+            if self.__isSection(line):
+                sections.append(self.__getSection(line))
+        return sections
+    
+    def getKeys(self, section):
+        keys = []
+        currentSection = None
+        for line in self.__getIni():
+            if self.__isSection(line):
+                currentSection = self.__getSection(line)
+                continue
+            if self.__isKey(line):
+                if currentSection is not None and currentSection == section:
+                    keys.append(self.__getKey(line))
+        return keys
+    
+    def set(self, section, key, value):
+        newLines = []
+        currentSection = None
+        for line in self.__getIni():
+            lineAdded = False
+            
+            if self.__isSection(line):
+                currentSection = self.__getSection(line)
+                
+            if self.__isKey(line):
+                currentKey = self.__getKey(line)
+                if currentSection is not None and currentSection == section and currentKey == key:
+                    newLines.append(key + " = " + value)
+                    lineAdded = True
+                    
+            if not lineAdded:
+                newLines.append(line)
+                
+        self.__iniLines = newLines
+        
+    def get(self, section, key):
+        currentSection = None
+        for line in self.__getIni():
+            if self.__isSection(line):
+                currentSection = self.__getSection(line)
+                continue
+            if self.__isKey(line):
+                currentKey = self.__getKey(line)
+                if currentSection is not None and currentSection == section and currentKey == key:
+                    return line[(line.index("=") + 1):].strip()
+        return None
+        
+    def save(self):
+        if self.__iniLines is not None:
+            with open(self.__path, 'w') as f:
+                for line in self.__iniLines:
+                    f.write(f"{line.strip()}\n")

net/hypki/Logger.py

@@ -0,0 +1,51 @@
+'''
+Created on Aug 13, 2025
+
+@author: ahypki
+'''
+from rich.console import Console
+
+console = Console()
+
+class Logger:
+    BLUE = '[blue]'#'\033[94m'
+    BLUE_END = '[/blue]'
+    GREEN = '[green]'#'\033[92m'
+    GREEN_END = '[/green]'
+    RED = '[red]'#'\033[91m'
+    RED_END = '[/red]'
+    WARNING = '[yellow]'#'\033[93m'
+    WARNING_END = '[/yellow]'
+    FAIL = RED#'\033[91m'
+    FAIL_END = RED_END
+    # ENDC = '\033[0m'
+    
+    DEBUG = False
+    INFO = True
+    
+    PRINT_LEVEL = True
+    
+    # def blue(s):  # @NoSelf
+    #     return Logger.BLUE + s + Logger.BLUE_END
+    
+    def debug(msg, printLogLevel = True, printNewLine = True):  # @NoSelf
+        if Logger.DEBUG:
+            console.print(("DEBUG " if printLogLevel and Logger.PRINT_LEVEL else '') 
+                  + str(msg), 
+                  end = ('\n' if printNewLine else ''))
+    
+    def info(msg, printLogLevel = True, printNewLine = True):  # @NoSelf
+        if Logger.INFO:
+            console.print(("INFO  " if printLogLevel and Logger.PRINT_LEVEL else '') 
+                  + str(msg), 
+                  end = ('\n' if printNewLine else ''))
+        
+    def warn(msg, printLogLevel = True, printNewLine = True):  # @NoSelf
+        console.print(("WARN  " if printLogLevel and Logger.PRINT_LEVEL else '') 
+              + str(Logger.WARNING + msg + Logger.WARNING_END), 
+              end = ('\n' if printNewLine else ''))
+    
+    def error(msg, printLogLevel = True, printNewLine = True):  # @NoSelf
+        console.print(("ERROR " if printLogLevel and Logger.PRINT_LEVEL else '') 
+              + str(msg), 
+              end = ('\n' if printNewLine else ''))

net/hypki/StringUtils.py

@@ -0,0 +1,23 @@
+'''
+Created on Oct 14, 2025
+
+@author: ahypki
+'''
+
+class StringUtils(object):
+    '''
+    classdocs
+    '''
+
+
+    def __init__(self, params):
+        '''
+        Constructor
+        '''
+        pass
+    
+    def notEmpty(s):  # @NoSelf
+        return not StringUtils.isEmpty(s) #len(s.strip()) > 0 if s is not None else False 
+    
+    def isEmpty(s):  # @NoSelf
+        return s is None or len(s.strip()) == 0 if s is not None else False 

net/hypki/apiutils.py

@@ -0,0 +1,32 @@
+'''
+Created on Aug 15, 2025
+
+@author: ahypki
+'''
+import requests
+from requests_toolbelt.multipart.encoder import MultipartEncoder
+
+def streamGET(url, params):
+    s = requests.Session()
+    return s.get(url, params=params, headers=None, stream=True).iter_lines()
+#    with s.get(url, params=params, headers=None, stream=True) as resp:
+#        for line in resp.iter_lines():
+#            if line:
+#                print(line)
+
+def streamPOST(url, file, params):
+    # s = requests.Session()
+    # Specify the content type
+    # headers = {'Content-Type': 'application/octet-stream'}
+    # headers = {'Content-Type': 'multipart/form-data'}
+    # return s.post(url, headers=headers, data=stream)
+    fields = params
+    fields['file'] = (file, open(file, 'rb'), 'text/plain')
+    multipart_data = MultipartEncoder(
+        fields=fields
+    )
+    return requests.post(url,
+                         data=multipart_data,
+                         headers={'Content-Type': multipart_data.content_type}, 
+                         timeout=30)
+

net/hypki/args.py

@@ -0,0 +1,21 @@
+'''
+Created on 15-06-2013
+
+@author: ahypki
+'''
+
+import sys
+
+class ArgUtils:
+
+    def getString(name, defaultvalue):  # @NoSelf
+        for i in range(len(sys.argv)):
+            if sys.argv[i] == "--" + name or sys.argv[i] == "-" + name:
+                return sys.argv[i + 1].strip()
+        return defaultvalue
+    
+    def isArgPresent(name):  # @NoSelf
+        for i in range(len(sys.argv)):
+            if sys.argv[i] == "--" + name or sys.argv[i] == "-" + name:
+                return True
+        return False

net/hypki/asserts.py

@@ -0,0 +1,12 @@
+'''
+Created on 15-06-2013
+
+@author: ahypki
+'''
+
+import sys
+from builtins import Exception
+
+def assertTrue(condition, msg):
+    if condition is False:
+        raise Exception("Assertion failed", msg)

net/hypki/io.py

@@ -0,0 +1,71 @@
+'''
+Created on 08-05-2013
+
+@author: ahypki
+'''
+
+import os.path
+import re
+import subprocess
+import sys
+import time
+import urllib.request
+from net.hypki.regex import *
+from net.hypki.log import *
+
+def absolutePath(path):
+    return os.__path.abspath(path)
+
+def findFolders(path):
+    dirsAll = []
+    for dirName, dirs, files in os.walk(path):
+        dirsAll.append(dirName)
+#         path = root.split('/')
+#         print (len(path) - 1) *'---' , os.path.basename(root)       
+#         for file in files:
+#             print len(path)*'---', file
+    return dirsAll
+
+def readString( prompt ):
+    print(bcolors.OKBLUE + prompt + bcolors.ENDC, end = '', flush = True)
+    return sys.stdin.readline().strip()
+
+def linux_cmd(cmd, workingDir = None, printImmediatelly = False):
+    if printImmediatelly:
+        debug("linux cmd: " + cmd);
+    result = []
+    s = subprocess.Popen(cmd, cwd = workingDir, shell=True, stderr=subprocess.STDOUT, stdout=subprocess.PIPE)
+    for line in s.stdout.readlines():
+        strLine = str(line).rstrip()
+#         if printImmediatelly is False:
+        result.append(firstGroup(strLine, 'b\'(.*)\\\\n\''))
+#         else:
+        if printImmediatelly:
+            print(firstGroup(strLine, 'b\'(.*)\\\\n\''))
+    return result
+
+def linux_cmd_fork(cmd):
+#     info("linux cmd fork: " + cmd);
+    os.system(cmd)
+
+def is_file_exists(path):
+    return os.path.exists(path)
+
+def url_read(url):
+    resp = urllib.request.urlopen(url)
+    return resp.read().decode('utf-8')
+
+def read_file(path):
+    with open(path, 'r') as content_file:
+        return content_file.read()
+    
+def append_to_file(path, line):
+    with open(path, 'a') as content_file:
+        content_file.write('\n')
+        content_file.write(line)
+    
+def file_save(path, content):
+    myFile = open(path, 'w', encoding = 'UTF8')
+    myFile.write(content)
+    myFile.close()
+    

net/hypki/log.py

@@ -0,0 +1,76 @@
+'''
+Created on 08-05-2013
+
+@author: ahypki
+'''
+
+import os.path
+import re
+import subprocess
+import sys
+import time
+import urllib.request
+
+TRACE = False
+DEBUG = True
+INFO = True
+
+class bcolors:
+    HEADER = '\033[95m'
+    OKBLUE = '\033[94m'
+    OKGREEN = '\033[92m'
+    RED = '\033[91m'
+    WARNING = '\033[93m'
+    FAIL = '\033[91m'
+    ENDC = '\033[0m'
+
+    def disable(self):
+        self.HEADER = ''
+        self.OKBLUE = ''
+        self.OKGREEN = ''
+        self.RED = ''
+        self.WARNING = ''
+        self.FAIL = ''
+        self.ENDC = ''
+        
+def green( msg ):
+    print(bcolors.OKGREEN + str(msg) + bcolors.ENDC)
+    return
+
+def blue( msg ):
+    print(bcolors.OKBLUE + msg + bcolors.ENDC)
+    return
+
+def red( msg ):
+    print(bcolors.RED + msg + bcolors.ENDC)
+    return
+
+def msg( msg ):
+    print(bcolors.OKBLUE + msg + bcolors.ENDC)
+    return
+
+def trace( parameters ):
+    if (TRACE):
+        print("TRACE  " + parameters)
+    return
+
+def debug(parameters, printLogLevel = True, printNewLine = True):
+    if (DEBUG):
+        print(("DEBUG " if printLogLevel else '') + str(parameters), end = ('\n' if printNewLine else ''))
+    return
+
+def log( parameters ):
+    debug(parameters)
+
+def info( parameters, printLogLevel = True, printNewLine = True):
+    if (INFO):
+        print(("INFO  " if printLogLevel else '') + str(parameters), end = ('\n' if printNewLine else ''))
+    return
+
+def warn(parameters, printLogLevel = True, printNewLine = True):
+    print(bcolors.WARNING + "WARN   " + str(parameters) + bcolors.ENDC)
+    return
+
+def error( parameters, printLogLevel = True, printNewLine = True, end = '' ):
+    print(bcolors.FAIL + "ERROR   " + str(parameters) + bcolors.ENDC, end, file=sys.stderr)
+    return

net/hypki/math.py

@@ -0,0 +1,13 @@
+
+import subprocess
+import sys
+import re
+import os
+from net.hypki.io import *
+from net.hypki.args import *
+
+def math_eval(s):
+    try:
+        return eval(s)
+    except SyntaxError:
+        return 0

net/hypki/random.py

@@ -0,0 +1,14 @@
+'''
+Created on 18-07-2013
+
+@author: ahypki
+'''
+
+import os
+import uuid
+import base64
+
+def nextUuid():
+    return str(uuid.uuid4())
+#     r_uuid = base64.urlsafe_b64encode(uuid.uuid4().bytes)
+#     return str(r_uuid.replace('=', ''))

net/hypki/regex.py

@@ -0,0 +1,50 @@
+'''
+Created on 08-05-2013
+
+@author: ahypki
+'''
+
+from net.hypki import *
+import os.path
+import re
+import subprocess
+import sys
+import time
+# import urllib.request
+
+def replace(string, searchPattern, replacement):
+    return re.sub(searchPattern, replacement, string, count=0, flags=0)
+
+def firstGroup(s, regex):
+    r = re.compile(regex)
+    m = r.search(s)
+    if m:
+        return m.group(1)
+    return ""
+
+def allGroups(s, regex):
+    r = re.compile(regex)
+    m = r.search(s)
+    if m:
+        return m.groups()
+    return ""
+
+def isLong(strNumber):
+    return matches(strNumber, "^[\d]+$")
+    
+def isDouble(strNumber):
+    return matches(strNumber, "[\d]+\.[\d\+\-eE]+")
+
+def matches(s, regex):
+    r = re.compile(regex)
+    m = r.search(s)
+    if m:
+        return True
+    return False
+
+def matchesIgnoreCase(s, regex):
+    r = re.compile(regex, re.IGNORECASE)
+    m = r.search(s)
+    if m:
+        return True
+    return False

net/hypki/system.py

@@ -0,0 +1,28 @@
+'''
+Created on 18-07-2013
+
+@author: ahypki
+'''
+
+import inspect
+import os.path
+import re
+import subprocess
+import sys
+import time
+import urllib.request
+
+
+def getThisScriptPath():
+    return inspect.stack()[1][1] # inspect.getfile(inspect.currentframe()) # script filename (usually with path)
+
+
+def getThisScriptDirectory():
+    return os.__path.dirname(os.__path.abspath(inspect.stack()[1][1])) # script directory
+
+
+def getSystemVariable(name):
+    for param in os.environ.keys():
+        if param == name:
+            return os.environ[param]
+    return None

net/hypki/tests/kepler.py

@@ -0,0 +1,91 @@
+from astropy.constants import G, M_earth, R_earth
+from astropy import units as u
+import numpy as np
+
+mu = G.value*M_earth.value
+Re = R_earth.value
+
+#Test vectors
+r_test = np.array([Re + 600.0*1000, 0, 50])
+v_test = np.array([0, 6.5 * 1000, 0])
+t = 0
+
+
+def cart_2_kep(r_vec,v_vec):
+    #1
+    h_bar = np.cross(r_vec,v_vec)
+    h = np.linalg.norm(h_bar)
+    #2
+    r = np.linalg.norm(r_vec)
+    v = np.linalg.norm(v_vec)
+    #3
+    E = 0.5*(v**2) - mu/r
+    #4
+    a = -mu/(2*E)
+    #5
+    e = np.sqrt(1 - (h**2)/(a*mu))
+    #6
+    i = np.arccos(h_bar[2]/h)
+    #7
+    omega_LAN = np.arctan2(h_bar[0],-h_bar[1])
+    #8
+    #beware of division by zero here
+    lat = np.arctan2(np.divide(r_vec[2],(np.sin(i))),\
+    (r_vec[0]*np.cos(omega_LAN) + r_vec[1]*np.sin(omega_LAN)))
+    #9
+    p = a*(1-e**2)
+    nu = np.arctan2(np.sqrt(p/mu) * np.dot(r_vec,v_vec), p-r)
+    #10
+    omega_AP = lat - nu
+    #11
+    EA = 2*np.arctan(np.sqrt((1-e)/(1+e)) * np.tan(nu/2))
+    #12
+    n = np.sqrt(mu/(a**3))
+    T = t - (1/n)*(EA - e*np.sin(EA))
+
+    return a,e,i,omega_AP,omega_LAN,T, EA
+
+def kep_2_cart(a,e,i,omega_AP,omega_LAN,T, EA):
+
+    #1
+    n = np.sqrt(mu/(a**3))
+    M = n*(t - T)
+    #2
+    MA = EA - e*np.sin(EA)
+    #3
+    #
+    # ERROR WAS HERE
+   #nu = 2*np.arctan(np.sqrt((1-e)/(1+e)) * np.tan(EA/2))
+    nu = 2*np.arctan(np.sqrt((1+e)/(1-e)) * np.tan(EA/2))
+    #4
+    r = a*(1 - e*np.cos(EA))
+    #5
+    h = np.sqrt(mu*a * (1 - e**2))
+    #6
+    Om = omega_LAN
+    w =  omega_AP
+
+    X = r*(np.cos(Om)*np.cos(w+nu) - np.sin(Om)*np.sin(w+nu)*np.cos(i))
+    Y = r*(np.sin(Om)*np.cos(w+nu) + np.cos(Om)*np.sin(w+nu)*np.cos(i))
+    Z = r*(np.sin(i)*np.sin(w+nu))
+
+    #7
+    p = a*(1-e**2)
+
+    V_X = (X*h*e/(r*p))*np.sin(nu) - (h/r)*(np.cos(Om)*np.sin(w+nu) + \
+    np.sin(Om)*np.cos(w+nu)*np.cos(i))
+    V_Y = (Y*h*e/(r*p))*np.sin(nu) - (h/r)*(np.sin(Om)*np.sin(w+nu) - \
+    np.cos(Om)*np.cos(w+nu)*np.cos(i))
+    V_Z = (Z*h*e/(r*p))*np.sin(nu) + (h/r)*(np.cos(w+nu)*np.sin(i))
+
+    return [X,Y,Z],[V_X,V_Y,V_Z]
+
+a,e,i,omega_AP,omega_LAN,T, EA = cart_2_kep(r_test,v_test)
+r_test2, v_test2 = kep_2_cart(a,e,i,omega_AP,omega_LAN,T, EA)
+
+
+# print(r_test2 - r_test)
+# print(v_test2 - v_test)
+for x in range(100):
+    r_test2, v_test2 = kep_2_cart(a,e,i,omega_AP+x,omega_LAN,T, EA)
+    print(r_test2, v_test2)

net/hypki/tests/numpytest.py

@@ -0,0 +1,34 @@
+'''
+Created on 21-06-2013
+
+@author: ahypki
+'''
+
+import numpy as np
+import scipy.signal
+
+a = np.zeros(100)
+a[0] = 1
+a[1] = 2
+a[2] = 3
+a[3] = 4
+a[4] = 3
+a[5] = 1
+a[6] = 0
+a[7] = -1
+a[8] = 2
+a[9] = 3
+a[10] = 2
+a[11] = 1
+print(a)
+res = scipy.signal.argrelmax(a)
+print(res)
+
+b = np.zeros(100)
+b[0] = 1
+b[1] = 2
+b[2] = 1
+b[3] = 2
+bb = np.trim_zeros(b)
+print(bb.mean())
+print(bb.std())