




'''simple docstring'''


import argparse
import logging
import os
from datetime import datetime

import numpy as np
import torch
from torch import nn
from torch.utils.data import DataLoader, RandomSampler, TensorDataset
from tqdm import tqdm

from transformers import GPTaLMHeadModel


_snake_case				:				str       							=		logging.getLogger(__name__)
def 		snake_case_   (UpperCamelCase	:  int			,						UpperCamelCase	:  List[str]				):





							'''simple docstring'''


							if os.path.exists(UpperCamelCase				):
														if os.path.exists(os.path.join(UpperCamelCase			,						'''config.json'''				)				) and os.path.isfile(
														    os.path.join(UpperCamelCase			,						'''config.json'''				)				):
																					os.remove(os.path.join(UpperCamelCase			,						'''config.json'''				)				)
														if os.path.exists(os.path.join(UpperCamelCase			,						'''pytorch_model.bin'''				)				) and os.path.isfile(
														    os.path.join(UpperCamelCase			,						'''pytorch_model.bin'''				)				):
																					os.remove(os.path.join(UpperCamelCase			,						'''pytorch_model.bin'''				)				)
							else:
														os.makedirs(UpperCamelCase				)
							model.save_pretrained(UpperCamelCase				)
def 		snake_case_   (UpperCamelCase	:  Tuple			,						UpperCamelCase	:  Tuple=False				):





							'''simple docstring'''


							_a     					=    2
							if unlogit:
														_a     					=    torch.pow(UpperCamelCase			,						UpperCamelCase				)
							_a     					=    p * torch.log(UpperCamelCase				)
							_a     					=    0
							return -plogp.sum(dim=-1				)
def 		snake_case_   (UpperCamelCase	:  List[str]				):





							'''simple docstring'''


							logger.info('''lv, h >\t''' + '''\t'''.join(f'{x + 1}' for x in range(len(UpperCamelCase				)				)				)				)
							for row in range(len(UpperCamelCase				)				):
														if tensor.dtype != torch.long:
																					logger.info(f'layer {row + 1}:\t' + '''\t'''.join(f'{x:.5f}' for x in tensor[row].cpu().data				)				)
														else:
																					logger.info(f'layer {row + 1}:\t' + '''\t'''.join(f'{x:d}' for x in tensor[row].cpu().data				)				)
def 		snake_case_   (UpperCamelCase	:  Any			,						UpperCamelCase	:  Tuple			,						UpperCamelCase	:  Dict			,						UpperCamelCase	:  Union[str, Any]=True			,						UpperCamelCase	:  Tuple=True			,						UpperCamelCase	:  Optional[Any]=None			,						UpperCamelCase	:  Optional[Any]=False				):





							'''simple docstring'''


							_a					,							_a     					=    model.config.num_hidden_layers, model.config.num_attention_heads
							_a     					=    torch.zeros(UpperCamelCase			,						UpperCamelCase				).to(args.device				)
							_a     					=    torch.zeros(UpperCamelCase			,						UpperCamelCase				).to(args.device				)

							if head_mask is None:
														_a     					=    torch.ones(UpperCamelCase			,						UpperCamelCase				).to(args.device				)

							head_mask.requires_grad_(requires_grad=UpperCamelCase				)
							# If actually pruned attention multi-head, set head mask to None to avoid shape mismatch
							if actually_pruned:
														_a     					=    None

							_a     					=    0.0
							_a     					=    0.0
							for step, inputs in enumerate(tqdm(UpperCamelCase			,						desc='''Iteration'''			,						disable=args.local_rank not in [-1, 0]				)				):
														_a     					=    tuple(t.to(args.device				) for t in inputs				)
														((_a)					,							)     					=    inputs

														# Do a forward pass (not with torch.no_grad() since we need gradients for importance score - see below)
														_a     					=    model(UpperCamelCase			,						labels=UpperCamelCase			,						head_mask=UpperCamelCase				)
														#  (loss), lm_logits, presents, (all hidden_states), (attentions)
														_a					,							_a					,							_a     					=    (
														    outputs[0],
														    outputs[1],
														    outputs[-1],
														)  # Loss and logits are the first, attention the last
														loss.backward()  # Backpropagate to populate the gradients in the head mask
														total_loss += loss.detach().cpu().numpy()
														if compute_entropy:
																					for layer, attn in enumerate(UpperCamelCase				):
																												_a     					=    entropy(attn.detach()			,						UpperCamelCase				)
																												attn_entropy[layer] += masked_entropy.sum(-1				).sum(0				).sum(0				).detach()

														if compute_importance:
																					head_importance += head_mask.grad.abs().detach()
														tot_tokens += torch.ones_like(UpperCamelCase				).float().detach().sum().data

							# Normalize
							attn_entropy /= tot_tokens
							head_importance /= tot_tokens
							# Layerwise importance normalization
							if not args.dont_normalize_importance_by_layer:
														_a     					=    2
														_a     					=    torch.pow(torch.pow(UpperCamelCase			,						UpperCamelCase				).sum(-1				)			,						1 / exponent				)
														head_importance /= norm_by_layer.unsqueeze(-1				) + 1e-20

							if not args.dont_normalize_global_importance:
														_a     					=    (head_importance - head_importance.min()) / (head_importance.max() - head_importance.min())

							# Print matrices
							if compute_entropy:
														logger.info('''Attention entropies'''				)
														print_ad_tensor(UpperCamelCase				)
							if compute_importance:
														logger.info('''Head importance scores'''				)
														print_ad_tensor(UpperCamelCase				)
							logger.info('''Head ranked by importance scores'''				)
							_a     					=    torch.zeros(head_importance.numel()			,						dtype=torch.long			,						device=args.device				)
							_a     					=    torch.arange(
							    head_importance.numel()			,						device=args.device				)
							_a     					=    head_ranks.view_as(UpperCamelCase				)
							print_ad_tensor(UpperCamelCase				)
							return attn_entropy, head_importance, total_loss
