studiorailgun/Renderer
2
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

remove prints from multigrid
Some checks failed
studiorailgun/Renderer/pipeline/head There was a failure building this commit
Details

Browse Source
This commit is contained in:
austin 2024-12-11 18:17:56 -05:00
parent 1d75ad9af9
commit 96e03bea33
1 changed files with 0 additions and 19 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

19
src/main/c/src/math/ode/multigrid.c
View File

@ -66,8 +66,6 @@ void initialization_check();
float solver_multigrid_iterate_serial_recursive(float * phi, float * phi0, float a, float c, int GRIDDIM){
initialization_check();
float residual;
printf("Current dim: %d \n",GRIDDIM);
fflush(stdout);
int LOWERDIM = ((GRIDDIM - 2) / 2) + 2;
float * currResidual = get_current_residual(GRIDDIM);
@ -75,11 +73,6 @@ float solver_multigrid_iterate_serial_recursive(float * phi, float * phi0, float
float * lowerPhi0 = get_current_phi0(LOWERDIM);
float * lowerResidual = get_current_residual(LOWERDIM);
printf("Solved arrays: %p %p %p %p \n",currResidual,lowerPhi, lowerPhi0, lowerResidual);
fflush(stdout);
printf("smooth\n");
fflush(stdout);
//
//gauss-seidel at highest res
solver_gauss_seidel_iterate_parallel(phi,phi0,a,c,GRIDDIM);
@ -89,8 +82,6 @@ float solver_multigrid_iterate_serial_recursive(float * phi, float * phi0, float
//compute residuals
residual = solver_multigrid_store_residual_serial(phi,phi0,currResidual,a,c,GRIDDIM);
printf("restrict\n");
fflush(stdout);
//restrict
//(current operator is injection -- inject r^2 from this grid at phi0 of the smaller grid)
for(int x = 0; x < LOWERDIM; x++){
@ -111,8 +102,6 @@ float solver_multigrid_iterate_serial_recursive(float * phi, float * phi0, float
}
}
printf("recurse\n");
fflush(stdout);
//solve next-coarsest grid
if(GRIDDIM <= LOWEST_DIM){
//smooth
@ -120,8 +109,6 @@ float solver_multigrid_iterate_serial_recursive(float * phi, float * phi0, float
//compute residual
residual = solver_multigrid_store_residual_serial(quarterGridPhi,quarterGridPhi0,quarterGridResidual,a,c,quarterDim);
} else {
printf("solve\n");
fflush(stdout);
float solution =
(
phi0[solver_gauss_seidel_get_index(1,1,1,GRIDDIM)] +
@ -135,8 +122,6 @@ float solver_multigrid_iterate_serial_recursive(float * phi, float * phi0, float
) / 8.0f
;
printf("interpolate solution\n");
fflush(stdout);
//interpolate from the lower grid
for(int x = 1; x < GRIDDIM - 1; x++){
for(int y = 1; y < GRIDDIM - 1; y++){
@ -148,8 +133,6 @@ float solver_multigrid_iterate_serial_recursive(float * phi, float * phi0, float
}
}
printf("interpolate\n");
fflush(stdout);
//interpolate from the lower grid
for(int x = 1; x < GRIDDIM - 1; x++){
for(int y = 1; y < GRIDDIM - 1; y++){
@ -183,8 +166,6 @@ float solver_multigrid_iterate_serial_recursive(float * phi, float * phi0, float
printf("smooth\n");
fflush(stdout);
//
// full res
//