


"""simple docstring"""




import numpy as np

import datasets


lowercase_			     =  "\nCompute the Mahalanobis Distance\n\nMahalonobis distance is the distance between a point and a distribution.\nAnd not between two distinct points. It is effectively a multivariate equivalent of the Euclidean distance.\nIt was introduced by Prof. P. C. Mahalanobis in 1936\nand has been used in various statistical applications ever since\n[source: https://www.machinelearningplus.com/statistics/mahalanobis-distance/]\n"

lowercase_			     =  "\\n@article{de2000mahalanobis,\n  title={The mahalanobis distance},\n  author={De Maesschalck, Roy and Jouan-Rimbaud, Delphine and Massart, D{\'e}sir{\'e} L},\n  journal={Chemometrics and intelligent laboratory systems},\n  volume={50},\n  number={1},\n  pages={1--18},\n  year={2000},\n  publisher={Elsevier}\n}\n"

lowercase_			     =  "\nArgs:\n    X: List of datapoints to be compared with the `reference_distribution`.\n    reference_distribution: List of datapoints from the reference distribution we want to compare to.\nReturns:\n    mahalanobis: The Mahalonobis distance for each datapoint in `X`.\nExamples:\n\n    >>> mahalanobis_metric = datasets.load_metric(\"mahalanobis\")\n    >>> results = mahalanobis_metric.compute(reference_distribution=[[0, 1], [1, 0]], X=[[0, 1]])\n    >>> print(results)\n    {'mahalanobis': array([0.5])}\n"
@datasets.utils.file_utils.add_start_docstrings(_DESCRIPTION     ,			_KWARGS_DESCRIPTION       )
class 					__lowerCAmelCase					(   datasets.Metric       ):


		'''simple docstring'''



		def 			__UpperCAmelCase       (  self      ):
				return datasets.MetricInfo(
				    description=_DESCRIPTION							, citation=_CITATION							, inputs_description=_KWARGS_DESCRIPTION							, features=datasets.Features(
				        {
				            '''X''': datasets.Sequence(datasets.Value('''float'''							, id='''sequence'''      )							, id='''X'''      ),
				        }      )							, )





		def 			__UpperCAmelCase       (  self							, _a							, _a      ):
				# convert to numpy arrays
				__a           =			np.array(_a      )
				__a           =			np.array(_a      )

				# Assert that arrays are 2D
				if len(X.shape      ) != 2:
						raise ValueError('''Expected `X` to be a 2D vector'''      )
				if len(reference_distribution.shape      ) != 2:
						raise ValueError('''Expected `reference_distribution` to be a 2D vector'''      )
				if reference_distribution.shape[0] < 2:
						raise ValueError(
						    '''Expected `reference_distribution` to be a 2D vector with more than one element in the first dimension'''      )

				# Get mahalanobis distance for each prediction
				__a           =			X - np.mean(_a      )
				__a           =			np.cov(reference_distribution.T      )
				try:
						__a           =			np.linalg.inv(_a      )
				except np.linalg.LinAlgError:
						__a           =			np.linalg.pinv(_a      )
				__a           =			np.dot(_a							, _a      )
				__a           =			np.dot(_a							, X_minus_mu.T      ).diagonal()

				return {"mahalanobis": mahal_dist}










"""simple docstring"""




lowercase_			     =  "ABCDEFGHIJKLMNOPQRSTUVWXYZabcdefghijklmnopqrstuvwxyz0123456789+/"

def        lowercase							(		lowerCAmelCase__ :     bytes			) ->      bytes:
		# Make sure the supplied data is a bytes-like object
		if not isinstance(lowerCAmelCase__    ,					lowerCAmelCase__			):
				__a           =			f'''a bytes-like object is required, not \'{data.__class__.__name__}\''''
				raise TypeError(lowerCAmelCase__			)

		__a           =			''''''.join(bin(lowerCAmelCase__			)[2:].zfill(8			) for byte in data			)

		__a           =			len(lowerCAmelCase__			) % 6 != 0

		if padding_needed:
				# The padding that will be added later
				__a           =			b'''=''' * ((6 - len(lowerCAmelCase__			) % 6) // 2)

				# Append binary_stream with arbitrary binary digits (0's by default) to make its
				# length a multiple of 6.
				binary_stream += "0" * (6 - len(lowerCAmelCase__			) % 6)
		else:
				__a           =			b''''''

		# Encode every 6 binary digits to their corresponding Base64 character
		return (
		    "".join(
		        B64_CHARSET[int(binary_stream[index : index + 6]    ,					2			)]
		        for index in range(0    ,					len(lowerCAmelCase__			)    ,					6			)			).encode()
		    + padding
		)

def        lowercase							(		lowerCAmelCase__ :     str			) ->      bytes:
		# Make sure encoded_data is either a string or a bytes-like object
		if not isinstance(lowerCAmelCase__    ,					lowerCAmelCase__			) and not isinstance(lowerCAmelCase__    ,					lowerCAmelCase__			):
				__a           =			(
				    '''argument should be a bytes-like object or ASCII string, '''
				    f'''not \'{encoded_data.__class__.__name__}\''''
				)
				raise TypeError(lowerCAmelCase__			)

		# In case encoded_data is a bytes-like object, make sure it contains only
		# ASCII characters so we convert it to a string object
		if isinstance(lowerCAmelCase__    ,					lowerCAmelCase__			):
				try:
						__a           =			encoded_data.decode('''utf-8'''			)
				except UnicodeDecodeError:
						raise ValueError('''base64 encoded data should only contain ASCII characters'''			)

		__a           =			encoded_data.count('''='''			)

		# Check if the encoded string contains non base64 characters
		if padding:
				assert all(
				    char in B64_CHARSET for char in encoded_data[:-padding]			), "Invalid base64 character(s) found."
		else:
				assert all(
				    char in B64_CHARSET for char in encoded_data			), "Invalid base64 character(s) found."

		# Check the padding
		assert len(lowerCAmelCase__			) % 4 == 0 and padding < 3, "Incorrect padding"

		if padding:
				# Remove padding if there is one
				__a           =			encoded_data[:-padding]

				__a           =			''''''.join(
				    bin(B64_CHARSET.index(lowerCAmelCase__			)			)[2:].zfill(6			) for char in encoded_data			)[: -padding * 2]
		else:
				__a           =			''''''.join(
				    bin(B64_CHARSET.index(lowerCAmelCase__			)			)[2:].zfill(6			) for char in encoded_data			)

		__a           =			[
		    int(binary_stream[index : index + 8]    ,					2			)
		    for index in range(0    ,					len(lowerCAmelCase__			)    ,					8			)
		]

		return bytes(lowerCAmelCase__			)


if __name__ == "__main__":
				import doctest

				doctest.testmod()






