Guest User

Untitled

a guest
May 20th, 2018
148
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
C 1.01 KB | None | 0 0
  1. void linha_pivot_primal(Matriz *mtx){
  2.     int check = 0;
  3.     int y;
  4.     double comp = INT_MAX, elem = 0;
  5.    
  6.     for(y = 0; y < (mtx->j)-1; y++){
  7.         if(mtx->matriz[y][mtx->coluna_pivot] > 0 && mtx->matriz[y][mtx->coluna_pivot] != 0
  8.         && abs(mtx->matriz[y][(mtx->i)+(mtx->j)-2]/mtx->matriz[y][mtx->coluna_pivot]) < comp ){
  9.             comp = abs(mtx->matriz[y][(mtx->i)+(mtx->j)-2]/mtx->matriz[y][mtx->coluna_pivot]);
  10.             mtx->linha_pivot = y;
  11.             check = 1;
  12.             elem = mtx->matriz[y][mtx->coluna_pivot];
  13.         }
  14.         if(mtx->matriz[y][mtx->coluna_pivot] > 0 && mtx->matriz[y][mtx->coluna_pivot] != 0
  15.         && abs(mtx->matriz[y][(mtx->i)+(mtx->j)-2]/mtx->matriz[y][mtx->coluna_pivot]) == comp ){
  16.             if(mtx->matriz[y][mtx->coluna_pivot] > elem){
  17.                 comp = abs(mtx->matriz[y][(mtx->i)+(mtx->j)-2]/mtx->matriz[y][mtx->coluna_pivot]);
  18.                 mtx->linha_pivot = y;
  19.                 elem = mtx->matriz[y][mtx->coluna_pivot];
  20.                 mtx->degenerada_flag = 1;
  21.             }              
  22.         }
  23.     }
  24.     if (check == 0) // n ha elementos positivos na coluna pivot
  25.     else // encontro um elem_pivot
  26. }
Add Comment
Please, Sign In to add comment