Accelerating Data Generation for Neural Operators via Krylov Subspace Recycling