Metode Jacobian-Free Newton-Krylov (JFNK), dan metode Krylov secara umum, dapat sangat berguna karena mereka tidak memerlukan penyimpanan eksplisit atau konstruksi matriks, hanya hasil dari produk vektor-matriks. Jika Anda benar-benar membentuk sistem sparse, ada banyak prasyarat di luar sana untuk Anda.
Apa yang tersedia untuk metode bebas matriks yang benar? Googling menghasilkan beberapa referensi untuk "estimasi matriks" dan beberapa hal lain yang menunjukkan bahwa itu mungkin. Bagaimana cara kerja metode ini secara umum? Bagaimana mereka dibandingkan dengan prekondisi tradisional? Apakah prasyarat berbasis matriks bebas fisika adalah jalan yang harus ditempuh? Apakah ada metode yang tersedia secara terbuka di alam liar, katakanlah dalam PETSc atau paket lain?