def parallel_check(vec1, vec2):
    from ..const import PRM
    import numpy as np
    par = False
    for n,v in enumerate([vec1, vec2]):
        if not len(v.shape) == 1:
            raise ValueError("Bad shape for vector
    if not vec1.shape[0] == vec2.shape[0]:
        raise ValueError("Vector length mismatch")
    angle = vec_angle(vec1, vec2)
    if min([abs(angle), abs(angle - 180.)]) < PRM.NON_PARALLEL_TOL:
        par = True
    return par