t4a-faer 0.25.1

Temporary tensor4all fork of faer: an AI-driven experiment build for tenferro-rs.
Documentation 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68

// Private test adapters. Copyright (c) tensor4all contributors. SPDX-License-Identifier: MIT
fn schur_call(
	h: MatMut<'_, Scalar>,
	q: Option<MatMut<'_, Scalar>>,
	mut wr: ColMut<'_, Real>,
	mut wi: ColMut<'_, Real>,
	par: Par,
	stack: &mut MemStack,
	p: Spec<SchurParams, Scalar>,
) -> Result<(), EvdError> {
	let n = h.nrows();
	let mut w = Col::<Scalar>::full(n, scalar(Real::NAN, Real::NAN));
	hessenberg_schur_cplx(h, q, w.as_mut(), par, stack, p)?;
	for i in 0..n {
		wr[i] = real(&w[i]);
		wi[i] = imag(&w[i]);
	}
	Ok(())
}
fn vectors_call(
	t: MatRef<'_, Scalar>,
	q: Option<MatRef<'_, Scalar>>,
	x: MatMut<'_, Scalar>,
	par: Par,
	stack: &mut MemStack,
	p: Spec<EvdFromSchurParams, Scalar>,
) -> Result<(), EvdError> {
	eigenvectors_from_schur_cplx(t, q, x, par, stack, p)
}
fn evd_call(
	a: MatRef<'_, Scalar>,
	mut wr: ColMut<'_, Real>,
	mut wi: ColMut<'_, Real>,
	left: Option<MatMut<'_, Scalar>>,
	right: Option<MatMut<'_, Scalar>>,
	par: Par,
	stack: &mut MemStack,
	p: Spec<EvdParams, Scalar>,
) -> Result<(), EvdError> {
	let n = a.nrows();
	let mut w = Col::<Scalar>::full(n, scalar(Real::NAN, Real::NAN));
	evd_cplx(a, w.as_mut().as_diagonal_mut(), left, right, par, stack, p)?;
	for i in 0..n {
		wr[i] = real(&w[i]);
		wi[i] = imag(&w[i]);
	}
	Ok(())
}
fn reorder_call(
	t: MatMut<'_, Scalar>,
	q: Option<MatMut<'_, Scalar>>,
	from: usize,
	to: usize,
	par: Par,
	stack: &mut MemStack,
) -> Result<usize, SchurReorderError> {
	Ok(reorder_schur_cplx(t, q, from, to, par, stack))
}

fn adjacent_call(
	t: MatMut<'_, Scalar>,
	q: Option<MatMut<'_, Scalar>>,
	start: usize,
	_n1: usize,
	_n2: usize,
) -> isize {
	schur::complex_schur::schur_swap(t, q, start)
}