def 		snake_case_   (UpperCamelCase	:  str			,						UpperCamelCase	:  int			,						UpperCamelCase	:  Any				):





							'''simple docstring'''


							_a					,							_a					,							_a     					=    compute_heads_importance(UpperCamelCase			,						UpperCamelCase			,						UpperCamelCase			,						compute_entropy=UpperCamelCase				)
							_a     					=    1 / loss  # instead of downsteam score use the LM loss
							logger.info('''Pruning: original score: %f, threshold: %f'''			,						UpperCamelCase			,						original_score * args.masking_threshold				)

							_a     					=    torch.ones_like(UpperCamelCase				)
							_a     					=    max(1			,						int(new_head_mask.numel() * args.masking_amount				)				)

							_a     					=    original_score
							while current_score >= original_score * args.masking_threshold:
														_a     					=    new_head_mask.clone().detach()  # save current head mask
														# heads from least important to most - keep only not-masked heads
														_a     					=    float('''Inf'''				)
														_a     					=    head_importance.view(-1				).sort()[1]

														if len(UpperCamelCase				) <= num_to_mask:
																					print('''BREAK BY num_to_mask'''				)
																					break

														# mask heads
														_a     					=    current_heads_to_mask[:num_to_mask]
														logger.info('''Heads to mask: %s'''			,						str(current_heads_to_mask.tolist()				)				)
														_a     					=    new_head_mask.view(-1				)
														_a     					=    0.0
														_a     					=    new_head_mask.view_as(UpperCamelCase				)
														_a     					=    new_head_mask.clone().detach()
														print_ad_tensor(UpperCamelCase				)

														# Compute metric and head importance again
														_a					,							_a					,							_a     					=    compute_heads_importance(
														    UpperCamelCase			,						UpperCamelCase			,						UpperCamelCase			,						compute_entropy=UpperCamelCase			,						head_mask=UpperCamelCase				)
														_a     					=    1 / loss
														logger.info(
														    '''Masking: current score: %f, remaining heads %d (%.1f percents)'''			,						UpperCamelCase			,						new_head_mask.sum()			,						new_head_mask.sum() / new_head_mask.numel() * 100			,						)

							logger.info('''Final head mask'''				)
							print_ad_tensor(UpperCamelCase				)
							np.save(os.path.join(args.output_dir			,						'''head_mask.npy'''				)			,						head_mask.detach().cpu().numpy()				)

							return head_mask
