Методы Ньютона-Крылова без Якобиана (JFNK) и методы Крылова в целом могут быть очень полезными, поскольку они не требуют явного хранения или построения матрицы, а только результатов произведений матрицы-вектора. Если вы действительно формируете разреженную систему, для вас есть много предварительных условий.
Что доступно для истинных безматричных методов? Поиск в Google приводит некоторые ссылки на «матричную оценку» и некоторые другие вещи, указывающие, что это возможно. Как эти методы обычно работают? Как они сравниваются с традиционными предварительными кондиционерами? Являются ли основанные на физике предварительные кондиционеры без матрицы подходящими? Есть ли какие-либо общедоступные методы в дикой природе, скажем, в PETSc или каком-либо другом пакете?