Dear Ales,
Not hearing back from you, I have performed a few tests myslef, in 1D, 2D and 3D. Please the results quoted below. I see no problems with the code, and close this bug case.
Thanks, Andrew
+++++++++++++++++++++++++++++++++++
octave:52> [lambda,V,A] = laplacian([2,4,5],{'P' 'P' 'P'}, 20);
ans =
Warning: (m+1)th eigenvalue is nearly equal
to mth.
The eigenvectors take 6400 bytes
The Laplacian matrix takes 3524 bytes
octave:53> whos
Variables in the current scope:
Attr Name Size Bytes Class
==== ==== ==== ===== =====
A 40x40 3524 double
V 40x20 6400 double
lambda 20x1 160 double
Total is 1100 elements using 10084 bytes
octave:54> [lambda,V,A] = laplacian([2,4],{'P' 'P' }, 2);
ans =
Warning: (m+1)th eigenvalue is nearly equal
to mth.
The eigenvectors take 128 bytes
The Laplacian matrix takes 516 bytes
octave:55> whos
Variables in the current scope:
Attr Name Size Bytes Class
==== ==== ==== ===== =====
A 8x8 516 double
V 8x2 128 double
lambda 2x1 16 double
Total is 58 elements using 660 bytes
octave:56> [lambda,V,A] = laplacian([2],{'P'}, 2);
The eigenvectors take 32 bytes
The Laplacian matrix takes 60 bytes
octave:57> whos
Variables in the current scope:
Attr Name Size Bytes Class
==== ==== ==== ===== =====
A 2x2 60 double
V 2x2 32 double
lambda 2x1 16 double
Total is 10 elements using 108 bytes
Re: no solution when periodic BC are applied.
Dear Ales,
Could you please provide a specific complete example, quoting the function call, the full output, and MATLAB version? I cannot fix a problem without being able to reproduce the problem myslef.
Thanks, Andrew
Dear Andrew:
Thank you for a nice software. I ran into a problem -- there is no solution when periodic BC are applied. I think it is endemic to these types of problems. Do you know why?
Comment only