def 		snake_case_   (UpperCamelCase	:  int			,						UpperCamelCase	:  Dict			,						UpperCamelCase	:  Optional[Any]			,						UpperCamelCase	:  List[str]				):





							'''simple docstring'''


							_a     					=    datetime.now()
							_a					,							_a					,							_a     					=    compute_heads_importance(
							    UpperCamelCase			,						UpperCamelCase			,						UpperCamelCase			,						compute_entropy=UpperCamelCase			,						compute_importance=UpperCamelCase			,						head_mask=UpperCamelCase				)
							_a     					=    1 / loss
							_a     					=    datetime.now() - before_time

							_a     					=    sum(p.numel() for p in model.parameters()				)
							_a     					=    {
							    layer: (1 - head_mask[layer].long()).nonzero().squeeze().tolist() for layer in range(len(UpperCamelCase				)				)
							}

							for k, v in heads_to_prune.items():
														if isinstance(UpperCamelCase			,						UpperCamelCase				):
																					_a     					=    [
																					    v,
																					]

							assert sum(len(UpperCamelCase				) for h in heads_to_prune.values()				) == (1 - head_mask.long()).sum().item()
							model.prune_heads(UpperCamelCase				)
							_a     					=    sum(p.numel() for p in model.parameters()				)

							_a     					=    datetime.now()
							_a					,							_a					,							_a     					=    compute_heads_importance(
							    UpperCamelCase			,						UpperCamelCase			,						UpperCamelCase			,						compute_entropy=UpperCamelCase			,						compute_importance=UpperCamelCase			,						head_mask=UpperCamelCase			,						actually_pruned=UpperCamelCase			,						)

							_a     					=    1 / loss
							_a     					=    datetime.now() - before_time

							logger.info(
							    '''Pruning: original num of params: %.2e, after pruning %.2e (%.1f percents)'''			,						UpperCamelCase			,						UpperCamelCase			,						pruned_num_params / original_num_params * 100			,						)
							logger.info('''Pruning: score with masking: %f score with pruning: %f'''			,						UpperCamelCase			,						UpperCamelCase				)
							logger.info('''Pruning: speed ratio (original timing / new timing): %f percents'''			,						original_time / new_time * 100				)
							save_model(UpperCamelCase			,						args.output_dir				)
