def printRho(SGLaue, rho, rhoMax):
    """default doc string
    
    :param type name: description
    
    :returns: type name: description
    
    """
    dim = len(rho.shape)
    if dim == 2:
        ix, jy = rho.shape
        for j in range(jy):
            line = ''
            if SGLaue in ['3', '3m1', '31m', '6/m', '6/mmm']:
                line += (jy - j) * '  '
            for i in range(ix):
                r = int(100 * rho[i, j] / rhoMax)
                line += '%4d' % r
            print(line + '\n')
    else:
        ix, jy, kz = rho.shape
        for k in range(kz):
            print('k = %d' % k)
            for j in range(jy):
                line = ''
                if SGLaue in ['3', '3m1', '31m', '6/m', '6/mmm']:
                    line += (jy - j) * '  '
                for i in range(ix):
                    r = int(100 * rho[i, j, k] / rhoMax)
                    line += '%4d' % r
                print(line + '\n')
