code: 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