Changeset 4969


Ignore:
Timestamp:
Jan 23, 2008, 8:34:11 PM (16 years ago)
Author:
steve
Message:

Fix up error in advection_ext.c

Location:
anuga_core/source/anuga/advection
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • anuga_core/source/anuga/advection/advection.py

    r4967 r4969  
    259259
    260260
    261         print 'N = ',N
    262         print 'timestep = ',timestep
    263         print 'huge_timestep = ',huge_timestep
     261        #print 'N = ',N
     262        #print 'timestep = ',timestep
     263        #print 'huge_timestep = ',huge_timestep
    264264               
    265265        timestep = compute_fluxes(stage_edge,stage_bdry,stage_update,
     
    269269                                  huge_timestep,max_timestep,v,N)
    270270       
    271         print 'timestep out2 =',timestep
     271        #print 'timestep out2 =',timestep
    272272
    273273        self.timestep = timestep
  • anuga_core/source/anuga/advection/advection_ext.c

    r4967 r4969  
    2525        double optimal_timestep;
    2626        double timestep;
    27         int I,M,K;
     27        int k_i,n_m,k_i_j;
    2828        int k,i,j,n,m;
    2929       
    3030        timestep = max_timestep;
    3131
    32         printf("N = %i\n",N);
    33         printf("timestep = %g\n",timestep);
    34         printf("huge_timestep = %g\n",huge_timestep);
     32        //printf("N = %i\n",N);
     33        //printf("timestep = %g\n",timestep);
     34        //printf("huge_timestep = %g\n",huge_timestep);
    3535
    36                
    3736        for (k=0; k<N; k++){
    3837            optimal_timestep = huge_timestep;
    3938            flux = 0.0;  //Reset work array
    4039            for (i=0; i<3; i++){
    41                                 I = 6*k+i;
     40                k_i = 3*k+i;
    4241                //Quantities inside volume facing neighbour i
    43                 ql = stage_edge[I];
     42                ql = stage_edge[k_i];
     43                //printf("stage_edge[%i] = %g\n",k_i,stage_edge[k_i]);
    4444
    4545                //Quantities at neighbour on nearest face
    46                 n = neighbours[I];
     46                n = neighbours[k_i];
    4747                if (n < 0) {
    4848                    m = -n-1; //Convert neg flag to index
    4949                    qr = stage_bdry[m];
    5050                } else {
    51                     m = neighbour_edges[I];
    52                                         M = 6*n+m;
    53                     qr = stage_edge[M];
     51                    m = neighbour_edges[k_i];
     52                    n_m = 3*n+m;
     53                    qr = stage_edge[n_m];
    5454                }
    5555
     
    5757                //Outward pointing normal vector
    5858                for (j=0; j<2; j++){
    59                                         K = 6*k+2*i+j;
    60                     normal[j] = normals[K];
     59                    k_i_j = 6*k+2*i+j;
     60                    normal[j] = normals[k_i_j];
    6161                }
    6262
     
    7272
    7373                max_speed = fabs(normal_velocity);
    74                 flux = flux - edgeflux * edgelengths[K];
     74                flux = flux - edgeflux * edgelengths[k_i];
    7575
    7676                //Update optimal_timestep
     
    9191        }
    9292
    93         printf("timestep out = %g\n",timestep);
     93        //printf("timestep out = %g\n",timestep);
    9494
    9595        return timestep;
Note: See TracChangeset for help on using the changeset viewer.