Nils Werner's questions

Escrevi uma classe de filtro IIR com um estado de filtro interno e uma função elemento a elemento float filter(float x). Naturalmente, devido ao estado de filtro interno, ela filter()só pode ser chamada nas amostras na ordem em que estão organizadas na sequência.

Se eu agora aplicar essa função a um Eigen Array, como

input.unaryExpr(filter);

a ordem de execução será unaryExpr()sempre estritamente na ordem dos valores no array, ou pode ocorrer execução fora de ordem ou até mesmo paralelização?

Seria mais seguro escrever o loop explicitamente para garantir que a ordem seja sempre a esperada?

No momento, parece estar funcionando corretamente, mas não consigo encontrar nenhuma documentação explícita sobre seu comportamento.

Estou tentando comparar dois ndarray::Array1<Complex64>usando a approxcaixa. Posso comparar duas f64matrizes sem problemas:

#[test]
fn repro() {
    let a = ndarray::array![1.0, 0.0];
    let b = ndarray::array![1.0, 0.0];

    assert!(a.abs_diff_eq(&b, 0.1));
}

mas, durante toda a minha vida, não posso comparar duas matrizes complexas:

#[test]
fn repro() {
    let a = ndarray::array![
        Complex64::new(1.0, 0.0)
    ];
    let b = ndarray::array![
        Complex64::new(1.0, 0.0)
    ];

    assert!(a.abs_diff_eq(&b, 0.1));
}

com a mensagem de erro

  --> src/lib.rs:85:31
   |
85 |         assert!(a.abs_diff_eq(&b, 0.1));
   |                   ----------- ^^^^^^ the trait `approx::abs_diff_eq::AbsDiffEq<_>` is not implemented for `Complex<f64>`
   |                   |
   |                   required by a bound introduced by this call

Mas eu sei que essa característica está de fato implementada paraComplex<T> e o recurso num-complexestá habilitado.

o que estou perdendo?