melhorias sem otimização

parent 6c80af9a
......@@ -73,8 +73,6 @@ int gaussSeidel( pde_t *pde, int maxIter, FILE *output){
dI = pde->hx * (-2 * pde->hx + pde->hx * pde->hy), //baixo
dII = pde->hy * (-2 * pde->hy + pde->hx * pde->hy); //esquerda
real dPinv = 1/dP;
// insere as bordas
for(int j=1; j < pde->ny-1 ; j++)
U[j]= sin(2*PI*(PI-pde->hx*j))*sinh(_p2(PI));
......@@ -95,7 +93,7 @@ int gaussSeidel( pde_t *pde, int maxIter, FILE *output){
for(int i = 1; i < pde->nx-1; i++){
for(int j = 1; j < pde->ny-1; j++){
b = 4 * _p2(PI) * ( sin(2 * PI * pde->hx * i) * sinh(PI * pde->hy * j) + sin(2 * PI * ( PI - pde->hx * i)) * sinh( PI * ( PI - pde->hy * j) ));
U[i*pde->ny+j] = (b - N - O - L - S ) * dPinv;
U[i*pde->ny+j] = (b - N - O - L - S ) / dP;
U[i*pde->ny+j] == INFINITY && inf();
}
}
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment