Hi everyone

I’m studying a pde (swift-hohenberg model) converting it in the ODE form.

I’ve discretized the laplacian:

```
function Laplacian2D(Nx, lx)
hx = lx/Nx
D2x = CenteredDifference(2, 2, hx, Nx)
Qx = PeriodicBC(hx |> typeof)
A = kron(sparse(I, Nx, Nx), sparse(D2x * Qx)[1]) + kron(sparse(D2x * Qx)[1], sparse(I, Nx, Nx))
return A
end
```

After that I have built the function for the ODEProblem:

```
function F_sh2!(du,u,p,t)
@unpack l , d, d2 = p
du .= - u .- d*u .- d2*u .+ (l .* u .- u.^3)
return nothing
end
```

Where `l`

is a float, `d`

is the laplacian and `d2`

is laplacian^2.

For no large discretizations (64x64 grid) it’s all ok, but if I try to increase the grid, for instance 128x128, the solver returns error: `ReadOnlyMemoryError()`

I’m wondering if it’s only due to my computational limits or if I’m missing some errors.

Every help will be welcome, thank you so much!