ORI-Muchim commited on
Commit
295a0ef
1 Parent(s): b0c1e24

Upload 16 files

Browse files
Files changed (16) hide show
  1. __init__.py +0 -0
  2. app.py +140 -0
  3. attentions.py +303 -0
  4. commons.py +161 -0
  5. data_utils.py +392 -0
  6. inference.py +37 -0
  7. losses.py +71 -0
  8. mel_processing.py +112 -0
  9. models.py +731 -0
  10. modules.py +390 -0
  11. pqmf.py +116 -0
  12. requirements.txt +14 -0
  13. stft.py +209 -0
  14. stft_loss.py +136 -0
  15. transforms.py +193 -0
  16. utils.py +258 -0
__init__.py ADDED
File without changes
app.py ADDED
@@ -0,0 +1,140 @@
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
1
+ import json
2
+ import os
3
+ import re
4
+ import librosa
5
+ import numpy as np
6
+ import torch
7
+ from torch import no_grad, LongTensor
8
+ import commons
9
+ import utils
10
+ import gradio as gr
11
+ from models import SynthesizerTrn
12
+ from text import text_to_sequence
13
+ from text.symbols import symbols
14
+
15
+ limitation = os.getenv("SYSTEM") == "spaces" # limit text and audio length in huggingface spaces
16
+
17
+
18
+ def get_text(text, hps):
19
+ text_norm = text_to_sequence(text, hps.data.text_cleaners)
20
+ if hps.data.add_blank:
21
+ text_norm = commons.intersperse(text_norm, 0)
22
+ text_norm = torch.LongTensor(text_norm)
23
+ return text_norm
24
+
25
+
26
+ def create_tts_fn(net_g, hps, speaker_ids):
27
+ def tts_fn(text, speaker, speed):
28
+ if limitation:
29
+ text_len = len(text)
30
+ max_len = 700
31
+ if text_len > max_len:
32
+ return "Error: Text is too long", None
33
+
34
+ speaker_id = speaker_ids[speaker]
35
+ stn_tst = get_text(text, hps)
36
+
37
+ with no_grad():
38
+ x_tst = stn_tst.unsqueeze(0)
39
+ x_tst_lengths = LongTensor([stn_tst.size(0)])
40
+ sid = LongTensor([speaker_id])
41
+ audio = net_g.infer(x_tst, x_tst_lengths, sid=sid, noise_scale=.667, noise_scale_w=0.8,
42
+ length_scale=1.0 / speed)[0][0, 0].data.cpu().float().numpy()
43
+ del stn_tst, x_tst, x_tst_lengths, sid
44
+ return "Success", (hps.data.sampling_rate, audio)
45
+
46
+ return tts_fn
47
+
48
+
49
+ css = """
50
+ #advanced-btn {
51
+ color: white;
52
+ border-color: black;
53
+ background: black;
54
+ font-size: .7rem !important;
55
+ line-height: 19px;
56
+ margin-top: 24px;
57
+ margin-bottom: 12px;
58
+ padding: 2px 8px;
59
+ border-radius: 14px !important;
60
+ }
61
+ #advanced-options {
62
+ display: none;
63
+ margin-bottom: 20px;
64
+ }
65
+ """
66
+
67
+ if __name__ == '__main__':
68
+ models_tts = []
69
+ name = 'AronaTTS'
70
+ lang = '日本語 (Japanese)'
71
+ example = '내가 누군가를 좋아한다는 사실이 그 사람에게는 상처가 될 수 있잖아요.'
72
+ config_path = f"saved_model/config.json"
73
+ model_path = f"saved_model/model.pth"
74
+ cover_path = f"saved_model/cover.png"
75
+
76
+ hps = utils.get_hparams_from_file(config_path)
77
+
78
+ net_g = SynthesizerTrn(
79
+ len(symbols),
80
+ hps.data.filter_length // 2 + 1,
81
+ hps.train.segment_size // hps.data.hop_length,
82
+ n_speakers=hps.data.n_speakers,
83
+ **hps.model).cuda()
84
+ _ = net_g.eval()
85
+
86
+ utils.load_checkpoint(model_path, net_g, None)
87
+
88
+ net_g.eval()
89
+
90
+ speaker_ids = [0]
91
+ speakers = [name]
92
+
93
+ t = 'vits'
94
+ models_tts.append((name, cover_path, speakers, lang, example,
95
+ hps.symbols, create_tts_fn(net_g, hps, speaker_ids)))
96
+
97
+ app = gr.Blocks(css=css)
98
+
99
+ with app:
100
+ gr.Markdown("# BlueArchive AronaTTS Using VITS Model\n"
101
+ "![visitor badge](https://visitor-badge.glitch.me/badge?page_id=openduckparty.AronaTTS)\n\n")
102
+
103
+ for i, (name, cover_path, speakers, lang, example, symbols, tts_fn
104
+ ) in enumerate(models_tts):
105
+
106
+ with gr.Column():
107
+ gr.Markdown(f"## {name}\n\n"
108
+ f"![cover](file/{cover_path})\n\n"
109
+ f"lang: {lang}")
110
+ tts_input1 = gr.TextArea(label="Text (700 words limitation)", value=example,
111
+ elem_id=f"tts-input{i}")
112
+ tts_input2 = gr.Dropdown(label="Speaker", choices=speakers,
113
+ type="index", value=speakers[0])
114
+ tts_input3 = gr.Slider(label="Speed", value=1, minimum=0.1, maximum=2, step=0.1)
115
+ tts_submit = gr.Button("Generate", variant="primary")
116
+ tts_output1 = gr.Textbox(label="Output Message")
117
+ tts_output2 = gr.Audio(label="Output Audio")
118
+ tts_submit.click(tts_fn, [tts_input1, tts_input2, tts_input3],
119
+ [tts_output1, tts_output2])
120
+ _js=f"""
121
+ (i,phonemes) => {{
122
+ let root = document.querySelector("body > gradio-app");
123
+ if (root.shadowRoot != null)
124
+ root = root.shadowRoot;
125
+ let text_input = root.querySelector("#tts-input{i}").querySelector("textarea");
126
+ let startPos = text_input.selectionStart;
127
+ let endPos = text_input.selectionEnd;
128
+ let oldTxt = text_input.value;
129
+ let result = oldTxt.substring(0, startPos) + phonemes[i] + oldTxt.substring(endPos);
130
+ text_input.value = result;
131
+ let x = window.scrollX, y = window.scrollY;
132
+ text_input.focus();
133
+ text_input.selectionStart = startPos + phonemes[i].length;
134
+ text_input.selectionEnd = startPos + phonemes[i].length;
135
+ text_input.blur();
136
+ window.scrollTo(x, y);
137
+ return [];
138
+ }}"""
139
+
140
+ app.queue(concurrency_count=3).launch(show_api=False)
attentions.py ADDED
@@ -0,0 +1,303 @@
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
1
+ import copy
2
+ import math
3
+ import numpy as np
4
+ import torch
5
+ from torch import nn
6
+ from torch.nn import functional as F
7
+
8
+ import commons
9
+ import modules
10
+ from modules import LayerNorm
11
+
12
+
13
+ class Encoder(nn.Module):
14
+ def __init__(self, hidden_channels, filter_channels, n_heads, n_layers, kernel_size=1, p_dropout=0., window_size=4, **kwargs):
15
+ super().__init__()
16
+ self.hidden_channels = hidden_channels
17
+ self.filter_channels = filter_channels
18
+ self.n_heads = n_heads
19
+ self.n_layers = n_layers
20
+ self.kernel_size = kernel_size
21
+ self.p_dropout = p_dropout
22
+ self.window_size = window_size
23
+
24
+ self.drop = nn.Dropout(p_dropout)
25
+ self.attn_layers = nn.ModuleList()
26
+ self.norm_layers_1 = nn.ModuleList()
27
+ self.ffn_layers = nn.ModuleList()
28
+ self.norm_layers_2 = nn.ModuleList()
29
+ for i in range(self.n_layers):
30
+ self.attn_layers.append(MultiHeadAttention(hidden_channels, hidden_channels, n_heads, p_dropout=p_dropout, window_size=window_size))
31
+ self.norm_layers_1.append(LayerNorm(hidden_channels))
32
+ self.ffn_layers.append(FFN(hidden_channels, hidden_channels, filter_channels, kernel_size, p_dropout=p_dropout))
33
+ self.norm_layers_2.append(LayerNorm(hidden_channels))
34
+
35
+ def forward(self, x, x_mask):
36
+ attn_mask = x_mask.unsqueeze(2) * x_mask.unsqueeze(-1)
37
+ x = x * x_mask
38
+ for i in range(self.n_layers):
39
+ y = self.attn_layers[i](x, x, attn_mask)
40
+ y = self.drop(y)
41
+ x = self.norm_layers_1[i](x + y)
42
+
43
+ y = self.ffn_layers[i](x, x_mask)
44
+ y = self.drop(y)
45
+ x = self.norm_layers_2[i](x + y)
46
+ x = x * x_mask
47
+ return x
48
+
49
+
50
+ class Decoder(nn.Module):
51
+ def __init__(self, hidden_channels, filter_channels, n_heads, n_layers, kernel_size=1, p_dropout=0., proximal_bias=False, proximal_init=True, **kwargs):
52
+ super().__init__()
53
+ self.hidden_channels = hidden_channels
54
+ self.filter_channels = filter_channels
55
+ self.n_heads = n_heads
56
+ self.n_layers = n_layers
57
+ self.kernel_size = kernel_size
58
+ self.p_dropout = p_dropout
59
+ self.proximal_bias = proximal_bias
60
+ self.proximal_init = proximal_init
61
+
62
+ self.drop = nn.Dropout(p_dropout)
63
+ self.self_attn_layers = nn.ModuleList()
64
+ self.norm_layers_0 = nn.ModuleList()
65
+ self.encdec_attn_layers = nn.ModuleList()
66
+ self.norm_layers_1 = nn.ModuleList()
67
+ self.ffn_layers = nn.ModuleList()
68
+ self.norm_layers_2 = nn.ModuleList()
69
+ for i in range(self.n_layers):
70
+ self.self_attn_layers.append(MultiHeadAttention(hidden_channels, hidden_channels, n_heads, p_dropout=p_dropout, proximal_bias=proximal_bias, proximal_init=proximal_init))
71
+ self.norm_layers_0.append(LayerNorm(hidden_channels))
72
+ self.encdec_attn_layers.append(MultiHeadAttention(hidden_channels, hidden_channels, n_heads, p_dropout=p_dropout))
73
+ self.norm_layers_1.append(LayerNorm(hidden_channels))
74
+ self.ffn_layers.append(FFN(hidden_channels, hidden_channels, filter_channels, kernel_size, p_dropout=p_dropout, causal=True))
75
+ self.norm_layers_2.append(LayerNorm(hidden_channels))
76
+
77
+ def forward(self, x, x_mask, h, h_mask):
78
+ """
79
+ x: decoder input
80
+ h: encoder output
81
+ """
82
+ self_attn_mask = commons.subsequent_mask(x_mask.size(2)).to(device=x.device, dtype=x.dtype)
83
+ encdec_attn_mask = h_mask.unsqueeze(2) * x_mask.unsqueeze(-1)
84
+ x = x * x_mask
85
+ for i in range(self.n_layers):
86
+ y = self.self_attn_layers[i](x, x, self_attn_mask)
87
+ y = self.drop(y)
88
+ x = self.norm_layers_0[i](x + y)
89
+
90
+ y = self.encdec_attn_layers[i](x, h, encdec_attn_mask)
91
+ y = self.drop(y)
92
+ x = self.norm_layers_1[i](x + y)
93
+
94
+ y = self.ffn_layers[i](x, x_mask)
95
+ y = self.drop(y)
96
+ x = self.norm_layers_2[i](x + y)
97
+ x = x * x_mask
98
+ return x
99
+
100
+
101
+ class MultiHeadAttention(nn.Module):
102
+ def __init__(self, channels, out_channels, n_heads, p_dropout=0., window_size=None, heads_share=True, block_length=None, proximal_bias=False, proximal_init=False):
103
+ super().__init__()
104
+ assert channels % n_heads == 0
105
+
106
+ self.channels = channels
107
+ self.out_channels = out_channels
108
+ self.n_heads = n_heads
109
+ self.p_dropout = p_dropout
110
+ self.window_size = window_size
111
+ self.heads_share = heads_share
112
+ self.block_length = block_length
113
+ self.proximal_bias = proximal_bias
114
+ self.proximal_init = proximal_init
115
+ self.attn = None
116
+
117
+ self.k_channels = channels // n_heads
118
+ self.conv_q = nn.Conv1d(channels, channels, 1)
119
+ self.conv_k = nn.Conv1d(channels, channels, 1)
120
+ self.conv_v = nn.Conv1d(channels, channels, 1)
121
+ self.conv_o = nn.Conv1d(channels, out_channels, 1)
122
+ self.drop = nn.Dropout(p_dropout)
123
+
124
+ if window_size is not None:
125
+ n_heads_rel = 1 if heads_share else n_heads
126
+ rel_stddev = self.k_channels**-0.5
127
+ self.emb_rel_k = nn.Parameter(torch.randn(n_heads_rel, window_size * 2 + 1, self.k_channels) * rel_stddev)
128
+ self.emb_rel_v = nn.Parameter(torch.randn(n_heads_rel, window_size * 2 + 1, self.k_channels) * rel_stddev)
129
+
130
+ nn.init.xavier_uniform_(self.conv_q.weight)
131
+ nn.init.xavier_uniform_(self.conv_k.weight)
132
+ nn.init.xavier_uniform_(self.conv_v.weight)
133
+ if proximal_init:
134
+ with torch.no_grad():
135
+ self.conv_k.weight.copy_(self.conv_q.weight)
136
+ self.conv_k.bias.copy_(self.conv_q.bias)
137
+
138
+ def forward(self, x, c, attn_mask=None):
139
+ q = self.conv_q(x)
140
+ k = self.conv_k(c)
141
+ v = self.conv_v(c)
142
+
143
+ x, self.attn = self.attention(q, k, v, mask=attn_mask)
144
+
145
+ x = self.conv_o(x)
146
+ return x
147
+
148
+ def attention(self, query, key, value, mask=None):
149
+ # reshape [b, d, t] -> [b, n_h, t, d_k]
150
+ b, d, t_s, t_t = (*key.size(), query.size(2))
151
+ query = query.view(b, self.n_heads, self.k_channels, t_t).transpose(2, 3)
152
+ key = key.view(b, self.n_heads, self.k_channels, t_s).transpose(2, 3)
153
+ value = value.view(b, self.n_heads, self.k_channels, t_s).transpose(2, 3)
154
+
155
+ scores = torch.matmul(query / math.sqrt(self.k_channels), key.transpose(-2, -1))
156
+ if self.window_size is not None:
157
+ assert t_s == t_t, "Relative attention is only available for self-attention."
158
+ key_relative_embeddings = self._get_relative_embeddings(self.emb_rel_k, t_s)
159
+ rel_logits = self._matmul_with_relative_keys(query /math.sqrt(self.k_channels), key_relative_embeddings)
160
+ scores_local = self._relative_position_to_absolute_position(rel_logits)
161
+ scores = scores + scores_local
162
+ if self.proximal_bias:
163
+ assert t_s == t_t, "Proximal bias is only available for self-attention."
164
+ scores = scores + self._attention_bias_proximal(t_s).to(device=scores.device, dtype=scores.dtype)
165
+ if mask is not None:
166
+ scores = scores.masked_fill(mask == 0, -1e4)
167
+ if self.block_length is not None:
168
+ assert t_s == t_t, "Local attention is only available for self-attention."
169
+ block_mask = torch.ones_like(scores).triu(-self.block_length).tril(self.block_length)
170
+ scores = scores.masked_fill(block_mask == 0, -1e4)
171
+ p_attn = F.softmax(scores, dim=-1) # [b, n_h, t_t, t_s]
172
+ p_attn = self.drop(p_attn)
173
+ output = torch.matmul(p_attn, value)
174
+ if self.window_size is not None:
175
+ relative_weights = self._absolute_position_to_relative_position(p_attn)
176
+ value_relative_embeddings = self._get_relative_embeddings(self.emb_rel_v, t_s)
177
+ output = output + self._matmul_with_relative_values(relative_weights, value_relative_embeddings)
178
+ output = output.transpose(2, 3).contiguous().view(b, d, t_t) # [b, n_h, t_t, d_k] -> [b, d, t_t]
179
+ return output, p_attn
180
+
181
+ def _matmul_with_relative_values(self, x, y):
182
+ """
183
+ x: [b, h, l, m]
184
+ y: [h or 1, m, d]
185
+ ret: [b, h, l, d]
186
+ """
187
+ ret = torch.matmul(x, y.unsqueeze(0))
188
+ return ret
189
+
190
+ def _matmul_with_relative_keys(self, x, y):
191
+ """
192
+ x: [b, h, l, d]
193
+ y: [h or 1, m, d]
194
+ ret: [b, h, l, m]
195
+ """
196
+ ret = torch.matmul(x, y.unsqueeze(0).transpose(-2, -1))
197
+ return ret
198
+
199
+ def _get_relative_embeddings(self, relative_embeddings, length):
200
+ max_relative_position = 2 * self.window_size + 1
201
+ # Pad first before slice to avoid using cond ops.
202
+ pad_length = max(length - (self.window_size + 1), 0)
203
+ slice_start_position = max((self.window_size + 1) - length, 0)
204
+ slice_end_position = slice_start_position + 2 * length - 1
205
+ if pad_length > 0:
206
+ padded_relative_embeddings = F.pad(
207
+ relative_embeddings,
208
+ commons.convert_pad_shape([[0, 0], [pad_length, pad_length], [0, 0]]))
209
+ else:
210
+ padded_relative_embeddings = relative_embeddings
211
+ used_relative_embeddings = padded_relative_embeddings[:,slice_start_position:slice_end_position]
212
+ return used_relative_embeddings
213
+
214
+ def _relative_position_to_absolute_position(self, x):
215
+ """
216
+ x: [b, h, l, 2*l-1]
217
+ ret: [b, h, l, l]
218
+ """
219
+ batch, heads, length, _ = x.size()
220
+ # Concat columns of pad to shift from relative to absolute indexing.
221
+ x = F.pad(x, commons.convert_pad_shape([[0,0],[0,0],[0,0],[0,1]]))
222
+
223
+ # Concat extra elements so to add up to shape (len+1, 2*len-1).
224
+ x_flat = x.view([batch, heads, length * 2 * length])
225
+ x_flat = F.pad(x_flat, commons.convert_pad_shape([[0,0],[0,0],[0,length-1]]))
226
+
227
+ # Reshape and slice out the padded elements.
228
+ x_final = x_flat.view([batch, heads, length+1, 2*length-1])[:, :, :length, length-1:]
229
+ return x_final
230
+
231
+ def _absolute_position_to_relative_position(self, x):
232
+ """
233
+ x: [b, h, l, l]
234
+ ret: [b, h, l, 2*l-1]
235
+ """
236
+ batch, heads, length, _ = x.size()
237
+ # padd along column
238
+ x = F.pad(x, commons.convert_pad_shape([[0, 0], [0, 0], [0, 0], [0, length-1]]))
239
+ x_flat = x.view([batch, heads, length**2 + length*(length -1)])
240
+ # add 0's in the beginning that will skew the elements after reshape
241
+ x_flat = F.pad(x_flat, commons.convert_pad_shape([[0, 0], [0, 0], [length, 0]]))
242
+ x_final = x_flat.view([batch, heads, length, 2*length])[:,:,:,1:]
243
+ return x_final
244
+
245
+ def _attention_bias_proximal(self, length):
246
+ """Bias for self-attention to encourage attention to close positions.
247
+ Args:
248
+ length: an integer scalar.
249
+ Returns:
250
+ a Tensor with shape [1, 1, length, length]
251
+ """
252
+ r = torch.arange(length, dtype=torch.float32)
253
+ diff = torch.unsqueeze(r, 0) - torch.unsqueeze(r, 1)
254
+ return torch.unsqueeze(torch.unsqueeze(-torch.log1p(torch.abs(diff)), 0), 0)
255
+
256
+
257
+ class FFN(nn.Module):
258
+ def __init__(self, in_channels, out_channels, filter_channels, kernel_size, p_dropout=0., activation=None, causal=False):
259
+ super().__init__()
260
+ self.in_channels = in_channels
261
+ self.out_channels = out_channels
262
+ self.filter_channels = filter_channels
263
+ self.kernel_size = kernel_size
264
+ self.p_dropout = p_dropout
265
+ self.activation = activation
266
+ self.causal = causal
267
+
268
+ if causal:
269
+ self.padding = self._causal_padding
270
+ else:
271
+ self.padding = self._same_padding
272
+
273
+ self.conv_1 = nn.Conv1d(in_channels, filter_channels, kernel_size)
274
+ self.conv_2 = nn.Conv1d(filter_channels, out_channels, kernel_size)
275
+ self.drop = nn.Dropout(p_dropout)
276
+
277
+ def forward(self, x, x_mask):
278
+ x = self.conv_1(self.padding(x * x_mask))
279
+ if self.activation == "gelu":
280
+ x = x * torch.sigmoid(1.702 * x)
281
+ else:
282
+ x = torch.relu(x)
283
+ x = self.drop(x)
284
+ x = self.conv_2(self.padding(x * x_mask))
285
+ return x * x_mask
286
+
287
+ def _causal_padding(self, x):
288
+ if self.kernel_size == 1:
289
+ return x
290
+ pad_l = self.kernel_size - 1
291
+ pad_r = 0
292
+ padding = [[0, 0], [0, 0], [pad_l, pad_r]]
293
+ x = F.pad(x, commons.convert_pad_shape(padding))
294
+ return x
295
+
296
+ def _same_padding(self, x):
297
+ if self.kernel_size == 1:
298
+ return x
299
+ pad_l = (self.kernel_size - 1) // 2
300
+ pad_r = self.kernel_size // 2
301
+ padding = [[0, 0], [0, 0], [pad_l, pad_r]]
302
+ x = F.pad(x, commons.convert_pad_shape(padding))
303
+ return x
commons.py ADDED
@@ -0,0 +1,161 @@
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
1
+ import math
2
+ import numpy as np
3
+ import torch
4
+ from torch import nn
5
+ from torch.nn import functional as F
6
+
7
+
8
+ def init_weights(m, mean=0.0, std=0.01):
9
+ classname = m.__class__.__name__
10
+ if classname.find("Conv") != -1:
11
+ m.weight.data.normal_(mean, std)
12
+
13
+
14
+ def get_padding(kernel_size, dilation=1):
15
+ return int((kernel_size*dilation - dilation)/2)
16
+
17
+
18
+ def convert_pad_shape(pad_shape):
19
+ l = pad_shape[::-1]
20
+ pad_shape = [item for sublist in l for item in sublist]
21
+ return pad_shape
22
+
23
+
24
+ def intersperse(lst, item):
25
+ result = [item] * (len(lst) * 2 + 1)
26
+ result[1::2] = lst
27
+ return result
28
+
29
+
30
+ def kl_divergence(m_p, logs_p, m_q, logs_q):
31
+ """KL(P||Q)"""
32
+ kl = (logs_q - logs_p) - 0.5
33
+ kl += 0.5 * (torch.exp(2. * logs_p) + ((m_p - m_q)**2)) * torch.exp(-2. * logs_q)
34
+ return kl
35
+
36
+
37
+ def rand_gumbel(shape):
38
+ """Sample from the Gumbel distribution, protect from overflows."""
39
+ uniform_samples = torch.rand(shape) * 0.99998 + 0.00001
40
+ return -torch.log(-torch.log(uniform_samples))
41
+
42
+
43
+ def rand_gumbel_like(x):
44
+ g = rand_gumbel(x.size()).to(dtype=x.dtype, device=x.device)
45
+ return g
46
+
47
+
48
+ def slice_segments(x, ids_str, segment_size=4):
49
+ ret = torch.zeros_like(x[:, :, :segment_size])
50
+ for i in range(x.size(0)):
51
+ idx_str = ids_str[i]
52
+ idx_end = idx_str + segment_size
53
+ ret[i] = x[i, :, idx_str:idx_end]
54
+ return ret
55
+
56
+
57
+ def rand_slice_segments(x, x_lengths=None, segment_size=4):
58
+ b, d, t = x.size()
59
+ if x_lengths is None:
60
+ x_lengths = t
61
+ ids_str_max = x_lengths - segment_size + 1
62
+ ids_str = (torch.rand([b]).to(device=x.device) * ids_str_max).to(dtype=torch.long)
63
+ ret = slice_segments(x, ids_str, segment_size)
64
+ return ret, ids_str
65
+
66
+
67
+ def get_timing_signal_1d(
68
+ length, channels, min_timescale=1.0, max_timescale=1.0e4):
69
+ position = torch.arange(length, dtype=torch.float)
70
+ num_timescales = channels // 2
71
+ log_timescale_increment = (
72
+ math.log(float(max_timescale) / float(min_timescale)) /
73
+ (num_timescales - 1))
74
+ inv_timescales = min_timescale * torch.exp(
75
+ torch.arange(num_timescales, dtype=torch.float) * -log_timescale_increment)
76
+ scaled_time = position.unsqueeze(0) * inv_timescales.unsqueeze(1)
77
+ signal = torch.cat([torch.sin(scaled_time), torch.cos(scaled_time)], 0)
78
+ signal = F.pad(signal, [0, 0, 0, channels % 2])
79
+ signal = signal.view(1, channels, length)
80
+ return signal
81
+
82
+
83
+ def add_timing_signal_1d(x, min_timescale=1.0, max_timescale=1.0e4):
84
+ b, channels, length = x.size()
85
+ signal = get_timing_signal_1d(length, channels, min_timescale, max_timescale)
86
+ return x + signal.to(dtype=x.dtype, device=x.device)
87
+
88
+
89
+ def cat_timing_signal_1d(x, min_timescale=1.0, max_timescale=1.0e4, axis=1):
90
+ b, channels, length = x.size()
91
+ signal = get_timing_signal_1d(length, channels, min_timescale, max_timescale)
92
+ return torch.cat([x, signal.to(dtype=x.dtype, device=x.device)], axis)
93
+
94
+
95
+ def subsequent_mask(length):
96
+ mask = torch.tril(torch.ones(length, length)).unsqueeze(0).unsqueeze(0)
97
+ return mask
98
+
99
+
100
+ @torch.jit.script
101
+ def fused_add_tanh_sigmoid_multiply(input_a, input_b, n_channels):
102
+ n_channels_int = n_channels[0]
103
+ in_act = input_a + input_b
104
+ t_act = torch.tanh(in_act[:, :n_channels_int, :])
105
+ s_act = torch.sigmoid(in_act[:, n_channels_int:, :])
106
+ acts = t_act * s_act
107
+ return acts
108
+
109
+
110
+ def convert_pad_shape(pad_shape):
111
+ l = pad_shape[::-1]
112
+ pad_shape = [item for sublist in l for item in sublist]
113
+ return pad_shape
114
+
115
+
116
+ def shift_1d(x):
117
+ x = F.pad(x, convert_pad_shape([[0, 0], [0, 0], [1, 0]]))[:, :, :-1]
118
+ return x
119
+
120
+
121
+ def sequence_mask(length, max_length=None):
122
+ if max_length is None:
123
+ max_length = length.max()
124
+ x = torch.arange(max_length, dtype=length.dtype, device=length.device)
125
+ return x.unsqueeze(0) < length.unsqueeze(1)
126
+
127
+
128
+ def generate_path(duration, mask):
129
+ """
130
+ duration: [b, 1, t_x]
131
+ mask: [b, 1, t_y, t_x]
132
+ """
133
+ device = duration.device
134
+
135
+ b, _, t_y, t_x = mask.shape
136
+ cum_duration = torch.cumsum(duration, -1)
137
+
138
+ cum_duration_flat = cum_duration.view(b * t_x)
139
+ path = sequence_mask(cum_duration_flat, t_y).to(mask.dtype)
140
+ path = path.view(b, t_x, t_y)
141
+ path = path - F.pad(path, convert_pad_shape([[0, 0], [1, 0], [0, 0]]))[:, :-1]
142
+ path = path.unsqueeze(1).transpose(2,3) * mask
143
+ return path
144
+
145
+
146
+ def clip_grad_value_(parameters, clip_value, norm_type=2):
147
+ if isinstance(parameters, torch.Tensor):
148
+ parameters = [parameters]
149
+ parameters = list(filter(lambda p: p.grad is not None, parameters))
150
+ norm_type = float(norm_type)
151
+ if clip_value is not None:
152
+ clip_value = float(clip_value)
153
+
154
+ total_norm = 0
155
+ for p in parameters:
156
+ param_norm = p.grad.data.norm(norm_type)
157
+ total_norm += param_norm.item() ** norm_type
158
+ if clip_value is not None:
159
+ p.grad.data.clamp_(min=-clip_value, max=clip_value)
160
+ total_norm = total_norm ** (1. / norm_type)
161
+ return total_norm
data_utils.py ADDED
@@ -0,0 +1,392 @@
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
1
+ import time
2
+ import os
3
+ import random
4
+ import numpy as np
5
+ import torch
6
+ import torch.utils.data
7
+
8
+ from ms_istft_vits import commons
9
+ from ms_istft_vits.mel_processing import spectrogram_torch
10
+ from ms_istft_vits.utils import load_wav_to_torch, load_filepaths_and_text
11
+ from ms_istft_vits.text import text_to_sequence, cleaned_text_to_sequence
12
+
13
+
14
+ class TextAudioLoader(torch.utils.data.Dataset):
15
+ """
16
+ 1) loads audio, text pairs
17
+ 2) normalizes text and converts them to sequences of integers
18
+ 3) computes spectrograms from audio files.
19
+ """
20
+ def __init__(self, audiopaths_and_text, hparams):
21
+ self.audiopaths_and_text = load_filepaths_and_text(audiopaths_and_text)
22
+ self.text_cleaners = hparams.text_cleaners
23
+ self.max_wav_value = hparams.max_wav_value
24
+ self.sampling_rate = hparams.sampling_rate
25
+ self.filter_length = hparams.filter_length
26
+ self.hop_length = hparams.hop_length
27
+ self.win_length = hparams.win_length
28
+ self.sampling_rate = hparams.sampling_rate
29
+
30
+ self.cleaned_text = getattr(hparams, "cleaned_text", False)
31
+
32
+ self.add_blank = hparams.add_blank
33
+ self.min_text_len = getattr(hparams, "min_text_len", 1)
34
+ self.max_text_len = getattr(hparams, "max_text_len", 190)
35
+
36
+ random.seed(1234)
37
+ random.shuffle(self.audiopaths_and_text)
38
+ self._filter()
39
+
40
+
41
+ def _filter(self):
42
+ """
43
+ Filter text & store spec lengths
44
+ """
45
+ # Store spectrogram lengths for Bucketing
46
+ # wav_length ~= file_size / (wav_channels * Bytes per dim) = file_size / (1 * 2)
47
+ # spec_length = wav_length // hop_length
48
+
49
+ audiopaths_and_text_new = []
50
+ lengths = []
51
+ for audiopath, text in self.audiopaths_and_text:
52
+ if self.min_text_len <= len(text) and len(text) <= self.max_text_len:
53
+ audiopaths_and_text_new.append([audiopath, text])
54
+ lengths.append(os.path.getsize(audiopath) // (2 * self.hop_length))
55
+ self.audiopaths_and_text = audiopaths_and_text_new
56
+ self.lengths = lengths
57
+
58
+ def get_audio_text_pair(self, audiopath_and_text):
59
+ # separate filename and text
60
+ audiopath, text = audiopath_and_text[0], audiopath_and_text[1]
61
+ text = self.get_text(text)
62
+ spec, wav = self.get_audio(audiopath)
63
+ return (text, spec, wav)
64
+
65
+ def get_audio(self, filename):
66
+ audio, sampling_rate = load_wav_to_torch(filename)
67
+ if sampling_rate != self.sampling_rate:
68
+ raise ValueError("{} {} SR doesn't match target {} SR".format(
69
+ sampling_rate, self.sampling_rate))
70
+ audio_norm = audio / self.max_wav_value
71
+ audio_norm = audio_norm.unsqueeze(0)
72
+ spec_filename = filename.replace(".wav", ".spec.pt")
73
+ if os.path.exists(spec_filename):
74
+ spec = torch.load(spec_filename)
75
+ else:
76
+ spec = spectrogram_torch(audio_norm, self.filter_length,
77
+ self.sampling_rate, self.hop_length, self.win_length,
78
+ center=False)
79
+ spec = torch.squeeze(spec, 0)
80
+ torch.save(spec, spec_filename)
81
+ return spec, audio_norm
82
+
83
+ def get_text(self, text):
84
+ if self.cleaned_text:
85
+ text_norm = cleaned_text_to_sequence(text)
86
+ else:
87
+ text_norm = text_to_sequence(text, self.text_cleaners)
88
+ if self.add_blank:
89
+ text_norm = commons.intersperse(text_norm, 0)
90
+ text_norm = torch.LongTensor(text_norm)
91
+ return text_norm
92
+
93
+ def __getitem__(self, index):
94
+ return self.get_audio_text_pair(self.audiopaths_and_text[index])
95
+
96
+ def __len__(self):
97
+ return len(self.audiopaths_and_text)
98
+
99
+
100
+ class TextAudioCollate():
101
+ """ Zero-pads model inputs and targets
102
+ """
103
+ def __init__(self, return_ids=False):
104
+ self.return_ids = return_ids
105
+
106
+ def __call__(self, batch):
107
+ """Collate's training batch from normalized text and aduio
108
+ PARAMS
109
+ ------
110
+ batch: [text_normalized, spec_normalized, wav_normalized]
111
+ """
112
+ # Right zero-pad all one-hot text sequences to max input length
113
+ _, ids_sorted_decreasing = torch.sort(
114
+ torch.LongTensor([x[1].size(1) for x in batch]),
115
+ dim=0, descending=True)
116
+
117
+ max_text_len = max([len(x[0]) for x in batch])
118
+ max_spec_len = max([x[1].size(1) for x in batch])
119
+ max_wav_len = max([x[2].size(1) for x in batch])
120
+
121
+ text_lengths = torch.LongTensor(len(batch))
122
+ spec_lengths = torch.LongTensor(len(batch))
123
+ wav_lengths = torch.LongTensor(len(batch))
124
+
125
+ text_padded = torch.LongTensor(len(batch), max_text_len)
126
+ spec_padded = torch.FloatTensor(len(batch), batch[0][1].size(0), max_spec_len)
127
+ wav_padded = torch.FloatTensor(len(batch), 1, max_wav_len)
128
+ text_padded.zero_()
129
+ spec_padded.zero_()
130
+ wav_padded.zero_()
131
+ for i in range(len(ids_sorted_decreasing)):
132
+ row = batch[ids_sorted_decreasing[i]]
133
+
134
+ text = row[0]
135
+ text_padded[i, :text.size(0)] = text
136
+ text_lengths[i] = text.size(0)
137
+
138
+ spec = row[1]
139
+ spec_padded[i, :, :spec.size(1)] = spec
140
+ spec_lengths[i] = spec.size(1)
141
+
142
+ wav = row[2]
143
+ wav_padded[i, :, :wav.size(1)] = wav
144
+ wav_lengths[i] = wav.size(1)
145
+
146
+ if self.return_ids:
147
+ return text_padded, text_lengths, spec_padded, spec_lengths, wav_padded, wav_lengths, ids_sorted_decreasing
148
+ return text_padded, text_lengths, spec_padded, spec_lengths, wav_padded, wav_lengths
149
+
150
+
151
+ """Multi speaker version"""
152
+ class TextAudioSpeakerLoader(torch.utils.data.Dataset):
153
+ """
154
+ 1) loads audio, speaker_id, text pairs
155
+ 2) normalizes text and converts them to sequences of integers
156
+ 3) computes spectrograms from audio files.
157
+ """
158
+ def __init__(self, audiopaths_sid_text, hparams):
159
+ self.audiopaths_sid_text = load_filepaths_and_text(audiopaths_sid_text)
160
+ self.text_cleaners = hparams.text_cleaners
161
+ self.max_wav_value = hparams.max_wav_value
162
+ self.sampling_rate = hparams.sampling_rate
163
+ self.filter_length = hparams.filter_length
164
+ self.hop_length = hparams.hop_length
165
+ self.win_length = hparams.win_length
166
+ self.sampling_rate = hparams.sampling_rate
167
+
168
+ self.cleaned_text = getattr(hparams, "cleaned_text", False)
169
+
170
+ self.add_blank = hparams.add_blank
171
+ self.min_text_len = getattr(hparams, "min_text_len", 1)
172
+ self.max_text_len = getattr(hparams, "max_text_len", 190)
173
+
174
+ random.seed(1234)
175
+ random.shuffle(self.audiopaths_sid_text)
176
+ self._filter()
177
+
178
+ def _filter(self):
179
+ """
180
+ Filter text & store spec lengths
181
+ """
182
+ # Store spectrogram lengths for Bucketing
183
+ # wav_length ~= file_size / (wav_channels * Bytes per dim) = file_size / (1 * 2)
184
+ # spec_length = wav_length // hop_length
185
+
186
+ audiopaths_sid_text_new = []
187
+ lengths = []
188
+ for audiopath, sid, text in self.audiopaths_sid_text:
189
+ if self.min_text_len <= len(text) and len(text) <= self.max_text_len:
190
+ audiopaths_sid_text_new.append([audiopath, sid, text])
191
+ lengths.append(os.path.getsize(audiopath) // (2 * self.hop_length))
192
+ self.audiopaths_sid_text = audiopaths_sid_text_new
193
+ self.lengths = lengths
194
+
195
+ def get_audio_text_speaker_pair(self, audiopath_sid_text):
196
+ # separate filename, speaker_id and text
197
+ audiopath, sid, text = audiopath_sid_text[0], audiopath_sid_text[1], audiopath_sid_text[2]
198
+ text = self.get_text(text)
199
+ spec, wav = self.get_audio(audiopath)
200
+ sid = self.get_sid(sid)
201
+ return (text, spec, wav, sid)
202
+
203
+ def get_audio(self, filename):
204
+ audio, sampling_rate = load_wav_to_torch(filename)
205
+ if sampling_rate != self.sampling_rate:
206
+ raise ValueError("{} {} SR doesn't match target {} SR".format(
207
+ sampling_rate, self.sampling_rate))
208
+ audio_norm = audio / self.max_wav_value
209
+ audio_norm = audio_norm.unsqueeze(0)
210
+ spec_filename = filename.replace(".wav", ".spec.pt")
211
+ if os.path.exists(spec_filename):
212
+ spec = torch.load(spec_filename)
213
+ else:
214
+ spec = spectrogram_torch(audio_norm, self.filter_length,
215
+ self.sampling_rate, self.hop_length, self.win_length,
216
+ center=False)
217
+ spec = torch.squeeze(spec, 0)
218
+ torch.save(spec, spec_filename)
219
+ return spec, audio_norm
220
+
221
+ def get_text(self, text):
222
+ if self.cleaned_text:
223
+ text_norm = cleaned_text_to_sequence(text)
224
+ else:
225
+ text_norm = text_to_sequence(text, self.text_cleaners)
226
+ if self.add_blank:
227
+ text_norm = commons.intersperse(text_norm, 0)
228
+ text_norm = torch.LongTensor(text_norm)
229
+ return text_norm
230
+
231
+ def get_sid(self, sid):
232
+ sid = torch.LongTensor([int(sid)])
233
+ return sid
234
+
235
+ def __getitem__(self, index):
236
+ return self.get_audio_text_speaker_pair(self.audiopaths_sid_text[index])
237
+
238
+ def __len__(self):
239
+ return len(self.audiopaths_sid_text)
240
+
241
+
242
+ class TextAudioSpeakerCollate():
243
+ """ Zero-pads model inputs and targets
244
+ """
245
+ def __init__(self, return_ids=False):
246
+ self.return_ids = return_ids
247
+
248
+ def __call__(self, batch):
249
+ """Collate's training batch from normalized text, audio and speaker identities
250
+ PARAMS
251
+ ------
252
+ batch: [text_normalized, spec_normalized, wav_normalized, sid]
253
+ """
254
+ # Right zero-pad all one-hot text sequences to max input length
255
+ _, ids_sorted_decreasing = torch.sort(
256
+ torch.LongTensor([x[1].size(1) for x in batch]),
257
+ dim=0, descending=True)
258
+
259
+ max_text_len = max([len(x[0]) for x in batch])
260
+ max_spec_len = max([x[1].size(1) for x in batch])
261
+ max_wav_len = max([x[2].size(1) for x in batch])
262
+
263
+ text_lengths = torch.LongTensor(len(batch))
264
+ spec_lengths = torch.LongTensor(len(batch))
265
+ wav_lengths = torch.LongTensor(len(batch))
266
+ sid = torch.LongTensor(len(batch))
267
+
268
+ text_padded = torch.LongTensor(len(batch), max_text_len)
269
+ spec_padded = torch.FloatTensor(len(batch), batch[0][1].size(0), max_spec_len)
270
+ wav_padded = torch.FloatTensor(len(batch), 1, max_wav_len)
271
+ text_padded.zero_()
272
+ spec_padded.zero_()
273
+ wav_padded.zero_()
274
+ for i in range(len(ids_sorted_decreasing)):
275
+ row = batch[ids_sorted_decreasing[i]]
276
+
277
+ text = row[0]
278
+ text_padded[i, :text.size(0)] = text
279
+ text_lengths[i] = text.size(0)
280
+
281
+ spec = row[1]
282
+ spec_padded[i, :, :spec.size(1)] = spec
283
+ spec_lengths[i] = spec.size(1)
284
+
285
+ wav = row[2]
286
+ wav_padded[i, :, :wav.size(1)] = wav
287
+ wav_lengths[i] = wav.size(1)
288
+
289
+ sid[i] = row[3]
290
+
291
+ if self.return_ids:
292
+ return text_padded, text_lengths, spec_padded, spec_lengths, wav_padded, wav_lengths, sid, ids_sorted_decreasing
293
+ return text_padded, text_lengths, spec_padded, spec_lengths, wav_padded, wav_lengths, sid
294
+
295
+
296
+ class DistributedBucketSampler(torch.utils.data.distributed.DistributedSampler):
297
+ """
298
+ Maintain similar input lengths in a batch.
299
+ Length groups are specified by boundaries.
300
+ Ex) boundaries = [b1, b2, b3] -> any batch is included either {x | b1 < length(x) <=b2} or {x | b2 < length(x) <= b3}.
301
+
302
+ It removes samples which are not included in the boundaries.
303
+ Ex) boundaries = [b1, b2, b3] -> any x s.t. length(x) <= b1 or length(x) > b3 are discarded.
304
+ """
305
+ def __init__(self, dataset, batch_size, boundaries, num_replicas=None, rank=None, shuffle=True):
306
+ super().__init__(dataset, num_replicas=num_replicas, rank=rank, shuffle=shuffle)
307
+ self.lengths = dataset.lengths
308
+ self.batch_size = batch_size
309
+ self.boundaries = boundaries
310
+
311
+ self.buckets, self.num_samples_per_bucket = self._create_buckets()
312
+ self.total_size = sum(self.num_samples_per_bucket)
313
+ self.num_samples = self.total_size // self.num_replicas
314
+
315
+ def _create_buckets(self):
316
+ buckets = [[] for _ in range(len(self.boundaries) - 1)]
317
+ for i in range(len(self.lengths)):
318
+ length = self.lengths[i]
319
+ idx_bucket = self._bisect(length)
320
+ if idx_bucket != -1:
321
+ buckets[idx_bucket].append(i)
322
+
323
+ for i in range(len(buckets) - 1, 0, -1):
324
+ if len(buckets[i]) == 0:
325
+ buckets.pop(i)
326
+ self.boundaries.pop(i+1)
327
+
328
+ num_samples_per_bucket = []
329
+ for i in range(len(buckets)):
330
+ len_bucket = len(buckets[i])
331
+ total_batch_size = self.num_replicas * self.batch_size
332
+ rem = (total_batch_size - (len_bucket % total_batch_size)) % total_batch_size
333
+ num_samples_per_bucket.append(len_bucket + rem)
334
+ return buckets, num_samples_per_bucket
335
+
336
+ def __iter__(self):
337
+ # deterministically shuffle based on epoch
338
+ g = torch.Generator()
339
+ g.manual_seed(self.epoch)
340
+
341
+ indices = []
342
+ if self.shuffle:
343
+ for bucket in self.buckets:
344
+ indices.append(torch.randperm(len(bucket), generator=g).tolist())
345
+ else:
346
+ for bucket in self.buckets:
347
+ indices.append(list(range(len(bucket))))
348
+
349
+ batches = []
350
+ for i in range(len(self.buckets)):
351
+ bucket = self.buckets[i]
352
+ len_bucket = len(bucket)
353
+ ids_bucket = indices[i]
354
+ num_samples_bucket = self.num_samples_per_bucket[i]
355
+
356
+ # add extra samples to make it evenly divisible
357
+ rem = num_samples_bucket - len_bucket
358
+ ids_bucket = ids_bucket + ids_bucket * (rem // len_bucket) + ids_bucket[:(rem % len_bucket)]
359
+
360
+ # subsample
361
+ ids_bucket = ids_bucket[self.rank::self.num_replicas]
362
+
363
+ # batching
364
+ for j in range(len(ids_bucket) // self.batch_size):
365
+ batch = [bucket[idx] for idx in ids_bucket[j*self.batch_size:(j+1)*self.batch_size]]
366
+ batches.append(batch)
367
+
368
+ if self.shuffle:
369
+ batch_ids = torch.randperm(len(batches), generator=g).tolist()
370
+ batches = [batches[i] for i in batch_ids]
371
+ self.batches = batches
372
+
373
+ assert len(self.batches) * self.batch_size == self.num_samples
374
+ return iter(self.batches)
375
+
376
+ def _bisect(self, x, lo=0, hi=None):
377
+ if hi is None:
378
+ hi = len(self.boundaries) - 1
379
+
380
+ if hi > lo:
381
+ mid = (hi + lo) // 2
382
+ if self.boundaries[mid] < x and x <= self.boundaries[mid+1]:
383
+ return mid
384
+ elif x <= self.boundaries[mid]:
385
+ return self._bisect(x, lo, mid)
386
+ else:
387
+ return self._bisect(x, mid + 1, hi)
388
+ else:
389
+ return -1
390
+
391
+ def __len__(self):
392
+ return self.num_samples // self.batch_size
inference.py ADDED
@@ -0,0 +1,37 @@
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
1
+ # import sengiri
2
+ import re
3
+ import torch
4
+ import commons
5
+ import utils
6
+ from models import SynthesizerTrn
7
+ from text.symbols import symbols
8
+ from text import text_to_sequence
9
+
10
+
11
+ def get_text(text, hps):
12
+ text_norm = text_to_sequence(text, hps.data.text_cleaners)
13
+ if hps.data.add_blank:
14
+ text_norm = commons.intersperse(text_norm, 0)
15
+ text_norm = torch.LongTensor(text_norm)
16
+ return text_norm
17
+
18
+
19
+ hps = utils.get_hparams_from_file(f"pretrained_model/arona_ms_istft_vits_config.json")
20
+ net_g = SynthesizerTrn(
21
+ len(symbols),
22
+ hps.data.filter_length // 2 + 1,
23
+ hps.train.segment_size // hps.data.hop_length,
24
+ n_speakers=hps.data.n_speakers,
25
+ **hps.model).cuda()
26
+ _ = net_g.eval()
27
+
28
+ _ = utils.load_checkpoint(f"pretrained_model/arona_ms_istft_vits.pth", net_g, None)
29
+
30
+ text = '物語は嘘から始まる。'
31
+ SPEECH_SPEED = 1
32
+ stn_tst = get_text(text, hps)
33
+ with torch.no_grad():
34
+ x_tst = stn_tst.cuda().unsqueeze(0)
35
+ x_tst_lengths = torch.LongTensor([stn_tst.size(0)]).cuda()
36
+ audio = net_g.infer(x_tst, x_tst_lengths, noise_scale=.667, noise_scale_w=0.8, length_scale=1/SPEECH_SPEED)[0][0,0].data.cpu().float().numpy()
37
+
losses.py ADDED
@@ -0,0 +1,71 @@
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
1
+ import torch
2
+ from torch.nn import functional as F
3
+ from ms_istft_vits.stft_loss import MultiResolutionSTFTLoss
4
+
5
+
6
+ from ms_istft_vits import commons
7
+
8
+
9
+ def feature_loss(fmap_r, fmap_g):
10
+ loss = 0
11
+ for dr, dg in zip(fmap_r, fmap_g):
12
+ for rl, gl in zip(dr, dg):
13
+ rl = rl.float().detach()
14
+ gl = gl.float()
15
+ loss += torch.mean(torch.abs(rl - gl))
16
+
17
+ return loss * 2
18
+
19
+
20
+ def discriminator_loss(disc_real_outputs, disc_generated_outputs):
21
+ loss = 0
22
+ r_losses = []
23
+ g_losses = []
24
+ for dr, dg in zip(disc_real_outputs, disc_generated_outputs):
25
+ dr = dr.float()
26
+ dg = dg.float()
27
+ r_loss = torch.mean((1-dr)**2)
28
+ g_loss = torch.mean(dg**2)
29
+ loss += (r_loss + g_loss)
30
+ r_losses.append(r_loss.item())
31
+ g_losses.append(g_loss.item())
32
+
33
+ return loss, r_losses, g_losses
34
+
35
+
36
+ def generator_loss(disc_outputs):
37
+ loss = 0
38
+ gen_losses = []
39
+ for dg in disc_outputs:
40
+ dg = dg.float()
41
+ l = torch.mean((1-dg)**2)
42
+ gen_losses.append(l)
43
+ loss += l
44
+
45
+ return loss, gen_losses
46
+
47
+
48
+ def kl_loss(z_p, logs_q, m_p, logs_p, z_mask):
49
+ """
50
+ z_p, logs_q: [b, h, t_t]
51
+ m_p, logs_p: [b, h, t_t]
52
+ """
53
+ z_p = z_p.float()
54
+ logs_q = logs_q.float()
55
+ m_p = m_p.float()
56
+ logs_p = logs_p.float()
57
+ z_mask = z_mask.float()
58
+
59
+ kl = logs_p - logs_q - 0.5
60
+ kl += 0.5 * ((z_p - m_p)**2) * torch.exp(-2. * logs_p)
61
+ kl = torch.sum(kl * z_mask)
62
+ l = kl / torch.sum(z_mask)
63
+ return l
64
+
65
+ def subband_stft_loss(h, y_mb, y_hat_mb):
66
+ sub_stft_loss = MultiResolutionSTFTLoss(h.train.fft_sizes, h.train.hop_sizes, h.train.win_lengths)
67
+ y_mb = y_mb.view(-1, y_mb.size(2))
68
+ y_hat_mb = y_hat_mb.view(-1, y_hat_mb.size(2))
69
+ sub_sc_loss, sub_mag_loss = sub_stft_loss(y_hat_mb[:, :y_mb.size(-1)], y_mb)
70
+ return sub_sc_loss+sub_mag_loss
71
+
mel_processing.py ADDED
@@ -0,0 +1,112 @@
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
1
+ import math
2
+ import os
3
+ import random
4
+ import torch
5
+ from torch import nn
6
+ import torch.nn.functional as F
7
+ import torch.utils.data
8
+ import numpy as np
9
+ import librosa
10
+ import librosa.util as librosa_util
11
+ from librosa.util import normalize, pad_center, tiny
12
+ from scipy.signal import get_window
13
+ from scipy.io.wavfile import read
14
+ from librosa.filters import mel as librosa_mel_fn
15
+
16
+ MAX_WAV_VALUE = 32768.0
17
+
18
+
19
+ def dynamic_range_compression_torch(x, C=1, clip_val=1e-5):
20
+ """
21
+ PARAMS
22
+ ------
23
+ C: compression factor
24
+ """
25
+ return torch.log(torch.clamp(x, min=clip_val) * C)
26
+
27
+
28
+ def dynamic_range_decompression_torch(x, C=1):
29
+ """
30
+ PARAMS
31
+ ------
32
+ C: compression factor used to compress
33
+ """
34
+ return torch.exp(x) / C
35
+
36
+
37
+ def spectral_normalize_torch(magnitudes):
38
+ output = dynamic_range_compression_torch(magnitudes)
39
+ return output
40
+
41
+
42
+ def spectral_de_normalize_torch(magnitudes):
43
+ output = dynamic_range_decompression_torch(magnitudes)
44
+ return output
45
+
46
+
47
+ mel_basis = {}
48
+ hann_window = {}
49
+
50
+
51
+ def spectrogram_torch(y, n_fft, sampling_rate, hop_size, win_size, center=False):
52
+ if torch.min(y) < -1.:
53
+ print('min value is ', torch.min(y))
54
+ if torch.max(y) > 1.:
55
+ print('max value is ', torch.max(y))
56
+
57
+ global hann_window
58
+ dtype_device = str(y.dtype) + '_' + str(y.device)
59
+ wnsize_dtype_device = str(win_size) + '_' + dtype_device
60
+ if wnsize_dtype_device not in hann_window:
61
+ hann_window[wnsize_dtype_device] = torch.hann_window(win_size).to(dtype=y.dtype, device=y.device)
62
+
63
+ y = torch.nn.functional.pad(y.unsqueeze(1), (int((n_fft-hop_size)/2), int((n_fft-hop_size)/2)), mode='reflect')
64
+ y = y.squeeze(1)
65
+
66
+ spec = torch.stft(y, n_fft, hop_length=hop_size, win_length=win_size, window=hann_window[wnsize_dtype_device],
67
+ center=center, pad_mode='reflect', normalized=False, onesided=True)
68
+
69
+ spec = torch.sqrt(spec.pow(2).sum(-1) + 1e-6)
70
+ return spec
71
+
72
+
73
+ def spec_to_mel_torch(spec, n_fft, num_mels, sampling_rate, fmin, fmax):
74
+ global mel_basis
75
+ dtype_device = str(spec.dtype) + '_' + str(spec.device)
76
+ fmax_dtype_device = str(fmax) + '_' + dtype_device
77
+ if fmax_dtype_device not in mel_basis:
78
+ mel = librosa_mel_fn(sampling_rate, n_fft, num_mels, fmin, fmax)
79
+ mel_basis[fmax_dtype_device] = torch.from_numpy(mel).to(dtype=spec.dtype, device=spec.device)
80
+ spec = torch.matmul(mel_basis[fmax_dtype_device], spec)
81
+ spec = spectral_normalize_torch(spec)
82
+ return spec
83
+
84
+
85
+ def mel_spectrogram_torch(y, n_fft, num_mels, sampling_rate, hop_size, win_size, fmin, fmax, center=False):
86
+ if torch.min(y) < -1.:
87
+ print('min value is ', torch.min(y))
88
+ if torch.max(y) > 1.:
89
+ print('max value is ', torch.max(y))
90
+
91
+ global mel_basis, hann_window
92
+ dtype_device = str(y.dtype) + '_' + str(y.device)
93
+ fmax_dtype_device = str(fmax) + '_' + dtype_device
94
+ wnsize_dtype_device = str(win_size) + '_' + dtype_device
95
+ if fmax_dtype_device not in mel_basis:
96
+ mel = librosa_mel_fn(sampling_rate, n_fft, num_mels, fmin, fmax)
97
+ mel_basis[fmax_dtype_device] = torch.from_numpy(mel).to(dtype=y.dtype, device=y.device)
98
+ if wnsize_dtype_device not in hann_window:
99
+ hann_window[wnsize_dtype_device] = torch.hann_window(win_size).to(dtype=y.dtype, device=y.device)
100
+
101
+ y = torch.nn.functional.pad(y.unsqueeze(1), (int((n_fft-hop_size)/2), int((n_fft-hop_size)/2)), mode='reflect')
102
+ y = y.squeeze(1)
103
+
104
+ spec = torch.stft(y, n_fft, hop_length=hop_size, win_length=win_size, window=hann_window[wnsize_dtype_device],
105
+ center=center, pad_mode='reflect', normalized=False, onesided=True)
106
+
107
+ spec = torch.sqrt(spec.pow(2).sum(-1) + 1e-6)
108
+
109
+ spec = torch.matmul(mel_basis[fmax_dtype_device], spec)
110
+ spec = spectral_normalize_torch(spec)
111
+
112
+ return spec
models.py ADDED
@@ -0,0 +1,731 @@
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
1
+ import copy
2
+ import math
3
+ import torch
4
+ from torch import nn
5
+ from torch.nn import functional as F
6
+
7
+ import commons
8
+ import modules
9
+ import attentions
10
+ import monotonic_align
11
+
12
+ from torch.nn import Conv1d, ConvTranspose1d, AvgPool1d, Conv2d
13
+ from torch.nn.utils import weight_norm, remove_weight_norm, spectral_norm
14
+ from commons import init_weights, get_padding
15
+ from pqmf import PQMF
16
+ from stft import TorchSTFT
17
+ import math
18
+
19
+
20
+ class StochasticDurationPredictor(nn.Module):
21
+ def __init__(self, in_channels, filter_channels, kernel_size, p_dropout, n_flows=4, gin_channels=0):
22
+ super().__init__()
23
+ filter_channels = in_channels # it needs to be removed from future version.
24
+ self.in_channels = in_channels
25
+ self.filter_channels = filter_channels
26
+ self.kernel_size = kernel_size
27
+ self.p_dropout = p_dropout
28
+ self.n_flows = n_flows
29
+ self.gin_channels = gin_channels
30
+
31
+ self.log_flow = modules.Log()
32
+ self.flows = nn.ModuleList()
33
+ self.flows.append(modules.ElementwiseAffine(2))
34
+ for i in range(n_flows):
35
+ self.flows.append(modules.ConvFlow(2, filter_channels, kernel_size, n_layers=3))
36
+ self.flows.append(modules.Flip())
37
+
38
+ self.post_pre = nn.Conv1d(1, filter_channels, 1)
39
+ self.post_proj = nn.Conv1d(filter_channels, filter_channels, 1)
40
+ self.post_convs = modules.DDSConv(filter_channels, kernel_size, n_layers=3, p_dropout=p_dropout)
41
+ self.post_flows = nn.ModuleList()
42
+ self.post_flows.append(modules.ElementwiseAffine(2))
43
+ for i in range(4):
44
+ self.post_flows.append(modules.ConvFlow(2, filter_channels, kernel_size, n_layers=3))
45
+ self.post_flows.append(modules.Flip())
46
+
47
+ self.pre = nn.Conv1d(in_channels, filter_channels, 1)
48
+ self.proj = nn.Conv1d(filter_channels, filter_channels, 1)
49
+ self.convs = modules.DDSConv(filter_channels, kernel_size, n_layers=3, p_dropout=p_dropout)
50
+ if gin_channels != 0:
51
+ self.cond = nn.Conv1d(gin_channels, filter_channels, 1)
52
+
53
+ def forward(self, x, x_mask, w=None, g=None, reverse=False, noise_scale=1.0):
54
+ x = torch.detach(x)
55
+ x = self.pre(x)
56
+ if g is not None:
57
+ g = torch.detach(g)
58
+ x = x + self.cond(g)
59
+ x = self.convs(x, x_mask)
60
+ x = self.proj(x) * x_mask
61
+
62
+ if not reverse:
63
+ flows = self.flows
64
+ assert w is not None
65
+
66
+ logdet_tot_q = 0
67
+ h_w = self.post_pre(w)
68
+ h_w = self.post_convs(h_w, x_mask)
69
+ h_w = self.post_proj(h_w) * x_mask
70
+ e_q = torch.randn(w.size(0), 2, w.size(2)).to(device=x.device, dtype=x.dtype) * x_mask
71
+ z_q = e_q
72
+ for flow in self.post_flows:
73
+ z_q, logdet_q = flow(z_q, x_mask, g=(x + h_w))
74
+ logdet_tot_q += logdet_q
75
+ z_u, z1 = torch.split(z_q, [1, 1], 1)
76
+ u = torch.sigmoid(z_u) * x_mask
77
+ z0 = (w - u) * x_mask
78
+ logdet_tot_q += torch.sum((F.logsigmoid(z_u) + F.logsigmoid(-z_u)) * x_mask, [1,2])
79
+ logq = torch.sum(-0.5 * (math.log(2*math.pi) + (e_q**2)) * x_mask, [1,2]) - logdet_tot_q
80
+
81
+ logdet_tot = 0
82
+ z0, logdet = self.log_flow(z0, x_mask)
83
+ logdet_tot += logdet
84
+ z = torch.cat([z0, z1], 1)
85
+ for flow in flows:
86
+ z, logdet = flow(z, x_mask, g=x, reverse=reverse)
87
+ logdet_tot = logdet_tot + logdet
88
+ nll = torch.sum(0.5 * (math.log(2*math.pi) + (z**2)) * x_mask, [1,2]) - logdet_tot
89
+ return nll + logq # [b]
90
+ else:
91
+ flows = list(reversed(self.flows))
92
+ flows = flows[:-2] + [flows[-1]] # remove a useless vflow
93
+ z = torch.randn(x.size(0), 2, x.size(2)).to(device=x.device, dtype=x.dtype) * noise_scale
94
+ for flow in flows:
95
+ z = flow(z, x_mask, g=x, reverse=reverse)
96
+ z0, z1 = torch.split(z, [1, 1], 1)
97
+ logw = z0
98
+ return logw
99
+
100
+
101
+ class DurationPredictor(nn.Module):
102
+ def __init__(self, in_channels, filter_channels, kernel_size, p_dropout, gin_channels=0):
103
+ super().__init__()
104
+
105
+ self.in_channels = in_channels
106
+ self.filter_channels = filter_channels
107
+ self.kernel_size = kernel_size
108
+ self.p_dropout = p_dropout
109
+ self.gin_channels = gin_channels
110
+
111
+ self.drop = nn.Dropout(p_dropout)
112
+ self.conv_1 = nn.Conv1d(in_channels, filter_channels, kernel_size, padding=kernel_size//2)
113
+ self.norm_1 = modules.LayerNorm(filter_channels)
114
+ self.conv_2 = nn.Conv1d(filter_channels, filter_channels, kernel_size, padding=kernel_size//2)
115
+ self.norm_2 = modules.LayerNorm(filter_channels)
116
+ self.proj = nn.Conv1d(filter_channels, 1, 1)
117
+
118
+ if gin_channels != 0:
119
+ self.cond = nn.Conv1d(gin_channels, in_channels, 1)
120
+
121
+ def forward(self, x, x_mask, g=None):
122
+ x = torch.detach(x)
123
+ if g is not None:
124
+ g = torch.detach(g)
125
+ x = x + self.cond(g)
126
+ x = self.conv_1(x * x_mask)
127
+ x = torch.relu(x)
128
+ x = self.norm_1(x)
129
+ x = self.drop(x)
130
+ x = self.conv_2(x * x_mask)
131
+ x = torch.relu(x)
132
+ x = self.norm_2(x)
133
+ x = self.drop(x)
134
+ x = self.proj(x * x_mask)
135
+ return x * x_mask
136
+
137
+
138
+ class TextEncoder(nn.Module):
139
+ def __init__(self,
140
+ n_vocab,
141
+ out_channels,
142
+ hidden_channels,
143
+ filter_channels,
144
+ n_heads,
145
+ n_layers,
146
+ kernel_size,
147
+ p_dropout):
148
+ super().__init__()
149
+ self.n_vocab = n_vocab
150
+ self.out_channels = out_channels
151
+ self.hidden_channels = hidden_channels
152
+ self.filter_channels = filter_channels
153
+ self.n_heads = n_heads
154
+ self.n_layers = n_layers
155
+ self.kernel_size = kernel_size
156
+ self.p_dropout = p_dropout
157
+
158
+ self.emb = nn.Embedding(n_vocab, hidden_channels)
159
+ nn.init.normal_(self.emb.weight, 0.0, hidden_channels**-0.5)
160
+
161
+ self.encoder = attentions.Encoder(
162
+ hidden_channels,
163
+ filter_channels,
164
+ n_heads,
165
+ n_layers,
166
+ kernel_size,
167
+ p_dropout)
168
+ self.proj= nn.Conv1d(hidden_channels, out_channels * 2, 1)
169
+
170
+ def forward(self, x, x_lengths):
171
+ x = self.emb(x) * math.sqrt(self.hidden_channels) # [b, t, h]
172
+ x = torch.transpose(x, 1, -1) # [b, h, t]
173
+ x_mask = torch.unsqueeze(commons.sequence_mask(x_lengths, x.size(2)), 1).to(x.dtype)
174
+
175
+ x = self.encoder(x * x_mask, x_mask)
176
+ stats = self.proj(x) * x_mask
177
+
178
+ m, logs = torch.split(stats, self.out_channels, dim=1)
179
+ return x, m, logs, x_mask
180
+
181
+
182
+ class ResidualCouplingBlock(nn.Module):
183
+ def __init__(self,
184
+ channels,
185
+ hidden_channels,
186
+ kernel_size,
187
+ dilation_rate,
188
+ n_layers,
189
+ n_flows=4,
190
+ gin_channels=0):
191
+ super().__init__()
192
+ self.channels = channels
193
+ self.hidden_channels = hidden_channels
194
+ self.kernel_size = kernel_size
195
+ self.dilation_rate = dilation_rate
196
+ self.n_layers = n_layers
197
+ self.n_flows = n_flows
198
+ self.gin_channels = gin_channels
199
+
200
+ self.flows = nn.ModuleList()
201
+ for i in range(n_flows):
202
+ self.flows.append(modules.ResidualCouplingLayer(channels, hidden_channels, kernel_size, dilation_rate, n_layers, gin_channels=gin_channels, mean_only=True))
203
+ self.flows.append(modules.Flip())
204
+
205
+ def forward(self, x, x_mask, g=None, reverse=False):
206
+ if not reverse:
207
+ for flow in self.flows:
208
+ x, _ = flow(x, x_mask, g=g, reverse=reverse)
209
+ else:
210
+ for flow in reversed(self.flows):
211
+ x = flow(x, x_mask, g=g, reverse=reverse)
212
+ return x
213
+
214
+
215
+ class PosteriorEncoder(nn.Module):
216
+ def __init__(self,
217
+ in_channels,
218
+ out_channels,
219
+ hidden_channels,
220
+ kernel_size,
221
+ dilation_rate,
222
+ n_layers,
223
+ gin_channels=0):
224
+ super().__init__()
225
+ self.in_channels = in_channels
226
+ self.out_channels = out_channels
227
+ self.hidden_channels = hidden_channels
228
+ self.kernel_size = kernel_size
229
+ self.dilation_rate = dilation_rate
230
+ self.n_layers = n_layers
231
+ self.gin_channels = gin_channels
232
+
233
+ self.pre = nn.Conv1d(in_channels, hidden_channels, 1)
234
+ self.enc = modules.WN(hidden_channels, kernel_size, dilation_rate, n_layers, gin_channels=gin_channels)
235
+ self.proj = nn.Conv1d(hidden_channels, out_channels * 2, 1)
236
+
237
+ def forward(self, x, x_lengths, g=None):
238
+ x_mask = torch.unsqueeze(commons.sequence_mask(x_lengths, x.size(2)), 1).to(x.dtype)
239
+ x = self.pre(x) * x_mask
240
+ x = self.enc(x, x_mask, g=g)
241
+ stats = self.proj(x) * x_mask
242
+ m, logs = torch.split(stats, self.out_channels, dim=1)
243
+ z = (m + torch.randn_like(m) * torch.exp(logs)) * x_mask
244
+ return z, m, logs, x_mask
245
+
246
+ class iSTFT_Generator(torch.nn.Module):
247
+ def __init__(self, initial_channel, resblock, resblock_kernel_sizes, resblock_dilation_sizes, upsample_rates, upsample_initial_channel, upsample_kernel_sizes, gen_istft_n_fft, gen_istft_hop_size, gin_channels=0):
248
+ super(iSTFT_Generator, self).__init__()
249
+ # self.h = h
250
+ self.gen_istft_n_fft = gen_istft_n_fft
251
+ self.gen_istft_hop_size = gen_istft_hop_size
252
+
253
+ self.num_kernels = len(resblock_kernel_sizes)
254
+ self.num_upsamples = len(upsample_rates)
255
+ self.conv_pre = weight_norm(Conv1d(initial_channel, upsample_initial_channel, 7, 1, padding=3))
256
+ resblock = modules.ResBlock1 if resblock == '1' else modules.ResBlock2
257
+
258
+ self.ups = nn.ModuleList()
259
+ for i, (u, k) in enumerate(zip(upsample_rates, upsample_kernel_sizes)):
260
+ self.ups.append(weight_norm(
261
+ ConvTranspose1d(upsample_initial_channel//(2**i), upsample_initial_channel//(2**(i+1)),
262
+ k, u, padding=(k-u)//2)))
263
+
264
+ self.resblocks = nn.ModuleList()
265
+ for i in range(len(self.ups)):
266
+ ch = upsample_initial_channel//(2**(i+1))
267
+ for j, (k, d) in enumerate(zip(resblock_kernel_sizes, resblock_dilation_sizes)):
268
+ self.resblocks.append(resblock(ch, k, d))
269
+
270
+ self.post_n_fft = self.gen_istft_n_fft
271
+ self.conv_post = weight_norm(Conv1d(ch, self.post_n_fft + 2, 7, 1, padding=3))
272
+ self.ups.apply(init_weights)
273
+ self.conv_post.apply(init_weights)
274
+ self.reflection_pad = torch.nn.ReflectionPad1d((1, 0))
275
+ self.stft = TorchSTFT(filter_length=self.gen_istft_n_fft, hop_length=self.gen_istft_hop_size, win_length=self.gen_istft_n_fft)
276
+ def forward(self, x, g=None):
277
+
278
+ x = self.conv_pre(x)
279
+ for i in range(self.num_upsamples):
280
+ x = F.leaky_relu(x, modules.LRELU_SLOPE)
281
+ x = self.ups[i](x)
282
+ xs = None
283
+ for j in range(self.num_kernels):
284
+ if xs is None:
285
+ xs = self.resblocks[i*self.num_kernels+j](x)
286
+ else:
287
+ xs += self.resblocks[i*self.num_kernels+j](x)
288
+ x = xs / self.num_kernels
289
+ x = F.leaky_relu(x)
290
+ x = self.reflection_pad(x)
291
+ x = self.conv_post(x)
292
+ spec = torch.exp(x[:,:self.post_n_fft // 2 + 1, :])
293
+ phase = math.pi*torch.sin(x[:, self.post_n_fft // 2 + 1:, :])
294
+ out = self.stft.inverse(spec, phase).to(x.device)
295
+ return out, None
296
+
297
+ def remove_weight_norm(self):
298
+ print('Removing weight norm...')
299
+ for l in self.ups:
300
+ remove_weight_norm(l)
301
+ for l in self.resblocks:
302
+ l.remove_weight_norm()
303
+ remove_weight_norm(self.conv_pre)
304
+ remove_weight_norm(self.conv_post)
305
+
306
+
307
+ class Multiband_iSTFT_Generator(torch.nn.Module):
308
+ def __init__(self, initial_channel, resblock, resblock_kernel_sizes, resblock_dilation_sizes, upsample_rates, upsample_initial_channel, upsample_kernel_sizes, gen_istft_n_fft, gen_istft_hop_size, subbands, gin_channels=0):
309
+ super(Multiband_iSTFT_Generator, self).__init__()
310
+ # self.h = h
311
+ self.subbands = subbands
312
+ self.num_kernels = len(resblock_kernel_sizes)
313
+ self.num_upsamples = len(upsample_rates)
314
+ self.conv_pre = weight_norm(Conv1d(initial_channel, upsample_initial_channel, 7, 1, padding=3))
315
+ resblock = modules.ResBlock1 if resblock == '1' else modules.ResBlock2
316
+
317
+ self.ups = nn.ModuleList()
318
+ for i, (u, k) in enumerate(zip(upsample_rates, upsample_kernel_sizes)):
319
+ self.ups.append(weight_norm(
320
+ ConvTranspose1d(upsample_initial_channel//(2**i), upsample_initial_channel//(2**(i+1)),
321
+ k, u, padding=(k-u)//2)))
322
+
323
+ self.resblocks = nn.ModuleList()
324
+ for i in range(len(self.ups)):
325
+ ch = upsample_initial_channel//(2**(i+1))
326
+ for j, (k, d) in enumerate(zip(resblock_kernel_sizes, resblock_dilation_sizes)):
327
+ self.resblocks.append(resblock(ch, k, d))
328
+
329
+ self.post_n_fft = gen_istft_n_fft
330
+ self.ups.apply(init_weights)
331
+ self.reflection_pad = torch.nn.ReflectionPad1d((1, 0))
332
+ self.reshape_pixelshuffle = []
333
+
334
+ self.subband_conv_post = weight_norm(Conv1d(ch, self.subbands*(self.post_n_fft + 2), 7, 1, padding=3))
335
+
336
+ self.subband_conv_post.apply(init_weights)
337
+
338
+ self.gen_istft_n_fft = gen_istft_n_fft
339
+ self.gen_istft_hop_size = gen_istft_hop_size
340
+
341
+
342
+ def forward(self, x, g=None):
343
+ stft = TorchSTFT(filter_length=self.gen_istft_n_fft, hop_length=self.gen_istft_hop_size, win_length=self.gen_istft_n_fft).to(x.device)
344
+ pqmf = PQMF(x.device)
345
+
346
+ x = self.conv_pre(x)#[B, ch, length]
347
+
348
+ for i in range(self.num_upsamples):
349
+ x = F.leaky_relu(x, modules.LRELU_SLOPE)
350
+ x = self.ups[i](x)
351
+
352
+
353
+ xs = None
354
+ for j in range(self.num_kernels):
355
+ if xs is None:
356
+ xs = self.resblocks[i*self.num_kernels+j](x)
357
+ else:
358
+ xs += self.resblocks[i*self.num_kernels+j](x)
359
+ x = xs / self.num_kernels
360
+
361
+ x = F.leaky_relu(x)
362
+ x = self.reflection_pad(x)
363
+ x = self.subband_conv_post(x)
364
+ x = torch.reshape(x, (x.shape[0], self.subbands, x.shape[1]//self.subbands, x.shape[-1]))
365
+
366
+ spec = torch.exp(x[:,:,:self.post_n_fft // 2 + 1, :])
367
+ phase = math.pi*torch.sin(x[:,:, self.post_n_fft // 2 + 1:, :])
368
+
369
+ y_mb_hat = stft.inverse(torch.reshape(spec, (spec.shape[0]*self.subbands, self.gen_istft_n_fft // 2 + 1, spec.shape[-1])), torch.reshape(phase, (phase.shape[0]*self.subbands, self.gen_istft_n_fft // 2 + 1, phase.shape[-1])))
370
+ y_mb_hat = torch.reshape(y_mb_hat, (x.shape[0], self.subbands, 1, y_mb_hat.shape[-1]))
371
+ y_mb_hat = y_mb_hat.squeeze(-2)
372
+
373
+ y_g_hat = pqmf.synthesis(y_mb_hat)
374
+
375
+ return y_g_hat, y_mb_hat
376
+
377
+ def remove_weight_norm(self):
378
+ print('Removing weight norm...')
379
+ for l in self.ups:
380
+ remove_weight_norm(l)
381
+ for l in self.resblocks:
382
+ l.remove_weight_norm()
383
+
384
+
385
+ class Multistream_iSTFT_Generator(torch.nn.Module):
386
+ def __init__(self, initial_channel, resblock, resblock_kernel_sizes, resblock_dilation_sizes, upsample_rates, upsample_initial_channel, upsample_kernel_sizes, gen_istft_n_fft, gen_istft_hop_size, subbands, gin_channels=0):
387
+ super(Multistream_iSTFT_Generator, self).__init__()
388
+ # self.h = h
389
+ self.subbands = subbands
390
+ self.num_kernels = len(resblock_kernel_sizes)
391
+ self.num_upsamples = len(upsample_rates)
392
+ self.conv_pre = weight_norm(Conv1d(initial_channel, upsample_initial_channel, 7, 1, padding=3))
393
+ resblock = modules.ResBlock1 if resblock == '1' else modules.ResBlock2
394
+
395
+ self.ups = nn.ModuleList()
396
+ for i, (u, k) in enumerate(zip(upsample_rates, upsample_kernel_sizes)):
397
+ self.ups.append(weight_norm(
398
+ ConvTranspose1d(upsample_initial_channel//(2**i), upsample_initial_channel//(2**(i+1)),
399
+ k, u, padding=(k-u)//2)))
400
+
401
+ self.resblocks = nn.ModuleList()
402
+ for i in range(len(self.ups)):
403
+ ch = upsample_initial_channel//(2**(i+1))
404
+ for j, (k, d) in enumerate(zip(resblock_kernel_sizes, resblock_dilation_sizes)):
405
+ self.resblocks.append(resblock(ch, k, d))
406
+
407
+ self.post_n_fft = gen_istft_n_fft
408
+ self.ups.apply(init_weights)
409
+ self.reflection_pad = torch.nn.ReflectionPad1d((1, 0))
410
+ self.reshape_pixelshuffle = []
411
+
412
+ self.subband_conv_post = weight_norm(Conv1d(ch, self.subbands*(self.post_n_fft + 2), 7, 1, padding=3))
413
+
414
+ self.subband_conv_post.apply(init_weights)
415
+
416
+ self.gen_istft_n_fft = gen_istft_n_fft
417
+ self.gen_istft_hop_size = gen_istft_hop_size
418
+
419
+ updown_filter = torch.zeros((self.subbands, self.subbands, self.subbands)).float()
420
+ for k in range(self.subbands):
421
+ updown_filter[k, k, 0] = 1.0
422
+ self.register_buffer("updown_filter", updown_filter)
423
+ self.multistream_conv_post = weight_norm(Conv1d(4, 1, kernel_size=63, bias=False, padding=get_padding(63, 1)))
424
+ self.multistream_conv_post.apply(init_weights)
425
+
426
+
427
+
428
+ def forward(self, x, g=None):
429
+ stft = TorchSTFT(filter_length=self.gen_istft_n_fft, hop_length=self.gen_istft_hop_size, win_length=self.gen_istft_n_fft).to(x.device)
430
+ # pqmf = PQMF(x.device)
431
+
432
+ x = self.conv_pre(x)#[B, ch, length]
433
+
434
+ for i in range(self.num_upsamples):
435
+
436
+
437
+ x = F.leaky_relu(x, modules.LRELU_SLOPE)
438
+ x = self.ups[i](x)
439
+
440
+
441
+ xs = None
442
+ for j in range(self.num_kernels):
443
+ if xs is None:
444
+ xs = self.resblocks[i*self.num_kernels+j](x)
445
+ else:
446
+ xs += self.resblocks[i*self.num_kernels+j](x)
447
+ x = xs / self.num_kernels
448
+
449
+ x = F.leaky_relu(x)
450
+ x = self.reflection_pad(x)
451
+ x = self.subband_conv_post(x)
452
+ x = torch.reshape(x, (x.shape[0], self.subbands, x.shape[1]//self.subbands, x.shape[-1]))
453
+
454
+ spec = torch.exp(x[:,:,:self.post_n_fft // 2 + 1, :])
455
+ phase = math.pi*torch.sin(x[:,:, self.post_n_fft // 2 + 1:, :])
456
+
457
+ y_mb_hat = stft.inverse(torch.reshape(spec, (spec.shape[0]*self.subbands, self.gen_istft_n_fft // 2 + 1, spec.shape[-1])), torch.reshape(phase, (phase.shape[0]*self.subbands, self.gen_istft_n_fft // 2 + 1, phase.shape[-1])))
458
+ y_mb_hat = torch.reshape(y_mb_hat, (x.shape[0], self.subbands, 1, y_mb_hat.shape[-1]))
459
+ y_mb_hat = y_mb_hat.squeeze(-2)
460
+
461
+ # y_mb_hat = F.conv_transpose1d(y_mb_hat, self.updown_filter.cuda(x.device) * self.subbands, stride=self.subbands)
462
+ y_mb_hat = F.conv_transpose1d(y_mb_hat, self.updown_filter.to(x.device) * self.subbands, stride=self.subbands)
463
+
464
+ y_g_hat = self.multistream_conv_post(y_mb_hat)
465
+
466
+ return y_g_hat, y_mb_hat
467
+
468
+ def remove_weight_norm(self):
469
+ print('Removing weight norm...')
470
+ for l in self.ups:
471
+ remove_weight_norm(l)
472
+ for l in self.resblocks:
473
+ l.remove_weight_norm()
474
+
475
+
476
+ class DiscriminatorP(torch.nn.Module):
477
+ def __init__(self, period, kernel_size=5, stride=3, use_spectral_norm=False):
478
+ super(DiscriminatorP, self).__init__()
479
+ self.period = period
480
+ self.use_spectral_norm = use_spectral_norm
481
+ norm_f = weight_norm if use_spectral_norm == False else spectral_norm
482
+ self.convs = nn.ModuleList([
483
+ norm_f(Conv2d(1, 32, (kernel_size, 1), (stride, 1), padding=(get_padding(kernel_size, 1), 0))),
484
+ norm_f(Conv2d(32, 128, (kernel_size, 1), (stride, 1), padding=(get_padding(kernel_size, 1), 0))),
485
+ norm_f(Conv2d(128, 512, (kernel_size, 1), (stride, 1), padding=(get_padding(kernel_size, 1), 0))),
486
+ norm_f(Conv2d(512, 1024, (kernel_size, 1), (stride, 1), padding=(get_padding(kernel_size, 1), 0))),
487
+ norm_f(Conv2d(1024, 1024, (kernel_size, 1), 1, padding=(get_padding(kernel_size, 1), 0))),
488
+ ])
489
+ self.conv_post = norm_f(Conv2d(1024, 1, (3, 1), 1, padding=(1, 0)))
490
+
491
+ def forward(self, x):
492
+ fmap = []
493
+
494
+ # 1d to 2d
495
+ b, c, t = x.shape
496
+ if t % self.period != 0: # pad first
497
+ n_pad = self.period - (t % self.period)
498
+ x = F.pad(x, (0, n_pad), "reflect")
499
+ t = t + n_pad
500
+ x = x.view(b, c, t // self.period, self.period)
501
+
502
+ for l in self.convs:
503
+ x = l(x)
504
+ x = F.leaky_relu(x, modules.LRELU_SLOPE)
505
+ fmap.append(x)
506
+ x = self.conv_post(x)
507
+ fmap.append(x)
508
+ x = torch.flatten(x, 1, -1)
509
+
510
+ return x, fmap
511
+
512
+
513
+ class DiscriminatorS(torch.nn.Module):
514
+ def __init__(self, use_spectral_norm=False):
515
+ super(DiscriminatorS, self).__init__()
516
+ norm_f = weight_norm if use_spectral_norm == False else spectral_norm
517
+ self.convs = nn.ModuleList([
518
+ norm_f(Conv1d(1, 16, 15, 1, padding=7)),
519
+ norm_f(Conv1d(16, 64, 41, 4, groups=4, padding=20)),
520
+ norm_f(Conv1d(64, 256, 41, 4, groups=16, padding=20)),
521
+ norm_f(Conv1d(256, 1024, 41, 4, groups=64, padding=20)),
522
+ norm_f(Conv1d(1024, 1024, 41, 4, groups=256, padding=20)),
523
+ norm_f(Conv1d(1024, 1024, 5, 1, padding=2)),
524
+ ])
525
+ self.conv_post = norm_f(Conv1d(1024, 1, 3, 1, padding=1))
526
+
527
+ def forward(self, x):
528
+ fmap = []
529
+
530
+ for l in self.convs:
531
+ x = l(x)
532
+ x = F.leaky_relu(x, modules.LRELU_SLOPE)
533
+ fmap.append(x)
534
+ x = self.conv_post(x)
535
+ fmap.append(x)
536
+ x = torch.flatten(x, 1, -1)
537
+
538
+ return x, fmap
539
+
540
+
541
+ class MultiPeriodDiscriminator(torch.nn.Module):
542
+ def __init__(self, use_spectral_norm=False):
543
+ super(MultiPeriodDiscriminator, self).__init__()
544
+ periods = [2,3,5,7,11]
545
+
546
+ discs = [DiscriminatorS(use_spectral_norm=use_spectral_norm)]
547
+ discs = discs + [DiscriminatorP(i, use_spectral_norm=use_spectral_norm) for i in periods]
548
+ self.discriminators = nn.ModuleList(discs)
549
+
550
+ def forward(self, y, y_hat):
551
+ y_d_rs = []
552
+ y_d_gs = []
553
+ fmap_rs = []
554
+ fmap_gs = []
555
+ for i, d in enumerate(self.discriminators):
556
+ y_d_r, fmap_r = d(y)
557
+ y_d_g, fmap_g = d(y_hat)
558
+ y_d_rs.append(y_d_r)
559
+ y_d_gs.append(y_d_g)
560
+ fmap_rs.append(fmap_r)
561
+ fmap_gs.append(fmap_g)
562
+
563
+ return y_d_rs, y_d_gs, fmap_rs, fmap_gs
564
+
565
+
566
+
567
+ class SynthesizerTrn(nn.Module):
568
+ """
569
+ Synthesizer for Training
570
+ """
571
+
572
+ def __init__(self,
573
+ n_vocab,
574
+ spec_channels,
575
+ segment_size,
576
+ inter_channels,
577
+ hidden_channels,
578
+ filter_channels,
579
+ n_heads,
580
+ n_layers,
581
+ kernel_size,
582
+ p_dropout,
583
+ resblock,
584
+ resblock_kernel_sizes,
585
+ resblock_dilation_sizes,
586
+ upsample_rates,
587
+ upsample_initial_channel,
588
+ upsample_kernel_sizes,
589
+ gen_istft_n_fft,
590
+ gen_istft_hop_size,
591
+ n_speakers=0,
592
+ gin_channels=0,
593
+ use_sdp=False,
594
+ ms_istft_vits=False,
595
+ mb_istft_vits = False,
596
+ subbands = False,
597
+ istft_vits=False,
598
+ **kwargs):
599
+
600
+ super().__init__()
601
+ self.n_vocab = n_vocab
602
+ self.spec_channels = spec_channels
603
+ self.inter_channels = inter_channels
604
+ self.hidden_channels = hidden_channels
605
+ self.filter_channels = filter_channels
606
+ self.n_heads = n_heads
607
+ self.n_layers = n_layers
608
+ self.kernel_size = kernel_size
609
+ self.p_dropout = p_dropout
610
+ self.resblock = resblock
611
+ self.resblock_kernel_sizes = resblock_kernel_sizes
612
+ self.resblock_dilation_sizes = resblock_dilation_sizes
613
+ self.upsample_rates = upsample_rates
614
+ self.upsample_initial_channel = upsample_initial_channel
615
+ self.upsample_kernel_sizes = upsample_kernel_sizes
616
+ self.segment_size = segment_size
617
+ self.n_speakers = n_speakers
618
+ self.gin_channels = gin_channels
619
+ self.ms_istft_vits = ms_istft_vits
620
+ self.mb_istft_vits = mb_istft_vits
621
+ self.istft_vits = istft_vits
622
+
623
+ self.use_sdp = use_sdp
624
+
625
+ self.enc_p = TextEncoder(n_vocab,
626
+ inter_channels,
627
+ hidden_channels,
628
+ filter_channels,
629
+ n_heads,
630
+ n_layers,
631
+ kernel_size,
632
+ p_dropout)
633
+ if mb_istft_vits == True:
634
+ print('Mutli-band iSTFT VITS')
635
+ self.dec = Multiband_iSTFT_Generator(inter_channels, resblock, resblock_kernel_sizes, resblock_dilation_sizes, upsample_rates, upsample_initial_channel, upsample_kernel_sizes, gen_istft_n_fft, gen_istft_hop_size, subbands, gin_channels=gin_channels)
636
+ elif ms_istft_vits == True:
637
+ print('Mutli-stream iSTFT VITS')
638
+ self.dec = Multistream_iSTFT_Generator(inter_channels, resblock, resblock_kernel_sizes, resblock_dilation_sizes, upsample_rates, upsample_initial_channel, upsample_kernel_sizes, gen_istft_n_fft, gen_istft_hop_size, subbands, gin_channels=gin_channels)
639
+ elif istft_vits == True:
640
+ print('iSTFT-VITS')
641
+ self.dec = iSTFT_Generator(inter_channels, resblock, resblock_kernel_sizes, resblock_dilation_sizes, upsample_rates, upsample_initial_channel, upsample_kernel_sizes, gen_istft_n_fft, gen_istft_hop_size, gin_channels=gin_channels)
642
+ else:
643
+ print('Decoder Error in json file')
644
+
645
+ self.enc_q = PosteriorEncoder(spec_channels, inter_channels, hidden_channels, 5, 1, 16, gin_channels=gin_channels)
646
+ self.flow = ResidualCouplingBlock(inter_channels, hidden_channels, 5, 1, 4, gin_channels=gin_channels)
647
+
648
+ if use_sdp:
649
+ self.dp = StochasticDurationPredictor(hidden_channels, 192, 3, 0.5, 4, gin_channels=gin_channels)
650
+ else:
651
+ self.dp = DurationPredictor(hidden_channels, 256, 3, 0.5, gin_channels=gin_channels)
652
+
653
+ if n_speakers > 1:
654
+ self.emb_g = nn.Embedding(n_speakers, gin_channels)
655
+
656
+ def forward(self, x, x_lengths, y, y_lengths, sid=None):
657
+
658
+ x, m_p, logs_p, x_mask = self.enc_p(x, x_lengths)
659
+ if self.n_speakers > 0:
660
+ g = self.emb_g(sid).unsqueeze(-1) # [b, h, 1]
661
+ else:
662
+ g = None
663
+
664
+ z, m_q, logs_q, y_mask = self.enc_q(y, y_lengths, g=g)
665
+ z_p = self.flow(z, y_mask, g=g)
666
+
667
+ with torch.no_grad():
668
+ # negative cross-entropy
669
+ s_p_sq_r = torch.exp(-2 * logs_p) # [b, d, t]
670
+ neg_cent1 = torch.sum(-0.5 * math.log(2 * math.pi) - logs_p, [1], keepdim=True) # [b, 1, t_s]
671
+ neg_cent2 = torch.matmul(-0.5 * (z_p ** 2).transpose(1, 2), s_p_sq_r) # [b, t_t, d] x [b, d, t_s] = [b, t_t, t_s]
672
+ neg_cent3 = torch.matmul(z_p.transpose(1, 2), (m_p * s_p_sq_r)) # [b, t_t, d] x [b, d, t_s] = [b, t_t, t_s]
673
+ neg_cent4 = torch.sum(-0.5 * (m_p ** 2) * s_p_sq_r, [1], keepdim=True) # [b, 1, t_s]
674
+ neg_cent = neg_cent1 + neg_cent2 + neg_cent3 + neg_cent4
675
+
676
+ attn_mask = torch.unsqueeze(x_mask, 2) * torch.unsqueeze(y_mask, -1)
677
+ attn = monotonic_align.maximum_path(neg_cent, attn_mask.squeeze(1)).unsqueeze(1).detach()
678
+
679
+ w = attn.sum(2)
680
+ if self.use_sdp:
681
+ l_length = self.dp(x, x_mask, w, g=g)
682
+ l_length = l_length / torch.sum(x_mask)
683
+ else:
684
+ logw_ = torch.log(w + 1e-6) * x_mask
685
+ logw = self.dp(x, x_mask, g=g)
686
+ l_length = torch.sum((logw - logw_)**2, [1,2]) / torch.sum(x_mask) # for averaging
687
+
688
+ # expand prior
689
+ m_p = torch.matmul(attn.squeeze(1), m_p.transpose(1, 2)).transpose(1, 2)
690
+ logs_p = torch.matmul(attn.squeeze(1), logs_p.transpose(1, 2)).transpose(1, 2)
691
+
692
+ z_slice, ids_slice = commons.rand_slice_segments(z, y_lengths, self.segment_size)
693
+ o, o_mb = self.dec(z_slice, g=g)
694
+ return o, o_mb, l_length, attn, ids_slice, x_mask, y_mask, (z, z_p, m_p, logs_p, m_q, logs_q)
695
+
696
+ def infer(self, x, x_lengths, sid=None, noise_scale=1, length_scale=1, noise_scale_w=1., max_len=None):
697
+ x, m_p, logs_p, x_mask = self.enc_p(x, x_lengths)
698
+ if self.n_speakers > 0:
699
+ g = self.emb_g(sid).unsqueeze(-1) # [b, h, 1]
700
+ else:
701
+ g = None
702
+
703
+ if self.use_sdp:
704
+ logw = self.dp(x, x_mask, g=g, reverse=True, noise_scale=noise_scale_w)
705
+ else:
706
+ logw = self.dp(x, x_mask, g=g)
707
+ w = torch.exp(logw) * x_mask * length_scale
708
+ w_ceil = torch.ceil(w)
709
+ y_lengths = torch.clamp_min(torch.sum(w_ceil, [1, 2]), 1).long()
710
+ y_mask = torch.unsqueeze(commons.sequence_mask(y_lengths, None), 1).to(x_mask.dtype)
711
+ attn_mask = torch.unsqueeze(x_mask, 2) * torch.unsqueeze(y_mask, -1)
712
+ attn = commons.generate_path(w_ceil, attn_mask)
713
+
714
+ m_p = torch.matmul(attn.squeeze(1), m_p.transpose(1, 2)).transpose(1, 2) # [b, t', t], [b, t, d] -> [b, d, t']
715
+ logs_p = torch.matmul(attn.squeeze(1), logs_p.transpose(1, 2)).transpose(1, 2) # [b, t', t], [b, t, d] -> [b, d, t']
716
+
717
+ z_p = m_p + torch.randn_like(m_p) * torch.exp(logs_p) * noise_scale
718
+ z = self.flow(z_p, y_mask, g=g, reverse=True)
719
+ o, o_mb = self.dec((z * y_mask)[:,:,:max_len], g=g)
720
+ return o, o_mb, attn, y_mask, (z, z_p, m_p, logs_p)
721
+
722
+ def voice_conversion(self, y, y_lengths, sid_src, sid_tgt):
723
+ assert self.n_speakers > 0, "n_speakers have to be larger than 0."
724
+ g_src = self.emb_g(sid_src).unsqueeze(-1)
725
+ g_tgt = self.emb_g(sid_tgt).unsqueeze(-1)
726
+ z, m_q, logs_q, y_mask = self.enc_q(y, y_lengths, g=g_src)
727
+ z_p = self.flow(z, y_mask, g=g_src)
728
+ z_hat = self.flow(z_p, y_mask, g=g_tgt, reverse=True)
729
+ o_hat, o_hat_mb = self.dec(z_hat * y_mask, g=g_tgt)
730
+ return o_hat, o_hat_mb, y_mask, (z, z_p, z_hat)
731
+
modules.py ADDED
@@ -0,0 +1,390 @@
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
1
+ import copy
2
+ import math
3
+ import numpy as np
4
+ import scipy
5
+ import torch
6
+ from torch import nn
7
+ from torch.nn import functional as F
8
+
9
+ from torch.nn import Conv1d, ConvTranspose1d, AvgPool1d, Conv2d
10
+ from torch.nn.utils import weight_norm, remove_weight_norm
11
+
12
+ import commons
13
+ from commons import init_weights, get_padding
14
+ from transforms import piecewise_rational_quadratic_transform
15
+
16
+
17
+ LRELU_SLOPE = 0.1
18
+
19
+
20
+ class LayerNorm(nn.Module):
21
+ def __init__(self, channels, eps=1e-5):
22
+ super().__init__()
23
+ self.channels = channels
24
+ self.eps = eps
25
+
26
+ self.gamma = nn.Parameter(torch.ones(channels))
27
+ self.beta = nn.Parameter(torch.zeros(channels))
28
+
29
+ def forward(self, x):
30
+ x = x.transpose(1, -1)
31
+ x = F.layer_norm(x, (self.channels,), self.gamma, self.beta, self.eps)
32
+ return x.transpose(1, -1)
33
+
34
+
35
+ class ConvReluNorm(nn.Module):
36
+ def __init__(self, in_channels, hidden_channels, out_channels, kernel_size, n_layers, p_dropout):
37
+ super().__init__()
38
+ self.in_channels = in_channels
39
+ self.hidden_channels = hidden_channels
40
+ self.out_channels = out_channels
41
+ self.kernel_size = kernel_size
42
+ self.n_layers = n_layers
43
+ self.p_dropout = p_dropout
44
+ assert n_layers > 1, "Number of layers should be larger than 0."
45
+
46
+ self.conv_layers = nn.ModuleList()
47
+ self.norm_layers = nn.ModuleList()
48
+ self.conv_layers.append(nn.Conv1d(in_channels, hidden_channels, kernel_size, padding=kernel_size//2))
49
+ self.norm_layers.append(LayerNorm(hidden_channels))
50
+ self.relu_drop = nn.Sequential(
51
+ nn.ReLU(),
52
+ nn.Dropout(p_dropout))
53
+ for _ in range(n_layers-1):
54
+ self.conv_layers.append(nn.Conv1d(hidden_channels, hidden_channels, kernel_size, padding=kernel_size//2))
55
+ self.norm_layers.append(LayerNorm(hidden_channels))
56
+ self.proj = nn.Conv1d(hidden_channels, out_channels, 1)
57
+ self.proj.weight.data.zero_()
58
+ self.proj.bias.data.zero_()
59
+
60
+ def forward(self, x, x_mask):
61
+ x_org = x
62
+ for i in range(self.n_layers):
63
+ x = self.conv_layers[i](x * x_mask)
64
+ x = self.norm_layers[i](x)
65
+ x = self.relu_drop(x)
66
+ x = x_org + self.proj(x)
67
+ return x * x_mask
68
+
69
+
70
+ class DDSConv(nn.Module):
71
+ """
72
+ Dialted and Depth-Separable Convolution
73
+ """
74
+ def __init__(self, channels, kernel_size, n_layers, p_dropout=0.):
75
+ super().__init__()
76
+ self.channels = channels
77
+ self.kernel_size = kernel_size
78
+ self.n_layers = n_layers
79
+ self.p_dropout = p_dropout
80
+
81
+ self.drop = nn.Dropout(p_dropout)
82
+ self.convs_sep = nn.ModuleList()
83
+ self.convs_1x1 = nn.ModuleList()
84
+ self.norms_1 = nn.ModuleList()
85
+ self.norms_2 = nn.ModuleList()
86
+ for i in range(n_layers):
87
+ dilation = kernel_size ** i
88
+ padding = (kernel_size * dilation - dilation) // 2
89
+ self.convs_sep.append(nn.Conv1d(channels, channels, kernel_size,
90
+ groups=channels, dilation=dilation, padding=padding
91
+ ))
92
+ self.convs_1x1.append(nn.Conv1d(channels, channels, 1))
93
+ self.norms_1.append(LayerNorm(channels))
94
+ self.norms_2.append(LayerNorm(channels))
95
+
96
+ def forward(self, x, x_mask, g=None):
97
+ if g is not None:
98
+ x = x + g
99
+ for i in range(self.n_layers):
100
+ y = self.convs_sep[i](x * x_mask)
101
+ y = self.norms_1[i](y)
102
+ y = F.gelu(y)
103
+ y = self.convs_1x1[i](y)
104
+ y = self.norms_2[i](y)
105
+ y = F.gelu(y)
106
+ y = self.drop(y)
107
+ x = x + y
108
+ return x * x_mask
109
+
110
+
111
+ class WN(torch.nn.Module):
112
+ def __init__(self, hidden_channels, kernel_size, dilation_rate, n_layers, gin_channels=0, p_dropout=0):
113
+ super(WN, self).__init__()
114
+ assert(kernel_size % 2 == 1)
115
+ self.hidden_channels =hidden_channels
116
+ self.kernel_size = kernel_size,
117
+ self.dilation_rate = dilation_rate
118
+ self.n_layers = n_layers
119
+ self.gin_channels = gin_channels
120
+ self.p_dropout = p_dropout
121
+
122
+ self.in_layers = torch.nn.ModuleList()
123
+ self.res_skip_layers = torch.nn.ModuleList()
124
+ self.drop = nn.Dropout(p_dropout)
125
+
126
+ if gin_channels != 0:
127
+ cond_layer = torch.nn.Conv1d(gin_channels, 2*hidden_channels*n_layers, 1)
128
+ self.cond_layer = torch.nn.utils.weight_norm(cond_layer, name='weight')
129
+
130
+ for i in range(n_layers):
131
+ dilation = dilation_rate ** i
132
+ padding = int((kernel_size * dilation - dilation) / 2)
133
+ in_layer = torch.nn.Conv1d(hidden_channels, 2*hidden_channels, kernel_size,
134
+ dilation=dilation, padding=padding)
135
+ in_layer = torch.nn.utils.weight_norm(in_layer, name='weight')
136
+ self.in_layers.append(in_layer)
137
+
138
+ # last one is not necessary
139
+ if i < n_layers - 1:
140
+ res_skip_channels = 2 * hidden_channels
141
+ else:
142
+ res_skip_channels = hidden_channels
143
+
144
+ res_skip_layer = torch.nn.Conv1d(hidden_channels, res_skip_channels, 1)
145
+ res_skip_layer = torch.nn.utils.weight_norm(res_skip_layer, name='weight')
146
+ self.res_skip_layers.append(res_skip_layer)
147
+
148
+ def forward(self, x, x_mask, g=None, **kwargs):
149
+ output = torch.zeros_like(x)
150
+ n_channels_tensor = torch.IntTensor([self.hidden_channels])
151
+
152
+ if g is not None:
153
+ g = self.cond_layer(g)
154
+
155
+ for i in range(self.n_layers):
156
+ x_in = self.in_layers[i](x)
157
+ if g is not None:
158
+ cond_offset = i * 2 * self.hidden_channels
159
+ g_l = g[:,cond_offset:cond_offset+2*self.hidden_channels,:]
160
+ else:
161
+ g_l = torch.zeros_like(x_in)
162
+
163
+ acts = commons.fused_add_tanh_sigmoid_multiply(
164
+ x_in,
165
+ g_l,
166
+ n_channels_tensor)
167
+ acts = self.drop(acts)
168
+
169
+ res_skip_acts = self.res_skip_layers[i](acts)
170
+ if i < self.n_layers - 1:
171
+ res_acts = res_skip_acts[:,:self.hidden_channels,:]
172
+ x = (x + res_acts) * x_mask
173
+ output = output + res_skip_acts[:,self.hidden_channels:,:]
174
+ else:
175
+ output = output + res_skip_acts
176
+ return output * x_mask
177
+
178
+ def remove_weight_norm(self):
179
+ if self.gin_channels != 0:
180
+ torch.nn.utils.remove_weight_norm(self.cond_layer)
181
+ for l in self.in_layers:
182
+ torch.nn.utils.remove_weight_norm(l)
183
+ for l in self.res_skip_layers:
184
+ torch.nn.utils.remove_weight_norm(l)
185
+
186
+
187
+ class ResBlock1(torch.nn.Module):
188
+ def __init__(self, channels, kernel_size=3, dilation=(1, 3, 5)):
189
+ super(ResBlock1, self).__init__()
190
+ self.convs1 = nn.ModuleList([
191
+ weight_norm(Conv1d(channels, channels, kernel_size, 1, dilation=dilation[0],
192
+ padding=get_padding(kernel_size, dilation[0]))),
193
+ weight_norm(Conv1d(channels, channels, kernel_size, 1, dilation=dilation[1],
194
+ padding=get_padding(kernel_size, dilation[1]))),
195
+ weight_norm(Conv1d(channels, channels, kernel_size, 1, dilation=dilation[2],
196
+ padding=get_padding(kernel_size, dilation[2])))
197
+ ])
198
+ self.convs1.apply(init_weights)
199
+
200
+ self.convs2 = nn.ModuleList([
201
+ weight_norm(Conv1d(channels, channels, kernel_size, 1, dilation=1,
202
+ padding=get_padding(kernel_size, 1))),
203
+ weight_norm(Conv1d(channels, channels, kernel_size, 1, dilation=1,
204
+ padding=get_padding(kernel_size, 1))),
205
+ weight_norm(Conv1d(channels, channels, kernel_size, 1, dilation=1,
206
+ padding=get_padding(kernel_size, 1)))
207
+ ])
208
+ self.convs2.apply(init_weights)
209
+
210
+ def forward(self, x, x_mask=None):
211
+ for c1, c2 in zip(self.convs1, self.convs2):
212
+ xt = F.leaky_relu(x, LRELU_SLOPE)
213
+ if x_mask is not None:
214
+ xt = xt * x_mask
215
+ xt = c1(xt)
216
+ xt = F.leaky_relu(xt, LRELU_SLOPE)
217
+ if x_mask is not None:
218
+ xt = xt * x_mask
219
+ xt = c2(xt)
220
+ x = xt + x
221
+ if x_mask is not None:
222
+ x = x * x_mask
223
+ return x
224
+
225
+ def remove_weight_norm(self):
226
+ for l in self.convs1:
227
+ remove_weight_norm(l)
228
+ for l in self.convs2:
229
+ remove_weight_norm(l)
230
+
231
+
232
+ class ResBlock2(torch.nn.Module):
233
+ def __init__(self, channels, kernel_size=3, dilation=(1, 3)):
234
+ super(ResBlock2, self).__init__()
235
+ self.convs = nn.ModuleList([
236
+ weight_norm(Conv1d(channels, channels, kernel_size, 1, dilation=dilation[0],
237
+ padding=get_padding(kernel_size, dilation[0]))),
238
+ weight_norm(Conv1d(channels, channels, kernel_size, 1, dilation=dilation[1],
239
+ padding=get_padding(kernel_size, dilation[1])))
240
+ ])
241
+ self.convs.apply(init_weights)
242
+
243
+ def forward(self, x, x_mask=None):
244
+ for c in self.convs:
245
+ xt = F.leaky_relu(x, LRELU_SLOPE)
246
+ if x_mask is not None:
247
+ xt = xt * x_mask
248
+ xt = c(xt)
249
+ x = xt + x
250
+ if x_mask is not None:
251
+ x = x * x_mask
252
+ return x
253
+
254
+ def remove_weight_norm(self):
255
+ for l in self.convs:
256
+ remove_weight_norm(l)
257
+
258
+
259
+ class Log(nn.Module):
260
+ def forward(self, x, x_mask, reverse=False, **kwargs):
261
+ if not reverse:
262
+ y = torch.log(torch.clamp_min(x, 1e-5)) * x_mask
263
+ logdet = torch.sum(-y, [1, 2])
264
+ return y, logdet
265
+ else:
266
+ x = torch.exp(x) * x_mask
267
+ return x
268
+
269
+
270
+ class Flip(nn.Module):
271
+ def forward(self, x, *args, reverse=False, **kwargs):
272
+ x = torch.flip(x, [1])
273
+ if not reverse:
274
+ logdet = torch.zeros(x.size(0)).to(dtype=x.dtype, device=x.device)
275
+ return x, logdet
276
+ else:
277
+ return x
278
+
279
+
280
+ class ElementwiseAffine(nn.Module):
281
+ def __init__(self, channels):
282
+ super().__init__()
283
+ self.channels = channels
284
+ self.m = nn.Parameter(torch.zeros(channels,1))
285
+ self.logs = nn.Parameter(torch.zeros(channels,1))
286
+
287
+ def forward(self, x, x_mask, reverse=False, **kwargs):
288
+ if not reverse:
289
+ y = self.m + torch.exp(self.logs) * x
290
+ y = y * x_mask
291
+ logdet = torch.sum(self.logs * x_mask, [1,2])
292
+ return y, logdet
293
+ else:
294
+ x = (x - self.m) * torch.exp(-self.logs) * x_mask
295
+ return x
296
+
297
+
298
+ class ResidualCouplingLayer(nn.Module):
299
+ def __init__(self,
300
+ channels,
301
+ hidden_channels,
302
+ kernel_size,
303
+ dilation_rate,
304
+ n_layers,
305
+ p_dropout=0,
306
+ gin_channels=0,
307
+ mean_only=False):
308
+ assert channels % 2 == 0, "channels should be divisible by 2"
309
+ super().__init__()
310
+ self.channels = channels
311
+ self.hidden_channels = hidden_channels
312
+ self.kernel_size = kernel_size
313
+ self.dilation_rate = dilation_rate
314
+ self.n_layers = n_layers
315
+ self.half_channels = channels // 2
316
+ self.mean_only = mean_only
317
+
318
+ self.pre = nn.Conv1d(self.half_channels, hidden_channels, 1)
319
+ self.enc = WN(hidden_channels, kernel_size, dilation_rate, n_layers, p_dropout=p_dropout, gin_channels=gin_channels)
320
+ self.post = nn.Conv1d(hidden_channels, self.half_channels * (2 - mean_only), 1)
321
+ self.post.weight.data.zero_()
322
+ self.post.bias.data.zero_()
323
+
324
+ def forward(self, x, x_mask, g=None, reverse=False):
325
+ x0, x1 = torch.split(x, [self.half_channels]*2, 1)
326
+ h = self.pre(x0) * x_mask
327
+ h = self.enc(h, x_mask, g=g)
328
+ stats = self.post(h) * x_mask
329
+ if not self.mean_only:
330
+ m, logs = torch.split(stats, [self.half_channels]*2, 1)
331
+ else:
332
+ m = stats
333
+ logs = torch.zeros_like(m)
334
+
335
+ if not reverse:
336
+ x1 = m + x1 * torch.exp(logs) * x_mask
337
+ x = torch.cat([x0, x1], 1)
338
+ logdet = torch.sum(logs, [1,2])
339
+ return x, logdet
340
+ else:
341
+ x1 = (x1 - m) * torch.exp(-logs) * x_mask
342
+ x = torch.cat([x0, x1], 1)
343
+ return x
344
+
345
+
346
+ class ConvFlow(nn.Module):
347
+ def __init__(self, in_channels, filter_channels, kernel_size, n_layers, num_bins=10, tail_bound=5.0):
348
+ super().__init__()
349
+ self.in_channels = in_channels
350
+ self.filter_channels = filter_channels
351
+ self.kernel_size = kernel_size
352
+ self.n_layers = n_layers
353
+ self.num_bins = num_bins
354
+ self.tail_bound = tail_bound
355
+ self.half_channels = in_channels // 2
356
+
357
+ self.pre = nn.Conv1d(self.half_channels, filter_channels, 1)
358
+ self.convs = DDSConv(filter_channels, kernel_size, n_layers, p_dropout=0.)
359
+ self.proj = nn.Conv1d(filter_channels, self.half_channels * (num_bins * 3 - 1), 1)
360
+ self.proj.weight.data.zero_()
361
+ self.proj.bias.data.zero_()
362
+
363
+ def forward(self, x, x_mask, g=None, reverse=False):
364
+ x0, x1 = torch.split(x, [self.half_channels]*2, 1)
365
+ h = self.pre(x0)
366
+ h = self.convs(h, x_mask, g=g)
367
+ h = self.proj(h) * x_mask
368
+
369
+ b, c, t = x0.shape
370
+ h = h.reshape(b, c, -1, t).permute(0, 1, 3, 2) # [b, cx?, t] -> [b, c, t, ?]
371
+
372
+ unnormalized_widths = h[..., :self.num_bins] / math.sqrt(self.filter_channels)
373
+ unnormalized_heights = h[..., self.num_bins:2*self.num_bins] / math.sqrt(self.filter_channels)
374
+ unnormalized_derivatives = h[..., 2 * self.num_bins:]
375
+
376
+ x1, logabsdet = piecewise_rational_quadratic_transform(x1,
377
+ unnormalized_widths,
378
+ unnormalized_heights,
379
+ unnormalized_derivatives,
380
+ inverse=reverse,
381
+ tails='linear',
382
+ tail_bound=self.tail_bound
383
+ )
384
+
385
+ x = torch.cat([x0, x1], 1) * x_mask
386
+ logdet = torch.sum(logabsdet * x_mask, [1,2])
387
+ if not reverse:
388
+ return x, logdet
389
+ else:
390
+ return x
pqmf.py ADDED
@@ -0,0 +1,116 @@
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
1
+ # -*- coding: utf-8 -*-
2
+
3
+ # Copyright 2020 Tomoki Hayashi
4
+ # MIT License (https://opensource.org/licenses/MIT)
5
+
6
+ """Pseudo QMF modules."""
7
+
8
+ import numpy as np
9
+ import torch
10
+ import torch.nn.functional as F
11
+
12
+ from scipy.signal import kaiser
13
+
14
+
15
+ def design_prototype_filter(taps=62, cutoff_ratio=0.15, beta=9.0):
16
+ """Design prototype filter for PQMF.
17
+ This method is based on `A Kaiser window approach for the design of prototype
18
+ filters of cosine modulated filterbanks`_.
19
+ Args:
20
+ taps (int): The number of filter taps.
21
+ cutoff_ratio (float): Cut-off frequency ratio.
22
+ beta (float): Beta coefficient for kaiser window.
23
+ Returns:
24
+ ndarray: Impluse response of prototype filter (taps + 1,).
25
+ .. _`A Kaiser window approach for the design of prototype filters of cosine modulated filterbanks`:
26
+ https://ieeexplore.ieee.org/abstract/document/681427
27
+ """
28
+ # check the arguments are valid
29
+ assert taps % 2 == 0, "The number of taps mush be even number."
30
+ assert 0.0 < cutoff_ratio < 1.0, "Cutoff ratio must be > 0.0 and < 1.0."
31
+
32
+ # make initial filter
33
+ omega_c = np.pi * cutoff_ratio
34
+ with np.errstate(invalid='ignore'):
35
+ h_i = np.sin(omega_c * (np.arange(taps + 1) - 0.5 * taps)) \
36
+ / (np.pi * (np.arange(taps + 1) - 0.5 * taps))
37
+ h_i[taps // 2] = np.cos(0) * cutoff_ratio # fix nan due to indeterminate form
38
+
39
+ # apply kaiser window
40
+ w = kaiser(taps + 1, beta)
41
+ h = h_i * w
42
+
43
+ return h
44
+
45
+
46
+ class PQMF(torch.nn.Module):
47
+ """PQMF module.
48
+ This module is based on `Near-perfect-reconstruction pseudo-QMF banks`_.
49
+ .. _`Near-perfect-reconstruction pseudo-QMF banks`:
50
+ https://ieeexplore.ieee.org/document/258122
51
+ """
52
+
53
+ def __init__(self, device, subbands=4, taps=62, cutoff_ratio=0.15, beta=9.0):
54
+ """Initilize PQMF module.
55
+ Args:
56
+ subbands (int): The number of subbands.
57
+ taps (int): The number of filter taps.
58
+ cutoff_ratio (float): Cut-off frequency ratio.
59
+ beta (float): Beta coefficient for kaiser window.
60
+ """
61
+ super(PQMF, self).__init__()
62
+
63
+ # define filter coefficient
64
+ h_proto = design_prototype_filter(taps, cutoff_ratio, beta)
65
+ h_analysis = np.zeros((subbands, len(h_proto)))
66
+ h_synthesis = np.zeros((subbands, len(h_proto)))
67
+ for k in range(subbands):
68
+ h_analysis[k] = 2 * h_proto * np.cos(
69
+ (2 * k + 1) * (np.pi / (2 * subbands)) *
70
+ (np.arange(taps + 1) - ((taps - 1) / 2)) +
71
+ (-1) ** k * np.pi / 4)
72
+ h_synthesis[k] = 2 * h_proto * np.cos(
73
+ (2 * k + 1) * (np.pi / (2 * subbands)) *
74
+ (np.arange(taps + 1) - ((taps - 1) / 2)) -
75
+ (-1) ** k * np.pi / 4)
76
+
77
+ # convert to tensor
78
+ analysis_filter = torch.from_numpy(h_analysis).float().unsqueeze(1).to(device)#.cuda(device)
79
+ synthesis_filter = torch.from_numpy(h_synthesis).float().unsqueeze(0).to(device)#.cuda(device)
80
+
81
+ # register coefficients as beffer
82
+ self.register_buffer("analysis_filter", analysis_filter)
83
+ self.register_buffer("synthesis_filter", synthesis_filter)
84
+
85
+ # filter for downsampling & upsampling
86
+ updown_filter = torch.zeros((subbands, subbands, subbands)).float().to(device)#.cuda(device)
87
+ for k in range(subbands):
88
+ updown_filter[k, k, 0] = 1.0
89
+ self.register_buffer("updown_filter", updown_filter)
90
+ self.subbands = subbands
91
+
92
+ # keep padding info
93
+ self.pad_fn = torch.nn.ConstantPad1d(taps // 2, 0.0)
94
+
95
+ def analysis(self, x):
96
+ """Analysis with PQMF.
97
+ Args:
98
+ x (Tensor): Input tensor (B, 1, T).
99
+ Returns:
100
+ Tensor: Output tensor (B, subbands, T // subbands).
101
+ """
102
+ x = F.conv1d(self.pad_fn(x), self.analysis_filter)
103
+ return F.conv1d(x, self.updown_filter, stride=self.subbands)
104
+
105
+ def synthesis(self, x):
106
+ """Synthesis with PQMF.
107
+ Args:
108
+ x (Tensor): Input tensor (B, subbands, T // subbands).
109
+ Returns:
110
+ Tensor: Output tensor (B, 1, T).
111
+ """
112
+ # NOTE(kan-bayashi): Power will be dreased so here multipy by # subbands.
113
+ # Not sure this is the correct way, it is better to check again.
114
+ # TODO(kan-bayashi): Understand the reconstruction procedure
115
+ x = F.conv_transpose1d(x, self.updown_filter * self.subbands, stride=self.subbands)
116
+ return F.conv1d(self.pad_fn(x), self.synthesis_filter)
requirements.txt ADDED
@@ -0,0 +1,14 @@
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
1
+ Cython==0.29.21
2
+ librosa==0.8.0
3
+ matplotlib==3.3.1
4
+ numpy==1.18.5
5
+ phonemizer==2.2.1
6
+ scipy==1.5.2
7
+ tensorboard==2.3.0
8
+ torch==1.6.0
9
+ torchvision==0.7.0
10
+ Unidecode==1.1.1
11
+ pysoundfile==0.9.0.post1
12
+ pyopenjtalk==0.2.0
13
+ monotonic-align
14
+ gradio
stft.py ADDED
@@ -0,0 +1,209 @@
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
1
+ """
2
+ BSD 3-Clause License
3
+ Copyright (c) 2017, Prem Seetharaman
4
+ All rights reserved.
5
+ * Redistribution and use in source and binary forms, with or without
6
+ modification, are permitted provided that the following conditions are met:
7
+ * Redistributions of source code must retain the above copyright notice,
8
+ this list of conditions and the following disclaimer.
9
+ * Redistributions in binary form must reproduce the above copyright notice, this
10
+ list of conditions and the following disclaimer in the
11
+ documentation and/or other materials provided with the distribution.
12
+ * Neither the name of the copyright holder nor the names of its
13
+ contributors may be used to endorse or promote products derived from this
14
+ software without specific prior written permission.
15
+ THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" AND
16
+ ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED
17
+ WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE
18
+ DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT HOLDER OR CONTRIBUTORS BE LIABLE FOR
19
+ ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES
20
+ (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES;
21
+ LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON
22
+ ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT
23
+ (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS
24
+ SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
25
+ """
26
+
27
+ import torch
28
+ import numpy as np
29
+ import torch.nn.functional as F
30
+ from torch.autograd import Variable
31
+ from scipy.signal import get_window
32
+ from librosa.util import pad_center, tiny
33
+ import librosa.util as librosa_util
34
+
35
+ def window_sumsquare(window, n_frames, hop_length=200, win_length=800,
36
+ n_fft=800, dtype=np.float32, norm=None):
37
+ """
38
+ # from librosa 0.6
39
+ Compute the sum-square envelope of a window function at a given hop length.
40
+ This is used to estimate modulation effects induced by windowing
41
+ observations in short-time fourier transforms.
42
+ Parameters
43
+ ----------
44
+ window : string, tuple, number, callable, or list-like
45
+ Window specification, as in `get_window`
46
+ n_frames : int > 0
47
+ The number of analysis frames
48
+ hop_length : int > 0
49
+ The number of samples to advance between frames
50
+ win_length : [optional]
51
+ The length of the window function. By default, this matches `n_fft`.
52
+ n_fft : int > 0
53
+ The length of each analysis frame.
54
+ dtype : np.dtype
55
+ The data type of the output
56
+ Returns
57
+ -------
58
+ wss : np.ndarray, shape=`(n_fft + hop_length * (n_frames - 1))`
59
+ The sum-squared envelope of the window function
60
+ """
61
+ if win_length is None:
62
+ win_length = n_fft
63
+
64
+ n = n_fft + hop_length * (n_frames - 1)
65
+ x = np.zeros(n, dtype=dtype)
66
+
67
+ # Compute the squared window at the desired length
68
+ win_sq = get_window(window, win_length, fftbins=True)
69
+ win_sq = librosa_util.normalize(win_sq, norm=norm)**2
70
+ win_sq = librosa_util.pad_center(win_sq, n_fft)
71
+
72
+ # Fill the envelope
73
+ for i in range(n_frames):
74
+ sample = i * hop_length
75
+ x[sample:min(n, sample + n_fft)] += win_sq[:max(0, min(n_fft, n - sample))]
76
+ return x
77
+
78
+
79
+ class STFT(torch.nn.Module):
80
+ """adapted from Prem Seetharaman's https://github.com/pseeth/pytorch-stft"""
81
+ def __init__(self, filter_length=800, hop_length=200, win_length=800,
82
+ window='hann'):
83
+ super(STFT, self).__init__()
84
+ self.filter_length = filter_length
85
+ self.hop_length = hop_length
86
+ self.win_length = win_length
87
+ self.window = window
88
+ self.forward_transform = None
89
+ scale = self.filter_length / self.hop_length
90
+ fourier_basis = np.fft.fft(np.eye(self.filter_length))
91
+
92
+ cutoff = int((self.filter_length / 2 + 1))
93
+ fourier_basis = np.vstack([np.real(fourier_basis[:cutoff, :]),
94
+ np.imag(fourier_basis[:cutoff, :])])
95
+
96
+ forward_basis = torch.FloatTensor(fourier_basis[:, None, :])
97
+ inverse_basis = torch.FloatTensor(
98
+ np.linalg.pinv(scale * fourier_basis).T[:, None, :])
99
+
100
+ if window is not None:
101
+ assert(filter_length >= win_length)
102
+ # get window and zero center pad it to filter_length
103
+ fft_window = get_window(window, win_length, fftbins=True)
104
+ fft_window = pad_center(fft_window, filter_length)
105
+ fft_window = torch.from_numpy(fft_window).float()
106
+
107
+ # window the bases
108
+ forward_basis *= fft_window
109
+ inverse_basis *= fft_window
110
+
111
+ self.register_buffer('forward_basis', forward_basis.float())
112
+ self.register_buffer('inverse_basis', inverse_basis.float())
113
+
114
+ def transform(self, input_data):
115
+ num_batches = input_data.size(0)
116
+ num_samples = input_data.size(1)
117
+
118
+ self.num_samples = num_samples
119
+
120
+ # similar to librosa, reflect-pad the input
121
+ input_data = input_data.view(num_batches, 1, num_samples)
122
+ input_data = F.pad(
123
+ input_data.unsqueeze(1),
124
+ (int(self.filter_length / 2), int(self.filter_length / 2), 0, 0),
125
+ mode='reflect')
126
+ input_data = input_data.squeeze(1)
127
+
128
+ forward_transform = F.conv1d(
129
+ input_data,
130
+ Variable(self.forward_basis, requires_grad=False),
131
+ stride=self.hop_length,
132
+ padding=0)
133
+
134
+ cutoff = int((self.filter_length / 2) + 1)
135
+ real_part = forward_transform[:, :cutoff, :]
136
+ imag_part = forward_transform[:, cutoff:, :]
137
+
138
+ magnitude = torch.sqrt(real_part**2 + imag_part**2)
139
+ phase = torch.autograd.Variable(
140
+ torch.atan2(imag_part.data, real_part.data))
141
+
142
+ return magnitude, phase
143
+
144
+ def inverse(self, magnitude, phase):
145
+ recombine_magnitude_phase = torch.cat(
146
+ [magnitude*torch.cos(phase), magnitude*torch.sin(phase)], dim=1)
147
+
148
+ inverse_transform = F.conv_transpose1d(
149
+ recombine_magnitude_phase,
150
+ Variable(self.inverse_basis, requires_grad=False),
151
+ stride=self.hop_length,
152
+ padding=0)
153
+
154
+ if self.window is not None:
155
+ window_sum = window_sumsquare(
156
+ self.window, magnitude.size(-1), hop_length=self.hop_length,
157
+ win_length=self.win_length, n_fft=self.filter_length,
158
+ dtype=np.float32)
159
+ # remove modulation effects
160
+ approx_nonzero_indices = torch.from_numpy(
161
+ np.where(window_sum > tiny(window_sum))[0])
162
+ window_sum = torch.autograd.Variable(
163
+ torch.from_numpy(window_sum), requires_grad=False)
164
+ window_sum = window_sum.to(inverse_transform.device()) if magnitude.is_cuda else window_sum
165
+ inverse_transform[:, :, approx_nonzero_indices] /= window_sum[approx_nonzero_indices]
166
+
167
+ # scale by hop ratio
168
+ inverse_transform *= float(self.filter_length) / self.hop_length
169
+
170
+ inverse_transform = inverse_transform[:, :, int(self.filter_length/2):]
171
+ inverse_transform = inverse_transform[:, :, :-int(self.filter_length/2):]
172
+
173
+ return inverse_transform
174
+
175
+ def forward(self, input_data):
176
+ self.magnitude, self.phase = self.transform(input_data)
177
+ reconstruction = self.inverse(self.magnitude, self.phase)
178
+ return reconstruction
179
+
180
+
181
+ class TorchSTFT(torch.nn.Module):
182
+ def __init__(self, filter_length=800, hop_length=200, win_length=800, window='hann'):
183
+ super().__init__()
184
+ self.filter_length = filter_length
185
+ self.hop_length = hop_length
186
+ self.win_length = win_length
187
+ self.window = torch.from_numpy(get_window(window, win_length, fftbins=True).astype(np.float32))
188
+
189
+ def transform(self, input_data):
190
+ forward_transform = torch.stft(
191
+ input_data,
192
+ self.filter_length, self.hop_length, self.win_length, window=self.window,
193
+ return_complex=True)
194
+
195
+ return torch.abs(forward_transform), torch.angle(forward_transform)
196
+
197
+ def inverse(self, magnitude, phase):
198
+ inverse_transform = torch.istft(
199
+ magnitude * torch.exp(phase * 1j),
200
+ self.filter_length, self.hop_length, self.win_length, window=self.window.to(magnitude.device))
201
+
202
+ return inverse_transform.unsqueeze(-2) # unsqueeze to stay consistent with conv_transpose1d implementation
203
+
204
+ def forward(self, input_data):
205
+ self.magnitude, self.phase = self.transform(input_data)
206
+ reconstruction = self.inverse(self.magnitude, self.phase)
207
+ return reconstruction
208
+
209
+
stft_loss.py ADDED
@@ -0,0 +1,136 @@
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
1
+ # -*- coding: utf-8 -*-
2
+
3
+ # Copyright 2019 Tomoki Hayashi
4
+ # MIT License (https://opensource.org/licenses/MIT)
5
+
6
+ """STFT-based Loss modules."""
7
+
8
+ import torch
9
+ import torch.nn.functional as F
10
+
11
+
12
+ def stft(x, fft_size, hop_size, win_length, window):
13
+ """Perform STFT and convert to magnitude spectrogram.
14
+ Args:
15
+ x (Tensor): Input signal tensor (B, T).
16
+ fft_size (int): FFT size.
17
+ hop_size (int): Hop size.
18
+ win_length (int): Window length.
19
+ window (str): Window function type.
20
+ Returns:
21
+ Tensor: Magnitude spectrogram (B, #frames, fft_size // 2 + 1).
22
+ """
23
+ x_stft = torch.stft(x, fft_size, hop_size, win_length, window.to(x.device))
24
+ real = x_stft[..., 0]
25
+ imag = x_stft[..., 1]
26
+
27
+ # NOTE(kan-bayashi): clamp is needed to avoid nan or inf
28
+ return torch.sqrt(torch.clamp(real ** 2 + imag ** 2, min=1e-7)).transpose(2, 1)
29
+
30
+
31
+ class SpectralConvergengeLoss(torch.nn.Module):
32
+ """Spectral convergence loss module."""
33
+
34
+ def __init__(self):
35
+ """Initilize spectral convergence loss module."""
36
+ super(SpectralConvergengeLoss, self).__init__()
37
+
38
+ def forward(self, x_mag, y_mag):
39
+ """Calculate forward propagation.
40
+ Args:
41
+ x_mag (Tensor): Magnitude spectrogram of predicted signal (B, #frames, #freq_bins).
42
+ y_mag (Tensor): Magnitude spectrogram of groundtruth signal (B, #frames, #freq_bins).
43
+ Returns:
44
+ Tensor: Spectral convergence loss value.
45
+ """
46
+ return torch.norm(y_mag - x_mag, p="fro") / torch.norm(y_mag, p="fro")
47
+
48
+
49
+ class LogSTFTMagnitudeLoss(torch.nn.Module):
50
+ """Log STFT magnitude loss module."""
51
+
52
+ def __init__(self):
53
+ """Initilize los STFT magnitude loss module."""
54
+ super(LogSTFTMagnitudeLoss, self).__init__()
55
+
56
+ def forward(self, x_mag, y_mag):
57
+ """Calculate forward propagation.
58
+ Args:
59
+ x_mag (Tensor): Magnitude spectrogram of predicted signal (B, #frames, #freq_bins).
60
+ y_mag (Tensor): Magnitude spectrogram of groundtruth signal (B, #frames, #freq_bins).
61
+ Returns:
62
+ Tensor: Log STFT magnitude loss value.
63
+ """
64
+ return F.l1_loss(torch.log(y_mag), torch.log(x_mag))
65
+
66
+
67
+ class STFTLoss(torch.nn.Module):
68
+ """STFT loss module."""
69
+
70
+ def __init__(self, fft_size=1024, shift_size=120, win_length=600, window="hann_window"):
71
+ """Initialize STFT loss module."""
72
+ super(STFTLoss, self).__init__()
73
+ self.fft_size = fft_size
74
+ self.shift_size = shift_size
75
+ self.win_length = win_length
76
+ self.window = getattr(torch, window)(win_length)
77
+ self.spectral_convergenge_loss = SpectralConvergengeLoss()
78
+ self.log_stft_magnitude_loss = LogSTFTMagnitudeLoss()
79
+
80
+ def forward(self, x, y):
81
+ """Calculate forward propagation.
82
+ Args:
83
+ x (Tensor): Predicted signal (B, T).
84
+ y (Tensor): Groundtruth signal (B, T).
85
+ Returns:
86
+ Tensor: Spectral convergence loss value.
87
+ Tensor: Log STFT magnitude loss value.
88
+ """
89
+ x_mag = stft(x, self.fft_size, self.shift_size, self.win_length, self.window)
90
+ y_mag = stft(y, self.fft_size, self.shift_size, self.win_length, self.window)
91
+ sc_loss = self.spectral_convergenge_loss(x_mag, y_mag)
92
+ mag_loss = self.log_stft_magnitude_loss(x_mag, y_mag)
93
+
94
+ return sc_loss, mag_loss
95
+
96
+
97
+ class MultiResolutionSTFTLoss(torch.nn.Module):
98
+ """Multi resolution STFT loss module."""
99
+
100
+ def __init__(self,
101
+ fft_sizes=[1024, 2048, 512],
102
+ hop_sizes=[120, 240, 50],
103
+ win_lengths=[600, 1200, 240],
104
+ window="hann_window"):
105
+ """Initialize Multi resolution STFT loss module.
106
+ Args:
107
+ fft_sizes (list): List of FFT sizes.
108
+ hop_sizes (list): List of hop sizes.
109
+ win_lengths (list): List of window lengths.
110
+ window (str): Window function type.
111
+ """
112
+ super(MultiResolutionSTFTLoss, self).__init__()
113
+ assert len(fft_sizes) == len(hop_sizes) == len(win_lengths)
114
+ self.stft_losses = torch.nn.ModuleList()
115
+ for fs, ss, wl in zip(fft_sizes, hop_sizes, win_lengths):
116
+ self.stft_losses += [STFTLoss(fs, ss, wl, window)]
117
+
118
+ def forward(self, x, y):
119
+ """Calculate forward propagation.
120
+ Args:
121
+ x (Tensor): Predicted signal (B, T).
122
+ y (Tensor): Groundtruth signal (B, T).
123
+ Returns:
124
+ Tensor: Multi resolution spectral convergence loss value.
125
+ Tensor: Multi resolution log STFT magnitude loss value.
126
+ """
127
+ sc_loss = 0.0
128
+ mag_loss = 0.0
129
+ for f in self.stft_losses:
130
+ sc_l, mag_l = f(x, y)
131
+ sc_loss += sc_l
132
+ mag_loss += mag_l
133
+ sc_loss /= len(self.stft_losses)
134
+ mag_loss /= len(self.stft_losses)
135
+
136
+ return sc_loss, mag_loss
transforms.py ADDED
@@ -0,0 +1,193 @@
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
1
+ import torch
2
+ from torch.nn import functional as F
3
+
4
+ import numpy as np
5
+
6
+
7
+ DEFAULT_MIN_BIN_WIDTH = 1e-3
8
+ DEFAULT_MIN_BIN_HEIGHT = 1e-3
9
+ DEFAULT_MIN_DERIVATIVE = 1e-3
10
+
11
+
12
+ def piecewise_rational_quadratic_transform(inputs,
13
+ unnormalized_widths,
14
+ unnormalized_heights,
15
+ unnormalized_derivatives,
16
+ inverse=False,
17
+ tails=None,
18
+ tail_bound=1.,
19
+ min_bin_width=DEFAULT_MIN_BIN_WIDTH,
20
+ min_bin_height=DEFAULT_MIN_BIN_HEIGHT,
21
+ min_derivative=DEFAULT_MIN_DERIVATIVE):
22
+
23
+ if tails is None:
24
+ spline_fn = rational_quadratic_spline
25
+ spline_kwargs = {}
26
+ else:
27
+ spline_fn = unconstrained_rational_quadratic_spline
28
+ spline_kwargs = {
29
+ 'tails': tails,
30
+ 'tail_bound': tail_bound
31
+ }
32
+
33
+ outputs, logabsdet = spline_fn(
34
+ inputs=inputs,
35
+ unnormalized_widths=unnormalized_widths,
36
+ unnormalized_heights=unnormalized_heights,
37
+ unnormalized_derivatives=unnormalized_derivatives,
38
+ inverse=inverse,
39
+ min_bin_width=min_bin_width,
40
+ min_bin_height=min_bin_height,
41
+ min_derivative=min_derivative,
42
+ **spline_kwargs
43
+ )
44
+ return outputs, logabsdet
45
+
46
+
47
+ def searchsorted(bin_locations, inputs, eps=1e-6):
48
+ bin_locations[..., -1] += eps
49
+ return torch.sum(
50
+ inputs[..., None] >= bin_locations,
51
+ dim=-1
52
+ ) - 1
53
+
54
+
55
+ def unconstrained_rational_quadratic_spline(inputs,
56
+ unnormalized_widths,
57
+ unnormalized_heights,
58
+ unnormalized_derivatives,
59
+ inverse=False,
60
+ tails='linear',
61
+ tail_bound=1.,
62
+ min_bin_width=DEFAULT_MIN_BIN_WIDTH,
63
+ min_bin_height=DEFAULT_MIN_BIN_HEIGHT,
64
+ min_derivative=DEFAULT_MIN_DERIVATIVE):
65
+ inside_interval_mask = (inputs >= -tail_bound) & (inputs <= tail_bound)
66
+ outside_interval_mask = ~inside_interval_mask
67
+
68
+ outputs = torch.zeros_like(inputs)
69
+ logabsdet = torch.zeros_like(inputs)
70
+
71
+ if tails == 'linear':
72
+ unnormalized_derivatives = F.pad(unnormalized_derivatives, pad=(1, 1))
73
+ constant = np.log(np.exp(1 - min_derivative) - 1)
74
+ unnormalized_derivatives[..., 0] = constant
75
+ unnormalized_derivatives[..., -1] = constant
76
+
77
+ outputs[outside_interval_mask] = inputs[outside_interval_mask]
78
+ logabsdet[outside_interval_mask] = 0
79
+ else:
80
+ raise RuntimeError('{} tails are not implemented.'.format(tails))
81
+
82
+ outputs[inside_interval_mask], logabsdet[inside_interval_mask] = rational_quadratic_spline(
83
+ inputs=inputs[inside_interval_mask],
84
+ unnormalized_widths=unnormalized_widths[inside_interval_mask, :],
85
+ unnormalized_heights=unnormalized_heights[inside_interval_mask, :],
86
+ unnormalized_derivatives=unnormalized_derivatives[inside_interval_mask, :],
87
+ inverse=inverse,
88
+ left=-tail_bound, right=tail_bound, bottom=-tail_bound, top=tail_bound,
89
+ min_bin_width=min_bin_width,
90
+ min_bin_height=min_bin_height,
91
+ min_derivative=min_derivative
92
+ )
93
+
94
+ return outputs, logabsdet
95
+
96
+ def rational_quadratic_spline(inputs,
97
+ unnormalized_widths,
98
+ unnormalized_heights,
99
+ unnormalized_derivatives,
100
+ inverse=False,
101
+ left=0., right=1., bottom=0., top=1.,
102
+ min_bin_width=DEFAULT_MIN_BIN_WIDTH,
103
+ min_bin_height=DEFAULT_MIN_BIN_HEIGHT,
104
+ min_derivative=DEFAULT_MIN_DERIVATIVE):
105
+ if torch.min(inputs) < left or torch.max(inputs) > right:
106
+ raise ValueError('Input to a transform is not within its domain')
107
+
108
+ num_bins = unnormalized_widths.shape[-1]
109
+
110
+ if min_bin_width * num_bins > 1.0:
111
+ raise ValueError('Minimal bin width too large for the number of bins')
112
+ if min_bin_height * num_bins > 1.0:
113
+ raise ValueError('Minimal bin height too large for the number of bins')
114
+
115
+ widths = F.softmax(unnormalized_widths, dim=-1)
116
+ widths = min_bin_width + (1 - min_bin_width * num_bins) * widths
117
+ cumwidths = torch.cumsum(widths, dim=-1)
118
+ cumwidths = F.pad(cumwidths, pad=(1, 0), mode='constant', value=0.0)
119
+ cumwidths = (right - left) * cumwidths + left
120
+ cumwidths[..., 0] = left
121
+ cumwidths[..., -1] = right
122
+ widths = cumwidths[..., 1:] - cumwidths[..., :-1]
123
+
124
+ derivatives = min_derivative + F.softplus(unnormalized_derivatives)
125
+
126
+ heights = F.softmax(unnormalized_heights, dim=-1)
127
+ heights = min_bin_height + (1 - min_bin_height * num_bins) * heights
128
+ cumheights = torch.cumsum(heights, dim=-1)
129
+ cumheights = F.pad(cumheights, pad=(1, 0), mode='constant', value=0.0)
130
+ cumheights = (top - bottom) * cumheights + bottom
131
+ cumheights[..., 0] = bottom
132
+ cumheights[..., -1] = top
133
+ heights = cumheights[..., 1:] - cumheights[..., :-1]
134
+
135
+ if inverse:
136
+ bin_idx = searchsorted(cumheights, inputs)[..., None]
137
+ else:
138
+ bin_idx = searchsorted(cumwidths, inputs)[..., None]
139
+
140
+ input_cumwidths = cumwidths.gather(-1, bin_idx)[..., 0]
141
+ input_bin_widths = widths.gather(-1, bin_idx)[..., 0]
142
+
143
+ input_cumheights = cumheights.gather(-1, bin_idx)[..., 0]
144
+ delta = heights / widths
145
+ input_delta = delta.gather(-1, bin_idx)[..., 0]
146
+
147
+ input_derivatives = derivatives.gather(-1, bin_idx)[..., 0]
148
+ input_derivatives_plus_one = derivatives[..., 1:].gather(-1, bin_idx)[..., 0]
149
+
150
+ input_heights = heights.gather(-1, bin_idx)[..., 0]
151
+
152
+ if inverse:
153
+ a = (((inputs - input_cumheights) * (input_derivatives
154
+ + input_derivatives_plus_one
155
+ - 2 * input_delta)
156
+ + input_heights * (input_delta - input_derivatives)))
157
+ b = (input_heights * input_derivatives
158
+ - (inputs - input_cumheights) * (input_derivatives
159
+ + input_derivatives_plus_one
160
+ - 2 * input_delta))
161
+ c = - input_delta * (inputs - input_cumheights)
162
+
163
+ discriminant = b.pow(2) - 4 * a * c
164
+ assert (discriminant >= 0).all()
165
+
166
+ root = (2 * c) / (-b - torch.sqrt(discriminant))
167
+ outputs = root * input_bin_widths + input_cumwidths
168
+
169
+ theta_one_minus_theta = root * (1 - root)
170
+ denominator = input_delta + ((input_derivatives + input_derivatives_plus_one - 2 * input_delta)
171
+ * theta_one_minus_theta)
172
+ derivative_numerator = input_delta.pow(2) * (input_derivatives_plus_one * root.pow(2)
173
+ + 2 * input_delta * theta_one_minus_theta
174
+ + input_derivatives * (1 - root).pow(2))
175
+ logabsdet = torch.log(derivative_numerator) - 2 * torch.log(denominator)
176
+
177
+ return outputs, -logabsdet
178
+ else:
179
+ theta = (inputs - input_cumwidths) / input_bin_widths
180
+ theta_one_minus_theta = theta * (1 - theta)
181
+
182
+ numerator = input_heights * (input_delta * theta.pow(2)
183
+ + input_derivatives * theta_one_minus_theta)
184
+ denominator = input_delta + ((input_derivatives + input_derivatives_plus_one - 2 * input_delta)
185
+ * theta_one_minus_theta)
186
+ outputs = input_cumheights + numerator / denominator
187
+
188
+ derivative_numerator = input_delta.pow(2) * (input_derivatives_plus_one * theta.pow(2)
189
+ + 2 * input_delta * theta_one_minus_theta
190
+ + input_derivatives * (1 - theta).pow(2))
191
+ logabsdet = torch.log(derivative_numerator) - 2 * torch.log(denominator)
192
+
193
+ return outputs, logabsdet
utils.py ADDED
@@ -0,0 +1,258 @@
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
1
+ import os
2
+ import glob
3
+ import sys
4
+ import argparse
5
+ import logging
6
+ import json
7
+ import subprocess
8
+ import numpy as np
9
+ from scipy.io.wavfile import read
10
+ import torch
11
+
12
+ MATPLOTLIB_FLAG = False
13
+
14
+ logging.basicConfig(stream=sys.stdout, level=logging.WARNING)
15
+ logger = logging
16
+
17
+
18
+ def load_checkpoint(checkpoint_path, model, optimizer=None):
19
+ assert os.path.isfile(checkpoint_path)
20
+ checkpoint_dict = torch.load(checkpoint_path, map_location='cpu')
21
+ iteration = checkpoint_dict['iteration']
22
+ learning_rate = checkpoint_dict['learning_rate']
23
+ if optimizer is not None:
24
+ optimizer.load_state_dict(checkpoint_dict['optimizer'])
25
+ saved_state_dict = checkpoint_dict['model']
26
+ if hasattr(model, 'module'):
27
+ state_dict = model.module.state_dict()
28
+ else:
29
+ state_dict = model.state_dict()
30
+ new_state_dict= {}
31
+ for k, v in state_dict.items():
32
+ try:
33
+ new_state_dict[k] = saved_state_dict[k]
34
+ except:
35
+ logger.info("%s is not in the checkpoint" % k)
36
+ new_state_dict[k] = v
37
+ if hasattr(model, 'module'):
38
+ model.module.load_state_dict(new_state_dict)
39
+ else:
40
+ model.load_state_dict(new_state_dict)
41
+ logger.info("Loaded checkpoint '{}' (iteration {})" .format(
42
+ checkpoint_path, iteration))
43
+ return model, optimizer, learning_rate, iteration
44
+
45
+
46
+ def save_checkpoint(model, optimizer, learning_rate, iteration, checkpoint_path):
47
+ logger.info("Saving model and optimizer state at iteration {} to {}".format(
48
+ iteration, checkpoint_path))
49
+ if hasattr(model, 'module'):
50
+ state_dict = model.module.state_dict()
51
+ else:
52
+ state_dict = model.state_dict()
53
+ torch.save({'model': state_dict,
54
+ 'iteration': iteration,
55
+ 'optimizer': optimizer.state_dict(),
56
+ 'learning_rate': learning_rate}, checkpoint_path)
57
+
58
+
59
+ def summarize(writer, global_step, scalars={}, histograms={}, images={}, audios={}, audio_sampling_rate=22050):
60
+ for k, v in scalars.items():
61
+ writer.add_scalar(k, v, global_step)
62
+ for k, v in histograms.items():
63
+ writer.add_histogram(k, v, global_step)
64
+ for k, v in images.items():
65
+ writer.add_image(k, v, global_step, dataformats='HWC')
66
+ for k, v in audios.items():
67
+ writer.add_audio(k, v, global_step, audio_sampling_rate)
68
+
69
+
70
+ def latest_checkpoint_path(dir_path, regex="G_*.pth"):
71
+ f_list = glob.glob(os.path.join(dir_path, regex))
72
+ f_list.sort(key=lambda f: int("".join(filter(str.isdigit, f))))
73
+ x = f_list[-1]
74
+ print(x)
75
+ return x
76
+
77
+
78
+ def plot_spectrogram_to_numpy(spectrogram):
79
+ global MATPLOTLIB_FLAG
80
+ if not MATPLOTLIB_FLAG:
81
+ import matplotlib
82
+ matplotlib.use("Agg")
83
+ MATPLOTLIB_FLAG = True
84
+ mpl_logger = logging.getLogger('matplotlib')
85
+ mpl_logger.setLevel(logging.WARNING)
86
+ import matplotlib.pylab as plt
87
+ import numpy as np
88
+
89
+ fig, ax = plt.subplots(figsize=(10,2))
90
+ im = ax.imshow(spectrogram, aspect="auto", origin="lower",
91
+ interpolation='none')
92
+ plt.colorbar(im, ax=ax)
93
+ plt.xlabel("Frames")
94
+ plt.ylabel("Channels")
95
+ plt.tight_layout()
96
+
97
+ fig.canvas.draw()
98
+ data = np.fromstring(fig.canvas.tostring_rgb(), dtype=np.uint8, sep='')
99
+ data = data.reshape(fig.canvas.get_width_height()[::-1] + (3,))
100
+ plt.close()
101
+ return data
102
+
103
+
104
+ def plot_alignment_to_numpy(alignment, info=None):
105
+ global MATPLOTLIB_FLAG
106
+ if not MATPLOTLIB_FLAG:
107
+ import matplotlib
108
+ matplotlib.use("Agg")
109
+ MATPLOTLIB_FLAG = True
110
+ mpl_logger = logging.getLogger('matplotlib')
111
+ mpl_logger.setLevel(logging.WARNING)
112
+ import matplotlib.pylab as plt
113
+ import numpy as np
114
+
115
+ fig, ax = plt.subplots(figsize=(6, 4))
116
+ im = ax.imshow(alignment.transpose(), aspect='auto', origin='lower',
117
+ interpolation='none')
118
+ fig.colorbar(im, ax=ax)
119
+ xlabel = 'Decoder timestep'
120
+ if info is not None:
121
+ xlabel += '\n\n' + info
122
+ plt.xlabel(xlabel)
123
+ plt.ylabel('Encoder timestep')
124
+ plt.tight_layout()
125
+
126
+ fig.canvas.draw()
127
+ data = np.fromstring(fig.canvas.tostring_rgb(), dtype=np.uint8, sep='')
128
+ data = data.reshape(fig.canvas.get_width_height()[::-1] + (3,))
129
+ plt.close()
130
+ return data
131
+
132
+
133
+ def load_wav_to_torch(full_path):
134
+ sampling_rate, data = read(full_path)
135
+ return torch.FloatTensor(data.astype(np.float32)), sampling_rate
136
+
137
+
138
+ def load_filepaths_and_text(filename, split="|"):
139
+ with open(filename, encoding='utf-8') as f:
140
+ filepaths_and_text = [line.strip().split(split) for line in f]
141
+ return filepaths_and_text
142
+
143
+
144
+ def get_hparams(init=True):
145
+ parser = argparse.ArgumentParser()
146
+ parser.add_argument('-c', '--config', type=str, default="./configs/base.json",
147
+ help='JSON file for configuration')
148
+ parser.add_argument('-m', '--model', type=str, required=True,
149
+ help='Model name')
150
+
151
+ args = parser.parse_args()
152
+ model_dir = os.path.join("./logs", args.model)
153
+
154
+ if not os.path.exists(model_dir):
155
+ os.makedirs(model_dir)
156
+
157
+ config_path = args.config
158
+ config_save_path = os.path.join(model_dir, "config.json")
159
+ if init:
160
+ with open(config_path, "r") as f:
161
+ data = f.read()
162
+ with open(config_save_path, "w") as f:
163
+ f.write(data)
164
+ else:
165
+ with open(config_save_path, "r") as f:
166
+ data = f.read()
167
+ config = json.loads(data)
168
+
169
+ hparams = HParams(**config)
170
+ hparams.model_dir = model_dir
171
+ return hparams
172
+
173
+
174
+ def get_hparams_from_dir(model_dir):
175
+ config_save_path = os.path.join(model_dir, "config.json")
176
+ with open(config_save_path, "r") as f:
177
+ data = f.read()
178
+ config = json.loads(data)
179
+
180
+ hparams =HParams(**config)
181
+ hparams.model_dir = model_dir
182
+ return hparams
183
+
184
+
185
+ def get_hparams_from_file(config_path):
186
+ with open(config_path, "r") as f:
187
+ data = f.read()
188
+ config = json.loads(data)
189
+
190
+ hparams =HParams(**config)
191
+ return hparams
192
+
193
+
194
+ def check_git_hash(model_dir):
195
+ source_dir = os.path.dirname(os.path.realpath(__file__))
196
+ if not os.path.exists(os.path.join(source_dir, ".git")):
197
+ logger.warn("{} is not a git repository, therefore hash value comparison will be ignored.".format(
198
+ source_dir
199
+ ))
200
+ return
201
+
202
+ cur_hash = subprocess.getoutput("git rev-parse HEAD")
203
+
204
+ path = os.path.join(model_dir, "githash")
205
+ if os.path.exists(path):
206
+ saved_hash = open(path).read()
207
+ if saved_hash != cur_hash:
208
+ logger.warn("git hash values are different. {}(saved) != {}(current)".format(
209
+ saved_hash[:8], cur_hash[:8]))
210
+ else:
211
+ open(path, "w").write(cur_hash)
212
+
213
+
214
+ def get_logger(model_dir, filename="train.log"):
215
+ global logger
216
+ logger = logging.getLogger(os.path.basename(model_dir))
217
+ logger.setLevel(logging.DEBUG)
218
+
219
+ formatter = logging.Formatter("%(asctime)s\t%(name)s\t%(levelname)s\t%(message)s")
220
+ if not os.path.exists(model_dir):
221
+ os.makedirs(model_dir)
222
+ h = logging.FileHandler(os.path.join(model_dir, filename))
223
+ h.setLevel(logging.DEBUG)
224
+ h.setFormatter(formatter)
225
+ logger.addHandler(h)
226
+ return logger
227
+
228
+
229
+ class HParams():
230
+ def __init__(self, **kwargs):
231
+ for k, v in kwargs.items():
232
+ if type(v) == dict:
233
+ v = HParams(**v)
234
+ self[k] = v
235
+
236
+ def keys(self):
237
+ return self.__dict__.keys()
238
+
239
+ def items(self):
240
+ return self.__dict__.items()
241
+
242
+ def values(self):
243
+ return self.__dict__.values()
244
+
245
+ def __len__(self):
246
+ return len(self.__dict__)
247
+
248
+ def __getitem__(self, key):
249
+ return getattr(self, key)
250
+
251
+ def __setitem__(self, key, value):
252
+ return setattr(self, key, value)
253
+
254
+ def __contains__(self, key):
255
+ return key in self.__dict__
256
+
257
+ def __repr__(self):
258
+ return self.__dict__.__repr__()