feat: add C implementation for stats/base/dists/levy/median

[vivekmaurya001]

Resolves #3727.

Description

• adds C implementation for @stdlib/stats/base/dists/levy/median along with relevant tests, benchmarks and examples

This pull request:

• Progresses: [RFC]: Add C implementation for @stdlib/stats/base/dists/levy/median #3727

Questions

Any questions for reviewers of this pull request?

No.

Other

Any other information relevant to this pull request? This may include screenshots, references, and/or implementation notes.

No.

Checklist

Please ensure the following tasks are completed before submitting this pull request.

• Read, understood, and followed the contributing guidelines

---

@stdlib-js/reviewers

[stdlib-bot]

Coverage Report

Package	Statements	Branches	Functions	Lines
stats/base/dists/levy/median	$\color{green}182/182$ $\color{green}+100.00\%$	$\color{green}9/9$ $\color{green}+100.00\%$	$\color{green}2/2$ $\color{green}+100.00\%$	$\color{green}182/182$ $\color{green}+100.00\%$

The above coverage report was generated for the changes in this PR.

[Planeshifter]

/stdlib update-copyright-years

[Planeshifter]

@vivekmaurya001 We should pre-compute this constant as a static const in the outer scope since erfcinv(0.5) is constant and doesn't change with the inputs.

[vivekmaurya001]

ook

[Planeshifter]

Suggested change

	#include "stdlib/math/base/special/pow.h"
	/**
	#include "stdlib/math/base/special/pow.h"
	/**

[Planeshifter]

Suggested change

	Evaluates the [median][median] for a [Lévy][levy-distribution] distribution with location parameter `mu` and scale parameter `c`.
	Returns the [median][median] for a [Lévy][levy-distribution] distribution with location parameter `mu` and scale parameter `c`.

[Planeshifter]

Suggested change

	mu = random_uniform( 0.0, 10.0 ) - 5.0;
	mu = random_uniform( -5.0, 5.0 );

[Planeshifter]

Suggested change

	mu = random_uniform( 0.0, 10.0 ) - 5.0;
	mu = random_uniform( -5.0, 5.0 );

[Planeshifter]

Suggested change

	mu[ i ] = random_uniform( 0.0, 100.0 ) - 50.0;
	c[ i ] = random_uniform( 0.0, 20.0 ) + STDLIB_CONSTANT_FLOAT64_EPS;
	mu[ i ] = random_uniform( -50.0, 50.0 );
	c[ i ] = random_uniform( STDLIB_CONSTANT_FLOAT64_EPS, 20.0 );

[Planeshifter]

Left a comment on how to optimize the source implementation and pointed out a bunch of other small issues. Should be ready to land once these are addressed.