rcce/lu/jacu.c

14   double  tmp1, tmp2, tmp3;
26                tmp1 = 1.0 / u(1,i,j,k);
27                tmp2 = tmp1 * tmp1;
28                tmp3 = tmp1 * tmp2;
30                tmp1 = 1.0 / u(1,i,j,k);
31                tmp2 = tmp1 * tmp1;
32                tmp3 = tmp1 * tmp2;
49                 * (  tx1 * r43 * c34 * tmp1
50                    + ty1 *       c34 * tmp1
51                    + tz1 *       c34 * tmp1 )
66                     * (  tx1 *       c34 * tmp1
67                        + ty1 * r43 * c34 * tmp1
68                        + tz1 *       c34 * tmp1 )
83                     * (  tx1 *       c34 * tmp1
84                        + ty1 *       c34 * tmp1
85                        + tz1 * r43 * c34 * tmp1 )
117          + dt * 2.0 * ( tx1 * c1345 * tmp1
118                           + ty1 * c1345 * tmp1
119                           + tz1 * c1345 * tmp1 )
127                tmp1 = 1.0 / u(1,i+1,j,k);
128                tmp2 = tmp1 * tmp1;
129                tmp3 = tmp1 * tmp2;
138                 * ( - ( u(2,i+1,j,k) * tmp1 ) *( u(2,i+1,j,k) * tmp1 )
144                 * ( ( 2.0 - c2 ) * ( u(2,i+1,j,k) * tmp1 ) )
145                 - dt * tx1 * ( r43 * c34 * tmp1 )
148                     * ( - c2 * ( u(3,i+1,j,k) * tmp1 ) );
150                     * ( - c2 * ( u(4,i+1,j,k) * tmp1 ) );
156                a(3,2,i,j) =  dt * tx2 * ( u(3,i+1,j,k) * tmp1 );
157                a(3,3,i,j) =  dt * tx2 * ( u(2,i+1,j,k) * tmp1 )
158                 - dt * tx1 * ( c34 * tmp1 )
166                a(4,2,i,j) = dt * tx2 * ( u(4,i+1,j,k) * tmp1 );
168                a(4,4,i,j) = dt * tx2 * ( u(2,i+1,j,k) * tmp1 )
169                 - dt * tx1 * ( c34 * tmp1 )
177                     - c1 * ( u(5,i+1,j,k) * tmp1 ) )
178                 * ( u(2,i+1,j,k) * tmp1 ) )
185                 * ( c1 * ( u(5,i+1,j,k) * tmp1 )
201                  * ( c1 * ( u(2,i+1,j,k) * tmp1 ) )
202                  - dt * tx1 * c1345 * tmp1
208                tmp1 = 1.0 / u(1,i,j+1,k);
209                tmp2 = tmp1 * tmp1;
210                tmp3 = tmp1 * tmp2;
221                b(2,2,i,j) =  dt * ty2 * ( u(3,i,j+1,k) * tmp1 )
222                 - dt * ty1 * ( c34 * tmp1 )
224                b(2,3,i,j) =  dt * ty2 * ( u(2,i,j+1,k) * tmp1 );
229                  * ( - ( u(3,i,j+1,k) * tmp1 ) * ( u(3,i,j+1,k) * tmp1 )
236                          * ( - c2 * ( u(2,i,j+1,k) * tmp1 ) );
238                          * ( u(3,i,j+1,k) * tmp1 ) )
239              - dt * ty1 * ( r43 * c34 * tmp1 )
242                          * ( - c2 * ( u(4,i,j+1,k) * tmp1 ) );
249                b(4,3,i,j) =  dt * ty2 * ( u(4,i,j+1,k) * tmp1 );
250                b(4,4,i,j) =  dt * ty2 * ( u(3,i,j+1,k) * tmp1 )
251                               - dt * ty1 * ( c34 * tmp1 )
259                      - c1 * ( u(5,i,j+1,k) * tmp1 ) )
260                 * ( u(3,i,j+1,k) * tmp1 ) )
271                 * ( c1 * ( u(5,i,j+1,k) * tmp1 )
282                 * ( c1 * ( u(3,i,j+1,k) * tmp1 ) )
283                 - dt * ty1 * c1345 * tmp1
289                tmp1 = 1.0 / u(1,i,j,k+1);
290                tmp2 = tmp1 * tmp1;
291                tmp3 = tmp1 * tmp2;
302                c(2,2,i,j) = dt * tz2 * ( u(4,i,j,k+1) * tmp1 )
303                  - dt * tz1 * c34 * tmp1
306                c(2,4,i,j) = dt * tz2 * ( u(2,i,j,k+1) * tmp1 );
313                c(3,3,i,j) = dt * tz2 * ( u(4,i,j,k+1) * tmp1 )
314                  - dt * tz1 * ( c34 * tmp1 )
316                c(3,4,i,j) = dt * tz2 * ( u(3,i,j,k+1) * tmp1 );
320               * ( - ( u(4,i,j,k+1) * tmp1 ) *( u(4,i,j,k+1) * tmp1 )
327                    * ( - c2 * ( u(2,i,j,k+1) * tmp1 ) );
329                    * ( - c2 * ( u(3,i,j,k+1) * tmp1 ) );
331                    * ( u(4,i,j,k+1) * tmp1 )
332                    - dt * tz1 * ( r43 * c34 * tmp1 )
340              - c1 * ( u(5,i,j,k+1) * tmp1 ) )
341                   * ( u(4,i,j,k+1) * tmp1 ) )
354              * ( c1 * ( u(5,i,j,k+1) * tmp1 )
361              * ( c1 * ( u(4,i,j,k+1) * tmp1 ) )
362              - dt * tz1 * c1345 * tmp1