def 		snake_case_   ():





							'''simple docstring'''


							_a     					=    argparse.ArgumentParser()
							# Required parameters
							parser.add_argument(
							    '''--data_dir'''			,						default=UpperCamelCase			,						type=UpperCamelCase			,						required=UpperCamelCase			,						help='''The input data dir. Should contain the .tsv files (or other data files) for the task.'''			,						)
							parser.add_argument(
							    '''--model_name_or_path'''			,						default=UpperCamelCase			,						type=UpperCamelCase			,						required=UpperCamelCase			,						help='''Path to pretrained model or model identifier from huggingface.co/models'''			,						)
							parser.add_argument(
							    '''--output_dir'''			,						default=UpperCamelCase			,						type=UpperCamelCase			,						required=UpperCamelCase			,						help='''The output directory where the model predictions and checkpoints will be written.'''			,						)

							# Other parameters
							parser.add_argument(
							    '''--config_name'''			,						default=''''''			,						type=UpperCamelCase			,						help='''Pretrained config name or path if not the same as model_name_or_path'''			,						)
							parser.add_argument(
							    '''--tokenizer_name'''			,						default=''''''			,						type=UpperCamelCase			,						help='''Pretrained tokenizer name or path if not the same as model_name_or_path'''			,						)
							parser.add_argument(
							    '''--cache_dir'''			,						default=UpperCamelCase			,						type=UpperCamelCase			,						help='''Where do you want to store the pre-trained models downloaded from s3'''			,						)
							parser.add_argument(
							    '''--data_subset'''			,						type=UpperCamelCase			,						default=-1			,						help='''If > 0: limit the data to a subset of data_subset instances.'''				)
							parser.add_argument(
							    '''--overwrite_output_dir'''			,						action='''store_true'''			,						help='''Whether to overwrite data in output directory'''				)
							parser.add_argument(
							    '''--overwrite_cache'''			,						action='''store_true'''			,						help='''Overwrite the cached training and evaluation sets'''				)

							parser.add_argument(
							    '''--dont_normalize_importance_by_layer'''			,						action='''store_true'''			,						help='''Don\'t normalize importance score by layers'''				)
							parser.add_argument(
							    '''--dont_normalize_global_importance'''			,						action='''store_true'''			,						help='''Don\'t normalize all importance scores between 0 and 1'''			,						)

							parser.add_argument(
							    '''--try_masking'''			,						action='''store_true'''			,						help='''Whether to try to mask head until a threshold of accuracy.'''				)
							parser.add_argument(
							    '''--masking_threshold'''			,						default=0.9			,						type=UpperCamelCase			,						help='''masking threshold in term of metrics (stop masking when metric < threshold * original metric value).'''			,						)
							parser.add_argument(
							    '''--masking_amount'''			,						default=0.1			,						type=UpperCamelCase			,						help='''Amount to heads to masking at each masking step.'''				)
							parser.add_argument('''--metric_name'''			,						default='''acc'''			,						type=UpperCamelCase			,						help='''Metric to use for head masking.'''				)

							parser.add_argument(
							    '''--max_seq_length'''			,						default=128			,						type=UpperCamelCase			,						help=(
							        '''The maximum total input sequence length after WordPiece tokenization. \n'''
							        '''Sequences longer than this will be truncated, sequences shorter padded.'''
							    )			,						)
							parser.add_argument('''--batch_size'''			,						default=1			,						type=UpperCamelCase			,						help='''Batch size.'''				)

							parser.add_argument('''--seed'''			,						type=UpperCamelCase			,						default=42				)
							parser.add_argument('''--local_rank'''			,						type=UpperCamelCase			,						default=-1			,						help='''local_rank for distributed training on gpus'''				)
							parser.add_argument('''--no_cuda'''			,						action='''store_true'''			,						help='''Whether not to use CUDA when available'''				)
							parser.add_argument('''--server_ip'''			,						type=UpperCamelCase			,						default=''''''			,						help='''Can be used for distant debugging.'''				)
							parser.add_argument('''--server_port'''			,						type=UpperCamelCase			,						default=''''''			,						help='''Can be used for distant debugging.'''				)
							_a     					=    parser.parse_args()

							if args.server_ip and args.server_port:
														# Distant debugging - see https://code.visualstudio.com/docs/python/debugging#_attach-to-a-local-script
														import ptvsd

														print('''Waiting for debugger attach'''				)
														ptvsd.enable_attach(address=(args.server_ip, args.server_port)			,						redirect_output=UpperCamelCase				)
														ptvsd.wait_for_attach()

							# Setup devices and distributed training
							if args.local_rank == -1 or args.no_cuda:
														_a     					=    torch.device('''cuda''' if torch.cuda.is_available() and not args.no_cuda else '''cpu'''				)
														_a     					=    0 if args.no_cuda else torch.cuda.device_count()
							else:
														torch.cuda.set_device(args.local_rank				)
														_a     					=    torch.device('''cuda'''			,						args.local_rank				)
														_a     					=    1
														torch.distributed.init_process_group(backend='''nccl'''				)  # Initializes the distributed backend

							# Setup logging
							logging.basicConfig(level=logging.INFO if args.local_rank in [-1, 0] else logging.WARN				)
							logger.info('''device: {} n_gpu: {}, distributed: {}'''.format(args.device			,						args.n_gpu			,						bool(args.local_rank != -1				)				)				)

							_a     					=    GPTaLMHeadModel.from_pretrained(args.model_name_or_path				)

							# Distributed and parallel training
							model.to(args.device				)
							if args.local_rank != -1:
														_a     					=    nn.parallel.DistributedDataParallel(
														    UpperCamelCase			,						device_ids=[args.local_rank]			,						output_device=args.local_rank			,						find_unused_parameters=UpperCamelCase				)
							elif args.n_gpu > 1:
														_a     					=    nn.DataParallel(UpperCamelCase				)

							# Print/save training arguments
							os.makedirs(args.output_dir			,						exist_ok=UpperCamelCase				)
							torch.save(UpperCamelCase			,						os.path.join(args.output_dir			,						'''run_args.bin'''				)				)
							logger.info('''Training/evaluation parameters %s'''			,						UpperCamelCase				)

							# Prepare dataset
							_a     					=    np.concatenate(
							    [
							        np.loadtxt(args.data_dir			,						dtype=np.intaa				),
							    ]				)
							_a     					=    (torch.from_numpy(UpperCamelCase				),)
							_a     					=    TensorDataset(*UpperCamelCase				)
							_a     					=    RandomSampler(UpperCamelCase				)
							_a     					=    DataLoader(UpperCamelCase			,						sampler=UpperCamelCase			,						batch_size=args.batch_size				)

							# Compute head entropy and importance score
							compute_heads_importance(UpperCamelCase			,						UpperCamelCase			,						UpperCamelCase				)

							# Try head masking (set heads to zero until the score goes under a threshole)
							# and head pruning (remove masked heads and see the effect on the network)
							if args.try_masking and args.masking_threshold > 0.0 and args.masking_threshold < 1.0:
														_a     					=    mask_heads(UpperCamelCase			,						UpperCamelCase			,						UpperCamelCase				)
														prune_heads(UpperCamelCase			,						UpperCamelCase			,						UpperCamelCase			,						UpperCamelCase				)


