import torch
import numpy as np

positional_encoding = None

def get_positional_encoding(n_positions, n_dimensions, device='cpu'):
    global positional_encoding
    if positional_encoding is None:
        numerators = torch.tensor(range(n_positions)).repeat(n_dimensions, 1).T
        denominators = 10000 ** (torch.tensor(range(n_dimensions)) // 2 * 2 / n_dimensions)
        print('denoms:', denominators)
        positional_encoding = numerators / denominators
        positional_encoding[:, ::2] = torch.sin(positional_encoding[:, ::2])
        positional_encoding[:, 1::2] = torch.cos(positional_encoding[:, 1::2])
    # output shape: (seqlen, hiddendim)
    return torch.tensor(positional_encoding, dtype=torch.float, device=device)