mirror of
https://github.com/vladmandic/automatic
synced 2026-08-28 08:00:59 +02:00
Disty0
222 lines
10 KiB
Python
222 lines
10 KiB
Python
# pylint: disable=redefined-builtin
|
|
|
|
import math
|
|
import torch
|
|
|
|
from modules import devices
|
|
from .common import dtype_dict, use_contiguous_mm, conv_types, conv_transpose_types
|
|
|
|
|
|
@devices.inference_context()
|
|
def get_scale_asymmetric(weight: torch.FloatTensor, reduction_axes: int | list[int], weights_dtype: str) -> tuple[torch.FloatTensor, torch.FloatTensor]:
|
|
zero_point = torch.amin(weight, dim=reduction_axes, keepdims=True)
|
|
scale = torch.amax(weight, dim=reduction_axes, keepdims=True).sub_(zero_point).div_(dtype_dict[weights_dtype]["max"] - dtype_dict[weights_dtype]["min"])
|
|
if dtype_dict[weights_dtype]["min"] != 0:
|
|
zero_point.sub_(torch.mul(scale, dtype_dict[weights_dtype]["min"]))
|
|
return scale, zero_point
|
|
|
|
|
|
@devices.inference_context()
|
|
def get_scale_symmetric(weight: torch.FloatTensor, reduction_axes: int | list[int], weights_dtype: str) -> torch.FloatTensor:
|
|
return torch.amax(weight.abs(), dim=reduction_axes, keepdims=True).div_(dtype_dict[weights_dtype]["max"])
|
|
|
|
|
|
@devices.inference_context()
|
|
def quantize_weight(weight: torch.FloatTensor, reduction_axes: int | list[int], weights_dtype: str, dtype: torch.dtype = None, use_stochastic_rounding: bool = False) -> tuple[torch.Tensor, torch.FloatTensor, torch.FloatTensor]:
|
|
if weight.dtype != torch.float64:
|
|
weight = weight.to(dtype=torch.float32, copy=False)
|
|
|
|
if dtype_dict[weights_dtype]["is_unsigned"]:
|
|
scale, zero_point = get_scale_asymmetric(weight, reduction_axes, weights_dtype)
|
|
if dtype is not None:
|
|
scale = scale.to(dtype=dtype)
|
|
zero_point = zero_point.to(dtype=dtype)
|
|
quantized_weight = torch.sub(weight, zero_point).div_(scale)
|
|
else:
|
|
scale = get_scale_symmetric(weight, reduction_axes, weights_dtype)
|
|
zero_point = None
|
|
if dtype is not None:
|
|
scale = scale.to(dtype=dtype)
|
|
quantized_weight = torch.div(weight, scale)
|
|
|
|
if dtype_dict[weights_dtype]["is_integer"]:
|
|
if use_stochastic_rounding:
|
|
quantized_weight.add_(torch.randn_like(quantized_weight), alpha=0.1)
|
|
quantized_weight.round_()
|
|
else:
|
|
if use_stochastic_rounding:
|
|
mantissa_difference = 1 << (23 - dtype_dict[weights_dtype]["mantissa"])
|
|
quantized_weight = quantized_weight.to(dtype=torch.float32).view(dtype=torch.int32)
|
|
quantized_weight = quantized_weight.add_(torch.randint_like(quantized_weight, low=0, high=mantissa_difference, dtype=torch.int32)).bitwise_and_(-mantissa_difference).view(dtype=torch.float32)
|
|
quantized_weight.nan_to_num_()
|
|
quantized_weight = quantized_weight.clamp_(dtype_dict[weights_dtype]["min"], dtype_dict[weights_dtype]["max"]).to(dtype_dict[weights_dtype]["torch_dtype"])
|
|
return quantized_weight, scale, zero_point
|
|
|
|
|
|
@devices.inference_context()
|
|
def apply_svdquant(weight: torch.FloatTensor, rank: int = 32, niter: int = 8, dtype: torch.dtype = None) -> tuple[torch.FloatTensor, torch.FloatTensor, torch.FloatTensor]:
|
|
reshape_weight = False
|
|
if weight.ndim > 2: # convs
|
|
reshape_weight = True
|
|
weight_shape = weight.shape
|
|
weight = weight.flatten(1,-1)
|
|
if weight.dtype != torch.float64:
|
|
weight = weight.to(dtype=torch.float32)
|
|
U, S, svd_down = torch.svd_lowrank(weight, q=rank, niter=niter)
|
|
svd_up = torch.mul(U, S.unsqueeze(0))
|
|
svd_down = svd_down.t_()
|
|
if dtype is not None:
|
|
svd_up = svd_up.to(dtype=dtype)
|
|
svd_down = svd_down.to(dtype=dtype)
|
|
weight = weight.sub(torch.mm(svd_up, svd_down))
|
|
if reshape_weight:
|
|
weight = weight.unflatten(-1, (*weight_shape[1:],)) # pylint: disable=possibly-used-before-assignment
|
|
return weight, svd_up, svd_down
|
|
|
|
|
|
@devices.inference_context()
|
|
def build_hadamard_n2(n: int, dtype: torch.dtype | None = None, device: torch.device | None = None) -> torch.FloatTensor:
|
|
current_size = 2
|
|
H = H_N2 = torch.tensor([[1, 1], [1, -1]], dtype=dtype, device=device)
|
|
while current_size < n:
|
|
H = torch.kron(H, H_N2)
|
|
current_size *= 2
|
|
H = H.div_(n**0.5)
|
|
H = prepare_weight_for_matmul(H)
|
|
return H
|
|
|
|
|
|
@devices.inference_context()
|
|
def build_hadamard_n4(n: int, dtype: torch.dtype | None = None, device: torch.device | None = None) -> torch.FloatTensor:
|
|
current_size = 4
|
|
H = H_N4 = torch.tensor([[ 1, 1, 1, -1], [ 1, 1, -1, 1], [ 1, -1, 1, 1], [-1, 1, 1, 1]], dtype=dtype, device=device)
|
|
while current_size < n:
|
|
H = torch.kron(H, H_N4)
|
|
current_size *= 4
|
|
H = H.div_(n**0.5)
|
|
H = prepare_weight_for_matmul(H)
|
|
return H
|
|
|
|
|
|
@devices.inference_context()
|
|
def build_hadamard(n: int, dtype: torch.dtype | None = None, device: torch.device | None = None) -> torch.FloatTensor:
|
|
if math.log(n, 4).is_integer():
|
|
return build_hadamard_n4(n, device=device, dtype=dtype)
|
|
elif math.log(n, 2).is_integer():
|
|
return build_hadamard_n2(n, device=device, dtype=dtype)
|
|
else:
|
|
raise RuntimeError("Hadamard Group Size must be a power of 2.")
|
|
|
|
|
|
# 256x256 Hadamard matrix is just 256 KB at FP32
|
|
# And is the exact same matrix on all model layers
|
|
# So we can safely cache a single one
|
|
HADAMARD_MATRIX_CACHE: dict[tuple[int, torch.device, torch.dtype], torch.FloatTensor] = {}
|
|
@devices.inference_context()
|
|
def get_hadamard(n: int, dtype: torch.dtype | None = None, device: torch.device | None = None) -> torch.FloatTensor:
|
|
device = devices.normalize_device(device)
|
|
H_key = (n, device, dtype)
|
|
H = HADAMARD_MATRIX_CACHE.get(H_key, None)
|
|
if H is None:
|
|
H = build_hadamard(n, dtype=dtype, device=device)
|
|
HADAMARD_MATRIX_CACHE[H_key] = H
|
|
return H
|
|
|
|
|
|
@devices.inference_context()
|
|
def rotate_hadamard(weight: torch.Tensor, group_size: int = 256, hadamard: torch.FloatTensor | None = None, is_conv: bool = False) -> torch.Tensor:
|
|
if hadamard is None:
|
|
hadamard = get_hadamard(group_size, dtype=weight.dtype, device=weight.device)
|
|
else:
|
|
group_size = hadamard.shape[-1]
|
|
if is_conv:
|
|
weight_shape = list(weight.shape)[1:]
|
|
weight = weight.flatten(1,-1)
|
|
weight = weight.unflatten(-1, (-1,group_size))
|
|
result = torch.matmul(weight, hadamard).flatten(-2,-1)
|
|
del hadamard
|
|
if is_conv:
|
|
result = result.unflatten(-1, weight_shape)
|
|
return result
|
|
|
|
|
|
@devices.inference_context()
|
|
def apply_hadamard(weight: torch.Tensor, group_size: int = 256, hadamard: torch.FloatTensor | None = None, layer_class_name: str | None = None) -> torch.Tensor:
|
|
is_conv = False
|
|
use_hadamard = True
|
|
if layer_class_name in conv_types or layer_class_name in conv_transpose_types:
|
|
is_conv = True
|
|
channel_size = weight.shape[1]
|
|
else:
|
|
channel_size = weight.shape[-1]
|
|
if channel_size % group_size != 0:
|
|
hadamard_pow2 = int(math.log2(group_size))
|
|
while channel_size % group_size != 0:
|
|
hadamard_pow2 -= 1
|
|
group_size = 2 ** hadamard_pow2
|
|
if group_size < 4:
|
|
use_hadamard = False
|
|
if use_hadamard:
|
|
weight = rotate_hadamard(weight, group_size=group_size, hadamard=hadamard, is_conv=is_conv)
|
|
return weight, use_hadamard, group_size
|
|
|
|
|
|
@devices.inference_context()
|
|
def prepare_weight_for_matmul(weight: torch.Tensor) -> torch.Tensor:
|
|
if use_contiguous_mm:
|
|
weight = weight.contiguous()
|
|
elif weight.is_contiguous():
|
|
weight = weight.t_().contiguous().t_()
|
|
return weight
|
|
|
|
|
|
@devices.inference_context()
|
|
def prepare_svd_for_matmul(svd_up: torch.FloatTensor, svd_down: torch.FloatTensor, use_quantized_matmul: bool) -> tuple[torch.FloatTensor, torch.FloatTensor]:
|
|
if svd_up is not None:
|
|
if use_quantized_matmul:
|
|
svd_up = prepare_weight_for_matmul(svd_up)
|
|
else:
|
|
svd_up = svd_up.contiguous()
|
|
if svd_down is not None:
|
|
svd_down = prepare_weight_for_matmul(svd_down)
|
|
return svd_up, svd_down
|
|
|
|
|
|
@devices.inference_context()
|
|
def quantize_int_mm(input: torch.FloatTensor, dim: int = -1, hadamard: torch.FloatTensor | None = None, matmul_dtype: str = "int8") -> tuple[torch.Tensor, torch.FloatTensor]:
|
|
if hadamard is not None:
|
|
input = rotate_hadamard(input, hadamard=hadamard)
|
|
scale = torch.amax(input.abs(), dim=dim, keepdims=True).div_(dtype_dict[matmul_dtype]["max"])
|
|
input = torch.div(input, scale).round_().clamp_(dtype_dict[matmul_dtype]["min"], dtype_dict[matmul_dtype]["max"]).to(dtype=dtype_dict[matmul_dtype]["torch_dtype"])
|
|
return input, scale
|
|
|
|
|
|
@devices.inference_context()
|
|
def quantize_int_mm_sr(input: torch.FloatTensor, dim: int = -1, hadamard: torch.FloatTensor | None = None, matmul_dtype: str = "int8") -> tuple[torch.Tensor, torch.FloatTensor]:
|
|
if hadamard is not None:
|
|
input = rotate_hadamard(input, hadamard=hadamard)
|
|
scale = torch.amax(input.abs(), dim=dim, keepdims=True).div_(dtype_dict[matmul_dtype]["max"])
|
|
input = torch.div(input, scale).add_(torch.randn_like(input), alpha=0.1).round_().clamp_(dtype_dict[matmul_dtype]["min"], dtype_dict[matmul_dtype]["max"]).to(dtype=dtype_dict[matmul_dtype]["torch_dtype"])
|
|
return input, scale
|
|
|
|
|
|
@devices.inference_context()
|
|
def quantize_fp_mm(input: torch.FloatTensor, dim: int = -1, hadamard: torch.FloatTensor | None = None, matmul_dtype: str = "float8_e4m3fn") -> tuple[torch.Tensor, torch.FloatTensor]:
|
|
if hadamard is not None:
|
|
input = rotate_hadamard(input, hadamard=hadamard)
|
|
scale = torch.amax(input.abs(), dim=dim, keepdims=True).div_(dtype_dict[matmul_dtype]["max"])
|
|
input = torch.div(input, scale).nan_to_num_().clamp_(dtype_dict[matmul_dtype]["min"], dtype_dict[matmul_dtype]["max"]).to(dtype=dtype_dict[matmul_dtype]["torch_dtype"])
|
|
return input, scale
|
|
|
|
|
|
@devices.inference_context()
|
|
def quantize_fp_mm_sr(input: torch.FloatTensor, dim: int = -1, hadamard: torch.FloatTensor | None = None, matmul_dtype: str = "float8_e4m3fn") -> tuple[torch.Tensor, torch.FloatTensor]:
|
|
if hadamard is not None:
|
|
input = rotate_hadamard(input, hadamard=hadamard)
|
|
mantissa_difference = 1 << (23 - dtype_dict[matmul_dtype]["mantissa"])
|
|
scale = torch.amax(input.abs(), dim=dim, keepdims=True).div_(dtype_dict[matmul_dtype]["max"])
|
|
input = torch.div(input, scale).to(dtype=torch.float32).view(dtype=torch.int32)
|
|
input = input.add_(torch.randint_like(input, low=0, high=mantissa_difference, dtype=torch.int32)).bitwise_and_(-mantissa_difference).view(dtype=torch.float32)
|
|
input = input.nan_to_num_().clamp_(dtype_dict[matmul_dtype]["min"], dtype_dict[matmul_dtype]["max"]).to(dtype=dtype_dict[matmul_dtype]["torch_dtype"])
|
|
return input, scale
|