if __name__ == "__main__":
				main()







'''simple docstring'''


def 		snake_case_   (UpperCamelCase	:  str			,						UpperCamelCase	:  Any				):





							'''simple docstring'''


							return (pointa[0] - pointa[0]) ** 2 + (pointa[1] - pointa[1]) ** 2
def 		snake_case_   (UpperCamelCase	:  Any			,						UpperCamelCase	:  str=0				):





							'''simple docstring'''


							return sorted(UpperCamelCase			,						key=lambda UpperCamelCase				: x[column]				)
def 		snake_case_   (UpperCamelCase	:  Optional[int]			,						UpperCamelCase	:  Any			,						UpperCamelCase	:  Union[str, Any]=float('''inf'''				)				):





							'''simple docstring'''


							for i in range(points_counts - 1				):
														for j in range(i + 1			,						UpperCamelCase				):
																					_a     					=    euclidean_distance_sqr(points[i]			,						points[j]				)
																					if current_dis < min_dis:
																												_a     					=    current_dis
							return min_dis
def 		snake_case_   (UpperCamelCase	:  int			,						UpperCamelCase	:  Tuple			,						UpperCamelCase	:  List[str]=float('''inf'''				)				):





							'''simple docstring'''


							for i in range(min(6			,						points_counts - 1				)			,						UpperCamelCase				):
														for j in range(max(0			,						i - 6				)			,						UpperCamelCase				):
																					_a     					=    euclidean_distance_sqr(points[i]			,						points[j]				)
																					if current_dis < min_dis:
																												_a     					=    current_dis
							return min_dis
def 		snake_case_   (UpperCamelCase	:  int			,						UpperCamelCase	:  List[Any]			,						UpperCamelCase	:  int				):





							'''simple docstring'''


							if points_counts <= 3:
														return dis_between_closest_pair(UpperCamelCase			,						UpperCamelCase				)

							# recursion
							_a     					=    points_counts // 2
							_a     					=    closest_pair_of_points_sqr(
							    UpperCamelCase			,						points_sorted_on_y[:mid]			,						UpperCamelCase				)
							_a     					=    closest_pair_of_points_sqr(
							    UpperCamelCase			,						points_sorted_on_y[mid:]			,						points_counts - mid				)
							_a     					=    min(UpperCamelCase			,						UpperCamelCase				)
							_a     					=    []
							for point in points_sorted_on_x:
														if abs(point[0] - points_sorted_on_x[mid][0]				) < closest_pair_dis:
																					cross_strip.append(UpperCamelCase				)

							_a     					=    dis_between_closest_in_strip(
							    UpperCamelCase			,						len(UpperCamelCase				)			,						UpperCamelCase				)
							return min(UpperCamelCase			,						UpperCamelCase				)
def 		snake_case_   (UpperCamelCase	:  Optional[int]			,						UpperCamelCase	:  List[str]				):





							'''simple docstring'''


							_a     					=    column_based_sort(UpperCamelCase			,						column=0				)
							_a     					=    column_based_sort(UpperCamelCase			,						column=1				)
							return (
							    closest_pair_of_points_sqr(
							        UpperCamelCase			,						UpperCamelCase			,						UpperCamelCase				)
							) ** 0.5


if __name__ == "__main__":
				_snake_case				:				int       							=		[(2, 3), (12, 30), (40, 50), (5, 1), (12, 10), (3, 4)]
				print('Distance:', closest_pair_of_points(points, len(points)))

