Spaces:
Configuration error
Configuration error
import os | |
import sys | |
import traceback | |
import logging | |
logger = logging.getLogger(__name__) | |
from time import time as ttime | |
import fairseq | |
import faiss | |
import numpy as np | |
import parselmouth | |
import pyworld | |
import scipy.signal as signal | |
import torch | |
import torch.nn as nn | |
import torch.nn.functional as F | |
import torchcrepe | |
from infer.lib.infer_pack.models import ( | |
SynthesizerTrnMs256NSFsid, | |
SynthesizerTrnMs256NSFsid_nono, | |
SynthesizerTrnMs768NSFsid, | |
SynthesizerTrnMs768NSFsid_nono, | |
) | |
now_dir = os.getcwd() | |
sys.path.append(now_dir) | |
from multiprocessing import Manager as M | |
from configs.config import Config | |
config = Config() | |
mm = M() | |
if config.dml == True: | |
def forward_dml(ctx, x, scale): | |
ctx.scale = scale | |
res = x.clone().detach() | |
return res | |
fairseq.modules.grad_multiply.GradMultiply.forward = forward_dml | |
# config.device=torch.device("cpu")########强制cpu测试 | |
# config.is_half=False########强制cpu测试 | |
class RVC: | |
def __init__( | |
self, | |
key, | |
pth_path, | |
index_path, | |
index_rate, | |
n_cpu, | |
inp_q, | |
opt_q, | |
device, | |
last_rvc=None, | |
) -> None: | |
""" | |
初始化 | |
""" | |
try: | |
global config | |
self.inp_q = inp_q | |
self.opt_q = opt_q | |
# device="cpu"########强制cpu测试 | |
self.device = device | |
self.f0_up_key = key | |
self.time_step = 160 / 16000 * 1000 | |
self.f0_min = 50 | |
self.f0_max = 1100 | |
self.f0_mel_min = 1127 * np.log(1 + self.f0_min / 700) | |
self.f0_mel_max = 1127 * np.log(1 + self.f0_max / 700) | |
self.sr = 16000 | |
self.window = 160 | |
self.n_cpu = n_cpu | |
if index_rate != 0: | |
self.index = faiss.read_index(index_path) | |
self.big_npy = self.index.reconstruct_n(0, self.index.ntotal) | |
logger.info("Index search enabled") | |
self.pth_path = pth_path | |
self.index_path = index_path | |
self.index_rate = index_rate | |
if last_rvc is None: | |
models, _, _ = fairseq.checkpoint_utils.load_model_ensemble_and_task( | |
["assets/hubert/hubert_base.pt"], | |
suffix="", | |
) | |
hubert_model = models[0] | |
hubert_model = hubert_model.to(device) | |
if config.is_half: | |
hubert_model = hubert_model.half() | |
else: | |
hubert_model = hubert_model.float() | |
hubert_model.eval() | |
self.model = hubert_model | |
else: | |
self.model = last_rvc.model | |
if last_rvc is None or last_rvc.pth_path != self.pth_path: | |
cpt = torch.load(self.pth_path, map_location="cpu") | |
self.tgt_sr = cpt["config"][-1] | |
cpt["config"][-3] = cpt["weight"]["emb_g.weight"].shape[0] | |
self.if_f0 = cpt.get("f0", 1) | |
self.version = cpt.get("version", "v1") | |
if self.version == "v1": | |
if self.if_f0 == 1: | |
self.net_g = SynthesizerTrnMs256NSFsid( | |
*cpt["config"], is_half=config.is_half | |
) | |
else: | |
self.net_g = SynthesizerTrnMs256NSFsid_nono(*cpt["config"]) | |
elif self.version == "v2": | |
if self.if_f0 == 1: | |
self.net_g = SynthesizerTrnMs768NSFsid( | |
*cpt["config"], is_half=config.is_half | |
) | |
else: | |
self.net_g = SynthesizerTrnMs768NSFsid_nono(*cpt["config"]) | |
del self.net_g.enc_q | |
logger.debug(self.net_g.load_state_dict(cpt["weight"], strict=False)) | |
self.net_g.eval().to(device) | |
# print(2333333333,device,config.device,self.device)#net_g是device,hubert是config.device | |
if config.is_half: | |
self.net_g = self.net_g.half() | |
else: | |
self.net_g = self.net_g.float() | |
self.is_half = config.is_half | |
else: | |
self.tgt_sr = last_rvc.tgt_sr | |
self.if_f0 = last_rvc.if_f0 | |
self.version = last_rvc.version | |
self.net_g = last_rvc.net_g | |
self.is_half = last_rvc.is_half | |
if last_rvc is not None and hasattr(last_rvc, "model_rmvpe"): | |
self.model_rmvpe = last_rvc.model_rmvpe | |
except: | |
logger.warn(traceback.format_exc()) | |
def change_key(self, new_key): | |
self.f0_up_key = new_key | |
def change_index_rate(self, new_index_rate): | |
if new_index_rate != 0 and self.index_rate == 0: | |
self.index = faiss.read_index(self.index_path) | |
self.big_npy = self.index.reconstruct_n(0, self.index.ntotal) | |
logger.info("Index search enabled") | |
self.index_rate = new_index_rate | |
def get_f0_post(self, f0): | |
f0_min = self.f0_min | |
f0_max = self.f0_max | |
f0_mel_min = 1127 * np.log(1 + f0_min / 700) | |
f0_mel_max = 1127 * np.log(1 + f0_max / 700) | |
f0bak = f0.copy() | |
f0_mel = 1127 * np.log(1 + f0 / 700) | |
f0_mel[f0_mel > 0] = (f0_mel[f0_mel > 0] - f0_mel_min) * 254 / ( | |
f0_mel_max - f0_mel_min | |
) + 1 | |
f0_mel[f0_mel <= 1] = 1 | |
f0_mel[f0_mel > 255] = 255 | |
f0_coarse = np.rint(f0_mel).astype(np.int32) | |
return f0_coarse, f0bak | |
def get_f0(self, x, f0_up_key, n_cpu, method="harvest"): | |
n_cpu = int(n_cpu) | |
if method == "crepe": | |
return self.get_f0_crepe(x, f0_up_key) | |
if method == "rmvpe": | |
return self.get_f0_rmvpe(x, f0_up_key) | |
if method == "pm": | |
p_len = x.shape[0] // 160 + 1 | |
f0 = ( | |
parselmouth.Sound(x, 16000) | |
.to_pitch_ac( | |
time_step=0.01, | |
voicing_threshold=0.6, | |
pitch_floor=50, | |
pitch_ceiling=1100, | |
) | |
.selected_array["frequency"] | |
) | |
pad_size = (p_len - len(f0) + 1) // 2 | |
if pad_size > 0 or p_len - len(f0) - pad_size > 0: | |
# print(pad_size, p_len - len(f0) - pad_size) | |
f0 = np.pad( | |
f0, [[pad_size, p_len - len(f0) - pad_size]], mode="constant" | |
) | |
f0 *= pow(2, f0_up_key / 12) | |
return self.get_f0_post(f0) | |
if n_cpu == 1: | |
f0, t = pyworld.harvest( | |
x.astype(np.double), | |
fs=16000, | |
f0_ceil=1100, | |
f0_floor=50, | |
frame_period=10, | |
) | |
f0 = signal.medfilt(f0, 3) | |
f0 *= pow(2, f0_up_key / 12) | |
return self.get_f0_post(f0) | |
f0bak = np.zeros(x.shape[0] // 160 + 1, dtype=np.float64) | |
length = len(x) | |
part_length = 160 * ((length // 160 - 1) // n_cpu + 1) | |
n_cpu = (length // 160 - 1) // (part_length // 160) + 1 | |
ts = ttime() | |
res_f0 = mm.dict() | |
for idx in range(n_cpu): | |
tail = part_length * (idx + 1) + 320 | |
if idx == 0: | |
self.inp_q.put((idx, x[:tail], res_f0, n_cpu, ts)) | |
else: | |
self.inp_q.put( | |
(idx, x[part_length * idx - 320 : tail], res_f0, n_cpu, ts) | |
) | |
while 1: | |
res_ts = self.opt_q.get() | |
if res_ts == ts: | |
break | |
f0s = [i[1] for i in sorted(res_f0.items(), key=lambda x: x[0])] | |
for idx, f0 in enumerate(f0s): | |
if idx == 0: | |
f0 = f0[:-3] | |
elif idx != n_cpu - 1: | |
f0 = f0[2:-3] | |
else: | |
f0 = f0[2:] | |
f0bak[ | |
part_length * idx // 160 : part_length * idx // 160 + f0.shape[0] | |
] = f0 | |
f0bak = signal.medfilt(f0bak, 3) | |
f0bak *= pow(2, f0_up_key / 12) | |
return self.get_f0_post(f0bak) | |
def get_f0_crepe(self, x, f0_up_key): | |
if "privateuseone" in str(self.device): ###不支持dml,cpu又太慢用不成,拿pm顶替 | |
return self.get_f0(x, f0_up_key, 1, "pm") | |
audio = torch.tensor(np.copy(x))[None].float() | |
# print("using crepe,device:%s"%self.device) | |
f0, pd = torchcrepe.predict( | |
audio, | |
self.sr, | |
160, | |
self.f0_min, | |
self.f0_max, | |
"full", | |
batch_size=512, | |
# device=self.device if self.device.type!="privateuseone" else "cpu",###crepe不用半精度全部是全精度所以不愁###cpu延迟高到没法用 | |
device=self.device, | |
return_periodicity=True, | |
) | |
pd = torchcrepe.filter.median(pd, 3) | |
f0 = torchcrepe.filter.mean(f0, 3) | |
f0[pd < 0.1] = 0 | |
f0 = f0[0].cpu().numpy() | |
f0 *= pow(2, f0_up_key / 12) | |
return self.get_f0_post(f0) | |
def get_f0_rmvpe(self, x, f0_up_key): | |
if hasattr(self, "model_rmvpe") == False: | |
from infer.lib.rmvpe import RMVPE | |
logger.info("Loading rmvpe model") | |
self.model_rmvpe = RMVPE( | |
# "rmvpe.pt", is_half=self.is_half if self.device.type!="privateuseone" else False, device=self.device if self.device.type!="privateuseone"else "cpu"####dml时强制对rmvpe用cpu跑 | |
# "rmvpe.pt", is_half=False, device=self.device####dml配置 | |
# "rmvpe.pt", is_half=False, device="cpu"####锁定cpu配置 | |
"assets/rmvpe/rmvpe.pt", | |
is_half=self.is_half, | |
device=self.device, ####正常逻辑 | |
) | |
# self.model_rmvpe = RMVPE("aug2_58000_half.pt", is_half=self.is_half, device=self.device) | |
f0 = self.model_rmvpe.infer_from_audio(x, thred=0.03) | |
f0 *= pow(2, f0_up_key / 12) | |
return self.get_f0_post(f0) | |
def infer( | |
self, | |
feats: torch.Tensor, | |
indata: np.ndarray, | |
block_frame_16k, | |
rate, | |
cache_pitch, | |
cache_pitchf, | |
f0method, | |
) -> np.ndarray: | |
feats = feats.view(1, -1) | |
if config.is_half: | |
feats = feats.half() | |
else: | |
feats = feats.float() | |
feats = feats.to(self.device) | |
t1 = ttime() | |
with torch.no_grad(): | |
padding_mask = torch.BoolTensor(feats.shape).to(self.device).fill_(False) | |
inputs = { | |
"source": feats, | |
"padding_mask": padding_mask, | |
"output_layer": 9 if self.version == "v1" else 12, | |
} | |
logits = self.model.extract_features(**inputs) | |
feats = ( | |
self.model.final_proj(logits[0]) if self.version == "v1" else logits[0] | |
) | |
feats = F.pad(feats, (0, 0, 1, 0)) | |
t2 = ttime() | |
try: | |
if hasattr(self, "index") and self.index_rate != 0: | |
leng_replace_head = int(rate * feats[0].shape[0]) | |
npy = feats[0][-leng_replace_head:].cpu().numpy().astype("float32") | |
score, ix = self.index.search(npy, k=8) | |
weight = np.square(1 / score) | |
weight /= weight.sum(axis=1, keepdims=True) | |
npy = np.sum(self.big_npy[ix] * np.expand_dims(weight, axis=2), axis=1) | |
if config.is_half: | |
npy = npy.astype("float16") | |
feats[0][-leng_replace_head:] = ( | |
torch.from_numpy(npy).unsqueeze(0).to(self.device) * self.index_rate | |
+ (1 - self.index_rate) * feats[0][-leng_replace_head:] | |
) | |
else: | |
logger.warn("Index search FAILED or disabled") | |
except: | |
traceback.print_exc() | |
logger.warn("Index search FAILED") | |
feats = F.interpolate(feats.permute(0, 2, 1), scale_factor=2).permute(0, 2, 1) | |
t3 = ttime() | |
if self.if_f0 == 1: | |
pitch, pitchf = self.get_f0(indata, self.f0_up_key, self.n_cpu, f0method) | |
start_frame = block_frame_16k // 160 | |
end_frame = len(cache_pitch) - (pitch.shape[0] - 4) + start_frame | |
cache_pitch[:] = np.append(cache_pitch[start_frame:end_frame], pitch[3:-1]) | |
cache_pitchf[:] = np.append( | |
cache_pitchf[start_frame:end_frame], pitchf[3:-1] | |
) | |
p_len = min(feats.shape[1], 13000, cache_pitch.shape[0]) | |
else: | |
cache_pitch, cache_pitchf = None, None | |
p_len = min(feats.shape[1], 13000) | |
t4 = ttime() | |
feats = feats[:, :p_len, :] | |
if self.if_f0 == 1: | |
cache_pitch = cache_pitch[:p_len] | |
cache_pitchf = cache_pitchf[:p_len] | |
cache_pitch = torch.LongTensor(cache_pitch).unsqueeze(0).to(self.device) | |
cache_pitchf = torch.FloatTensor(cache_pitchf).unsqueeze(0).to(self.device) | |
p_len = torch.LongTensor([p_len]).to(self.device) | |
ii = 0 # sid | |
sid = torch.LongTensor([ii]).to(self.device) | |
with torch.no_grad(): | |
if self.if_f0 == 1: | |
# print(12222222222,feats.device,p_len.device,cache_pitch.device,cache_pitchf.device,sid.device,rate2) | |
infered_audio = ( | |
self.net_g.infer( | |
feats, p_len, cache_pitch, cache_pitchf, sid, rate | |
)[0][0, 0] | |
.data | |
.float() | |
) | |
else: | |
infered_audio = ( | |
self.net_g.infer(feats, p_len, sid, rate)[0][0, 0] | |
.data | |
.float() | |
) | |
t5 = ttime() | |
logger.info( | |
"Spent time: fea = %.2fs, index = %.2fs, f0 = %.2fs, model = %.2fs", | |
t2 - t1, | |
t3 - t2, | |
t4 - t3, | |
t5 - t4, | |
) | |
return infered_audio |