pub fn eigen_frule<T, C>(
ctx: &mut C,
tensor: &Tensor<T>,
tangent: &Tensor<T>,
) -> AdResult<(EigenResult<T, T::Real>, EigenResult<T, T::Real>)>where
T: KernelLinalgScalar + Conjugate,
T::Real: KernelLinalgScalar<Real = T::Real> + Float,
C: TensorLinalgContextFor<T>,
C::Backend: 'static,Expand description
Forward-mode AD rule for eigendecomposition (JVP / pushforward).
ยงExamples
use tenferro_linalg::eigen_frule;
use tenferro_prims::CpuContext;
use tenferro_tensor::{Tensor, MemoryOrder};
use tenferro_device::LogicalMemorySpace;
let col = MemoryOrder::ColumnMajor;
let mem = LogicalMemorySpace::MainMemory;
let mut ctx = CpuContext::new(1);
let a = Tensor::<f64>::zeros(&[3, 3], mem, col).unwrap();
let da = Tensor::<f64>::ones(&[3, 3], mem, col).unwrap();
let (result, dresult) = eigen_frule(&mut ctx, &a, &da).unwrap();