programming:transform_spectral_coeff_to_physical_data_points
Differences
This shows you the differences between two versions of the page.
| Both sides previous revision Previous revision Next revision | Previous revision Last revision Both sides next revision | ||
|
programming:transform_spectral_coeff_to_physical_data_points [2010/08/28 16:04] emmanuel |
programming:transform_spectral_coeff_to_physical_data_points [2010/08/28 16:06] emmanuel |
||
|---|---|---|---|
| Line 11: | Line 11: | ||
| We can determine the physical values from the spectral coefficients by evaluating the expansion sum (eqn 20 in the user's manual.) | We can determine the physical values from the spectral coefficients by evaluating the expansion sum (eqn 20 in the user's manual.) | ||
| - | I wanted to compared the physical values computed by this sum with the values returned by the method u.makephysical(). This method use FFTW. | + | I wanted to compared the physical values computed by this sum with the values returned by the method u.makephysical(). This method use the FFTW algorithm and |
| - | So I was just wondering if there is any difference. The fact is that, I do not obtained the same results. | + | I was expecting no difference. The fact is that, I do not obtained the same results. |
| I have attached part of the code I am using. Do you think that something is incorrect ? | I have attached part of the code I am using. Do you think that something is incorrect ? | ||
| Line 36: | Line 36: | ||
| //display gridpoint coordinates x,y,z | //display gridpoint coordinates x,y,z | ||
| std::cout << u.x(nx) << ' ' << u.y(ny) << ' ' << u.z(nz) << endl ; | std::cout << u.x(nx) << ' ' << u.y(ny) << ' ' << u.z(nz) << endl ; | ||
| + | | ||
| //compute physical value using expansion sum and display | //compute physical value using expansion sum and display | ||
| u.makeSpectral(); | u.makeSpectral(); | ||
| Line 45: | Line 45: | ||
| std::cout << 2.0 * sum_u.real()<< endl ; | std::cout << 2.0 * sum_u.real()<< endl ; | ||
| - | //display physical value from makePhysical() | + | //display physical value from makePhysical() |
| u.makePhysical(); | u.makePhysical(); | ||
| std::cout << u(nx,ny,nz,0) << endl ; | std::cout << u(nx,ny,nz,0) << endl ; | ||
| - | + | | |
| std::cin >> pause; | std::cin >> pause; | ||
| } | } | ||
programming/transform_spectral_coeff_to_physical_data_points.txt · Last modified: 2010/08/28 16:08 by emmanuel
Page Tools
Except where otherwise noted, content on this wiki is licensed under the following license: CC Attribution-Share Alike 3.0 Unported
