From bf3da9f200aaddfcf699048c485da1259adb7c81 Mon Sep 17 00:00:00 2001 From: Antonio De Lucreziis Date: Fri, 21 Mar 2025 13:49:34 +0000 Subject: [PATCH] fewojfpiewijfewpi cluster --- matrices/laplacian/julia-laplaciano.jl | 18 +++++++++--------- .../laplacian/laplacian-discretization-3d.mat | Bin 926216 -> 926216 bytes 2 files changed, 9 insertions(+), 9 deletions(-) diff --git a/matrices/laplacian/julia-laplaciano.jl b/matrices/laplacian/julia-laplaciano.jl index 6ca4fcb..e833cf5 100644 --- a/matrices/laplacian/julia-laplaciano.jl +++ b/matrices/laplacian/julia-laplaciano.jl @@ -8,17 +8,17 @@ nx = N ny = N nz = N -ex = fill(1, nx) -ey = fill(1, ny) -ez = fill(1, nz) +ex = fill(1, nx - 1) +ey = fill(1, ny - 1) +ez = fill(1, nz - 1) -Dxx = spdiagm(-1 => ex, 0 => -2 * ex, +1 => ex) -Dyy = spdiagm(-1 => ey, 0 => -2 * ey, +1 => ey) -Dzz = spdiagm(-1 => ez, 0 => -2 * ez, +1 => ez) +Dxx = spdiagm(-1 => ex, 0 => -2 * [ex; 1], +1 => ex) +Dyy = spdiagm(-1 => ey, 0 => -2 * [ey; 1], +1 => ey) +Dzz = spdiagm(-1 => ez, 0 => -2 * [ez; 1], +1 => ez) -Ix = spdiagm(0 => [ex; 1]) -Iy = spdiagm(0 => [ey; 1]) -Iz = spdiagm(0 => [ez; 1]) +Ix = spdiagm(0 => ex) +Iy = spdiagm(0 => ey) +Iz = spdiagm(0 => ez) L = kron(Dxx, Iy, Iz) + kron(Ix, Dyy, Iz) + kron(Ix, Iy, Dzz) diff --git a/matrices/laplacian/laplacian-discretization-3d.mat b/matrices/laplacian/laplacian-discretization-3d.mat index d33a882c1608a19dc1744134cebbe5209363748c..eac550e669ac1f0295fb0619fcce8645e1afa6da 100644 GIT binary patch delta 90 zcmeA;V$pHLVuKAQv!vnmW_!+ddrn3mW&&bnAZ7t#*6sG3Y`#4p1>5_lve`0%SncPh bu>mnV5OV-AClGT1F*gwNY(GDZ*RmS`KGh)X delta 90 zcmeA;V$pHLVuKAQ^H-M3&Gww__MD7B%ml>DK+FQftlRB5*?fCI3byx8WwT`jvD(j1 bV*_G#Am#vKP9Ww2Vs0Sj*?xW+uVps?tpFnJ