User Tools

Site Tools


programming:transform_spectral_coeff_to_physical_data_points

Differences

This shows you the differences between two versions of the page.

Link to this comparison view

Both sides previous revision Previous revision
Next revision
Previous revision
programming:transform_spectral_coeff_to_physical_data_points [2010/08/28 10:48]
emmanuel
programming:transform_spectral_coeff_to_physical_data_points [2010/08/28 16:08]
emmanuel
Line 10: Line 10:
 Dear Prof. Gibson, Dear Prof. Gibson,
  
-I would like to transform ​the spectral coefficients in physical values 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.)  
-Starting from an object Flowfield u, have compared the results of the expansion ​with the value given by the method u.makephysical(). I do not obtain ​the same results. ​+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  
 +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 17: Line 19:
  
 Example of results displayed : Example of results displayed :
- 0.00258953 ​     // u computing Fourier series +  ​0.00258953 ​     // u from expansion sum 
- ​0.00604268 ​     // return of u.makePhysical+  0.00604268 ​     // return of u.makePhysical ​(FFTW)
  
  
Line 34: 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 43: 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