Os métodos Newton-Krylov livre de Jacobiano (JFNK) e métodos de Krylov em geral podem ser muito úteis porque não requerem armazenamento explícito ou construção de uma matriz, apenas os resultados de produtos vetoriais de matriz. Se você realmente formar o sistema esparso, existem muitos pré-condicionadores para você.
O que está disponível para métodos verdadeiros sem matriz? O Google busca algumas referências à "estimativa de matriz" e outras coisas indicando que isso é possível. Como esses métodos geralmente funcionam? Como eles se comparam aos pré-condicionadores tradicionais? Os pré-condicionadores livres de matriz baseados em física são o caminho a percorrer? Existem métodos disponíveis na natureza, digamos no PETSc ou em algum outro pacote?