Merge pull request #1462 from xtalax/pdesystem

Vectorize pdeeqs

Author: ChrisRackauckas

docs/src/index.md

@@ -136,7 +136,7 @@ system:
     - Multi-package interface for numerical solving `OptimizationSystem`
 - [NeuralPDE.jl](https://github.com/SciML/NeuralPDE.jl)
     - Physics-Informed Neural Network (PINN) training on `PDESystem`
-- [DiffEqOperators.jl](https://github.com/SciML/DiffEqOperators.jl)
+- [MethodOfLines.jl](https://github.com/SciML/MethodOfLines.jl)
     - Automated finite difference method (FDM) discretization of `PDESystem`
 
 ## Contributing

src/systems/pde/pdesystem.jl

@@ -70,6 +70,7 @@ struct PDESystem <: ModelingToolkit.AbstractMultivariateSystem
         if checks
             all_dimensionless([dvs;ivs;ps]) ||check_units(eqs)
         end
+        eqs = eqs isa Vector ? eqs : [eqs]
         new(eqs, bcs, domain, ivs, dvs, ps, defaults, connector_type, name)
     end
 end