Abstract:
The present invention relates to a parallel device for solving linear equations over finite fields, including a processor, an input port, an output port, a pivot finding component, a partial inversion component, a normalization component and an elimination component. The processor is connected to each of the pivot finding component, the partial inversion component, the normalization component, the elimination component, and the input port and the output port. The partial inversion component is connected to the elimination component and the normalization component. The pivot finding component is connected to the elimination component. The present invention enables parallel computing to a certain extent with fast solving speed and simple design, and thus can be widely used in various engineering fields.

Description:
FIELD OF THE INVENTION 
       [0001]    The present invention relates generally to a device for solving linear equations, and more particularly to a parallel device for solving linear equations over finite fields. 
       BACKGROUND OF THE INVENTION 
       [0002]    Finite fields are fields containing only a finite number of elements, and linear equations over finite fields are equations of which each factor is an element in finite fields. Finding solutions to linear equations over finite fields is widely used in various engineering fields, such as the field of cryptography, and also the field of solving other mathematical problems. 
         [0003]    Methods adopted in solving linear equations include Gaussian elimination method and Gauss-Jordan method. By Gaussian elimination method, linear equations are multi-iterated into upper or lower triangular forms, wherein each iteration operation includes three sub-operations: finding pivot, normalization and elimination. If the equations are solvable, then the final solutions are obtained by substitution operation. Gauss-Jordan method, a variant of Gaussian elimination method, is able to solve linear equations by multiple iterations, but it consumes more resources than Gaussian elimination method. 
         [0004]    Solving linear equations is a highly computationally complex and time-consuming issue. At present, there is still large room for optimization of solving linear equations, especially for optimization of those over finite fields. So far, devices dedicated to solving linear equations over finite fields have not yet been reported. 
       SUMMARY OF THE INVENTION 
       [0005]    Therefore, in order to address the deficiencies and inadequacies in the art, the present invention aims to provide a parallel device for solving linear equations over finite fields. 
         [0006]    The object of the invention is achieved by the following technical solutions. 
         [0007]    A parallel device for solving linear equations over finite fields, including: 
         [0008]    an input port for inputting coefficient matrix B of linear equations over finite fields and irreducible polynomial p(x) selected over the field GF(2 n ); the coefficient matrix B is an m×(m+1) matrix; the element a(x) of the coefficient matrix B is an element in the field GF(2 n ); 
         [0009]    a first processor, including a scheduler and a memory interconnected with each other; the scheduler is configured to control a pivot finding component, a partial inversion component, a normalization component and an elimination component; the memory is configured to store the inputted coefficient matrix B, and to update the stored coefficient matrix B after each iteration operation; 
         [0010]    a pivot finding component, including a second processor for finding out the pivot β of the coefficient matrix B; 
         [0011]    a partial inversion component, including a third processor for implementing partial inversion calculation; 
         [0012]    a normalization component, including a fourth processor and m+1 normalization calculation units; the fourth processor is connected to each of the m+1 normalization calculation units to perform scheduling of the normalization calculation units; 
         [0013]    an elimination component, including a fifth processor and m×(m+1) elimination calculation units; the fifth processor is connected to each of the m×(m+1) elimination calculation units to perform scheduling of the elimination calculation units and data transmission; 
         [0014]    an output port for outputting the results of linear equations over finite fields being solved, 
         [0015]    wherein the first processor is connected to the pivot finding component, the partial inversion component, the normalization component, the elimination component, and the input port and the output port; the partial inversion component is connected to the elimination component and the normalization component; 
         [0016]    wherein the first processor receives the coefficient matrix B and the irreducible polynomial p(x) which is selected over the field GF(2 n ) and outputted from the input port, stores the coefficient matrix B into the memory, sets the counter at m, and sends the column containing the pivot to be found in present iteration process to the pivot finding component; the column containing the pivot to be found in present iteration process is set so that in the g th  iteration process, the column containing the pivot to be found is the g th  column, 0≦g≦m; 
         [0017]    wherein the pivot finding component implements a process of finding the pivot: determining whether an element with a row number of g and a column number of g is a none-zero element; if yes, this element is determined to be a pivot, and the pivot finding component sends no back feed into the first processor; 
         [0018]    if not, then finding none-zero elements one by one among the elements with row numbers of g+1 to m and a column number of g; the first none-zero element found is determined to be the pivot β, and the pivot finding component sends the row number of the pivot as feedback to the processor; 
         [0019]    the first processor sends the pivot β found to the partial inversion component, and sends the row where the pivot β is in to the normalization component, and sends the other lines of the input coefficient matrix B to the elimination component; 
         [0020]    the partial inversion component implements partial inversion calculation, and outputs the calculated results to the normalization component and the elimination component; 
         [0021]    the normalization component and the elimination component implement the normalization calculation and the elimination calculation respectively, and output the calculated results to the first processor; 
         [0022]    the first processor updates the coefficient matrix B according to the results of the normalization calculation and the elimination calculation, stores the updated coefficient matrix B into the memory, subtracts 1 from the value of the counter, and sets the present counter as j; a next iteration is implemented until the counter reaches 0; 
         [0023]    if the counter is 0, then outputting the last column of the updated coefficient matrix B into the port as the solved results. 
         [0024]    The normalization calculation units are logical gate circuits. 
         [0025]    The elimination calculation units are logical gate circuits. 
         [0026]    The partial inversion calculation specifically includes: 
         [0027]    for i=1, . . . , n−1, the third processor calculates β 2     i   ; 
         [0028]    let i′=(n−1)÷3, for k=0, 1, . . . , i′−1, calculating S k =MUL3(β 2     3k+1   , β 2     3k+2   , β 2     3k+3   ); MUL3 is the multiplication of three operands defined over GF(2 n ); 
         [0029]    outputting β 2     i    and S k  respectively to the normalization component and the elimination component, and h=3i′+1, 3i′+2 . . . n−1. 
         [0030]    The normalization calculation specifically includes: 
         [0031]    the l th  normalization calculation unit calculates 
         [0000]    
       
         
           
             
               
                 S 
                 
                   i 
                   ′ 
                 
               
               = 
               
                 
                   a 
                   tl 
                 
                 × 
                 
                   
                     ∏ 
                     
                       k 
                       = 
                       
                         
                           3 
                            
                           
                             i 
                             ′ 
                           
                         
                         + 
                         1 
                       
                     
                     
                       n 
                       - 
                       1 
                     
                   
                    
                   
                     β 
                     
                       2 
                       k 
                     
                   
                 
               
             
             , 
           
         
       
     
         [0000]    wherein i=(n−1)÷3, and a t1  is the element in row t and column l of the coefficient matrix B; t is the row where the pivot of the present iteration is in; then calculating 
         [0000]    
       
         
           
             
               
                 a 
                 tl 
               
               = 
               
                 
                   ∏ 
                   
                     q 
                     = 
                     0 
                   
                   
                     i 
                     ′ 
                   
                 
                  
                 
                   S 
                   q 
                 
               
             
             , 
           
         
       
     
         [0000]    wherein l=0, 1, 2, . . . , m. 
         [0032]    The elimination calculation specifically includes: 
         [0033]    the elimination calculation unit numbered (k′, l) calculates 
         [0000]    
       
         
           
             
               
                 S 
                 
                   
                     i 
                     ′ 
                   
                   + 
                   1 
                 
               
               = 
               
                 { 
                 
                   
                     
                       
                         
                           
                             a 
                             
                               
                                 k 
                                 ′ 
                               
                                
                               t 
                             
                           
                           × 
                           
                             a 
                             tl 
                           
                           × 
                           
                             
                               ∏ 
                               
                                 l 
                                 = 
                                 
                                   
                                     3 
                                      
                                     
                                       i 
                                       ′ 
                                     
                                   
                                   + 
                                   1 
                                 
                               
                               
                                 n 
                                 - 
                                 1 
                               
                             
                              
                             
                               β 
                               
                                 2 
                                 l 
                               
                             
                           
                         
                         , 
                         
                           
                             j 
                             ′ 
                           
                           &lt; 
                           2 
                         
                       
                     
                   
                   
                     
                       
                         
                           
                             α 
                             tl 
                           
                           × 
                           
                             
                               ∏ 
                               
                                 l 
                                 = 
                                 
                                   
                                     3 
                                      
                                     
                                       i 
                                       ′ 
                                     
                                   
                                   + 
                                   1 
                                 
                               
                               
                                 n 
                                 - 
                                 1 
                               
                             
                              
                             
                               β 
                               
                                 2 
                                 l 
                               
                             
                           
                         
                         , 
                         
                           
                             j 
                             ′ 
                           
                           = 
                           2 
                         
                       
                     
                   
                 
                 } 
               
             
             , 
           
         
       
     
         [0000]    wherein i′=(n−1)÷3, j′=(n−1)mod 3, and mod represents modular operation; 
         [0034]    then calculating 
         [0000]    
       
         
           
             
               
                 a 
                 
                   
                     k 
                     ′ 
                   
                    
                   l 
                 
               
               = 
               
                 { 
                 
                   
                     
                       
                         
                           
                             
                               S 
                               
                                 
                                   i 
                                   ′ 
                                 
                                 + 
                                 1 
                               
                             
                             × 
                             
                               
                                 ∏ 
                                 
                                   q 
                                   = 
                                   0 
                                 
                                 
                                   
                                     i 
                                     ′ 
                                   
                                   - 
                                   1 
                                 
                               
                                
                               
                                 S 
                                 q 
                               
                             
                           
                           + 
                           
                             a 
                             
                               
                                 k 
                                 ′ 
                               
                                
                               l 
                             
                           
                         
                         , 
                         
                           
                             j 
                             ′ 
                           
                           &lt; 
                           2 
                         
                       
                     
                   
                   
                     
                       
                         
                           
                             
                               a 
                               
                                 
                                   k 
                                   ′ 
                                 
                                  
                                 t 
                               
                             
                             × 
                             
                               S 
                               
                                 
                                   i 
                                   ′ 
                                 
                                 + 
                                 1 
                               
                             
                             × 
                             
                               
                                 ∏ 
                                 
                                   q 
                                   = 
                                   0 
                                 
                                 
                                   
                                     i 
                                     ′ 
                                   
                                   - 
                                   1 
                                 
                               
                                
                               
                                 S 
                                 q 
                               
                             
                           
                           + 
                           
                             a 
                             
                               
                                 k 
                                 ′ 
                               
                                
                               l 
                             
                           
                         
                         , 
                         
                           
                             j 
                             ′ 
                           
                           = 
                           2 
                         
                       
                     
                   
                 
                 } 
               
             
             , 
           
         
       
     
         [0000]    wherein k′=0, 1, 2, . . . , m−1; l=0, 1, 2, . . . , m. 
         [0035]    The irreducible polynomial p(x) selected over the field GF(2 n ) has the following form: 
         [0000]        p ( x )= x   n   +p   n−1   x   n−1   +p   n−2   x   n−2   + . . . +p   1   x+ 1. 
         [0036]    The element a(x) of the coefficient matrix B has the following form: 
         [0000]        a ( x )= a   n 1   x   n 1   +a   n 2   x   n 2   + . . . +a   0 . 
         [0037]    Compared with the prior art, the present invention has the following advantages and technical effects. 
         [0038]    By setting independent pivot finding component, partial inversion component, normalization component and elimination component, the present invention realizes parallel computing to a certain extent. The parallel device for solving linear equations over finite fields of the present invention is fast in solving, and simple in design; the partial inversion component, normalization component and elimination component, as well as the solving device for linear equations and other computing devices over finite fields can be widely used in various engineering fields, especially in the hardware implementation of cryptographic algorithms and in solving a variety of mathematical problems. 
     
    
     
       BRIEF DESCRIPTION OF THE DRAWINGS 
         [0039]      FIG. 1  is a schematic structure chart of the parallel device for solving linear equations over finite fields in one embodiment of the present invention; 
           [0040]      FIG. 2  is a schematic structure chart of the processor in one embodiment of the present invention; 
           [0041]      FIG. 3  is a schematic structure chart of the pivot finding component in one embodiment of the present invention; 
           [0042]      FIG. 4  is a schematic structure chart of the partial inversion component in one embodiment of the present invention; 
           [0043]      FIG. 5  is a schematic structure chart of the normalization component in one embodiment of the present invention; and 
           [0044]      FIG. 6  is a schematic structure chart of the elimination component in one embodiment of the present invention. 
       
    
    
     DETAILED DESCRIPTION OF THE INVENTION 
       [0045]    The invention will be better understood with reference to the following description taken in conjunction with the specific embodiments and the accompanying drawings. Accordingly, the scope of the present invention is defined by the appended claims rather than the foregoing description and the exemplary embodiments described therein. 
       EXAMPLES 
       [0046]    As illustrated in  FIG. 1 , the parallel device for solving linear equations over finite fields of the embodiment includes input ports, a first processor, a pivot finding component, a partial inversion component and an output port; the first processor is connected respectively to the pivot finding component, the partial inversion component, the normalization component, the elimination component, the input port and the output port; the partial inversion component is connected respectively to the elimination component and the normalization component. 
         [0047]    The following is a detailed description of the components of the parallel device of the embodiment. 
         [0048]    (1) The input ports: as illustrated in  FIG. 1 , there are two input ports in the parallel device of the embodiment of the invention, wherein input port P is configured to input the irreducible polynomial p(x) selected over the field GF(2 n ); input port A is configured to input the coefficient matrix B of the linear equations to be solved. 
         [0049]    The coefficient matrix B is a m×(m+1) matrix, and a(x) is an element of the coefficient matrix B; a(x) and p(x) can be in the form of: 
         [0000]        a ( x )= a   n 1   x   n 1   +a   n 2   x   n 2   + . . . +a   0 ; 
         [0000]        p ( x )= x   n   +p   n−1   x   n−1   +p   n−2   x   n−2   + . . . +p   1   x+ 1; 
         [0000]    wherein a n−1 , a n−2 , . . . , a 0  and p n−1 , p n−2 , . . . , p 1  are elements in GF(2). 
         [0050]    (2) The first processor: as the only component that can communicate with I/O port, the first processor is a core component of the device of the present invention, which can control the pivot finding component, the partial inversion component, the normalization component and the elimination component. 
         [0051]    As illustrated in  FIG. 2 , the first processor includes a scheduler and a memory interconnected with each other; the scheduler is configured to control the pivot finding component, the partial inversion component, the normalization component and the elimination component; the memory is configured to store the inputted coefficient matrix B, and to update the stored coefficient matrix B after each iteration operation. 
         [0052]    (3) The pivot finding component, as illustrated in  FIG. 3 , includes a second processor for finding out the pivot β of the coefficient matrix B. 
         [0053]    (4) The partial inversion component, as illustrated in  FIG. 4 , includes a third processor for implementing partial inversion calculation of the pivot. 
         [0054]    (5) The normalization component, as illustrated in  FIG. 5 , includes a fourth processor and m+1 normalization calculation units d 0 , d 2  . . . d m ; the fourth processor is connected respectively to the m+1 normalization calculation units to perform scheduling of the normalization calculation units; the normalization calculation units are logical gate circuits. 
         [0055]    (6) The elimination component, as illustrated in  FIG. 6 , includes a fifth processor and m×(m+1) elimination calculation units c 0, 0 , c 0, 1  . . . c m−1, m−1 , the fifth processor is connected respectively to the m×(m+1) elimination calculation units to perform scheduling of elimination calculation units and data transmission; the elimination calculation units are logical gate circuits. 
         [0056]    (7) The output port: as illustrated in  FIG. 1 , the output port V is configured to output the calculated results of linear equations over finite fields being solved. 
         [0057]    Working process of the parallel device for solving the linear equations over finite fields of this embodiment is as follows: 
         [0058]    the first processor receives the coefficient matrix B and the irreducible polynomial p(x) which is selected over the field GF(2 n ) and outputted from the input port, stores the coefficient matrix B into the memory, sets the counter as m, and sends the column containing the pivot to be found in present iteration process to the pivot finding component; the column containing the pivot to be found in present iteration process is set that: in the g th  iteration process, the column containing the pivot to be found is the g th  column, 0≦g≦m; 
         [0059]    the pivot finding component implements the operation of finding the pivot: determining whether an element with the row number of g and column number of g is a none-zero element; if yes, this element is determined to be the pivot, and the pivot finding component sends no back feed to the first processor; 
         [0060]    if not, then finding none-zero elements one by one among the elements with the row numbers of g+1 to m and the column number of g; the first none-zero element found is determined to be the pivot, and the pivot finding component sends the row number of the pivot as feedback to the processor; 
         [0061]    the first processor sends the pivot found to the partial inversion component, and sends the row where the pivot is in to the normalization component, and sends the other lines of the input coefficient matrix B to the elimination component; 
         [0062]    the partial inversion component implements partial inversion calculation, and outputs the calculated results to the normalization component and the elimination component; the partial inversion calculation specifically includes: 
         [0063]    for i=1, . . . , n−1, the third processor calculates ⊕ 2     i   ; 
         [0064]    let i′=(n−1)÷3, for k=0, 1, . . . , i′−1, calculating S k =MUL3(β 3k+1 , β 2     3k+2   , β 2     3k+1   ); MUL3 is the multiplication of three operands defined over GF(2 n ); 
         [0065]    outputting β 2     i    and S k  respectively to the normalization component and the elimination component by the partial inversion component, h=3i′+1, 3i′+2 . . . n−1. 
         [0066]    The normalization component and the elimination component respectively implement the normalization calculation and the elimination calculation, and output the results to the first processor; 
         [0067]    the normalization calculation specifically includes: 
         [0068]    l th  normalization calculation unit calculates 
         [0000]    
       
         
           
             
               
                 S 
                 
                   i 
                   ′ 
                 
               
               = 
               
                 
                   a 
                   tl 
                 
                 × 
                 
                   
                     ∏ 
                     
                       k 
                       = 
                       
                         
                           3 
                            
                           
                             i 
                             ′ 
                           
                         
                         + 
                         1 
                       
                     
                     
                       n 
                       - 
                       1 
                     
                   
                    
                   
                     β 
                     
                       2 
                       k 
                     
                   
                 
               
             
             , 
           
         
       
     
         [0000]    wherein i=(n−1)÷3, and a t1  is the element in row t and column l of the coefficient matrix; t is the row where the pivot of the present iteration is in; then calculating 
         [0000]    
       
         
           
             
               
                 a 
                 tl 
               
               = 
               
                 
                   ∏ 
                   
                     q 
                     = 
                     0 
                   
                   
                     i 
                     ′ 
                   
                 
                  
                 
                   S 
                   q 
                 
               
             
             , 
           
         
       
     
         [0000]    wherein l=0, 1, 2, . . . , m. 
         [0069]    The elimination calculation specifically includes: 
         [0070]    the elimination calculation unit numbered (k′, l) calculates 
         [0000]    
       
         
           
             
               
                 S 
                 
                   
                     i 
                     ′ 
                   
                   + 
                   1 
                 
               
               = 
               
                 { 
                 
                   
                     
                       
                         
                           
                             a 
                             
                               
                                 k 
                                 ′ 
                               
                                
                               t 
                             
                           
                           × 
                           
                             a 
                             tl 
                           
                           × 
                           
                             
                               ∏ 
                               
                                 l 
                                 = 
                                 
                                   
                                     3 
                                      
                                     
                                       i 
                                       ′ 
                                     
                                   
                                   + 
                                   1 
                                 
                               
                               
                                 n 
                                 - 
                                 1 
                               
                             
                              
                             
                               β 
                               
                                 2 
                                 l 
                               
                             
                           
                         
                         , 
                         
                           
                             j 
                             ′ 
                           
                           &lt; 
                           2 
                         
                       
                     
                   
                   
                     
                       
                         
                           
                             α 
                             tl 
                           
                           × 
                           
                             
                               ∏ 
                               
                                 l 
                                 = 
                                 
                                   
                                     3 
                                      
                                     
                                       i 
                                       ′ 
                                     
                                   
                                   + 
                                   1 
                                 
                               
                               
                                 n 
                                 - 
                                 1 
                               
                             
                              
                             
                               β 
                               
                                 2 
                                 l 
                               
                             
                           
                         
                         , 
                         
                           
                             j 
                             ′ 
                           
                           = 
                           2 
                         
                       
                     
                   
                 
                 } 
               
             
             , 
           
         
       
     
         [0000]    wherein i′=(n−1)÷3, j′=(n−1)mod3; and 
         [0071]    then calculates 
         [0000]    
       
         
           
             
               
                 a 
                 
                   
                     k 
                     ′ 
                   
                    
                   l 
                 
               
               = 
               
                 { 
                 
                   
                     
                       
                         
                           
                             
                               S 
                               
                                 
                                   i 
                                   ′ 
                                 
                                 + 
                                 1 
                               
                             
                             × 
                             
                               
                                 ∏ 
                                 
                                   q 
                                   = 
                                   0 
                                 
                                 
                                   
                                     i 
                                     ′ 
                                   
                                   - 
                                   1 
                                 
                               
                                
                               
                                 S 
                                 q 
                               
                             
                           
                           + 
                           
                             a 
                             
                               
                                 k 
                                 ′ 
                               
                                
                               l 
                             
                           
                         
                         , 
                         
                           
                             j 
                             ′ 
                           
                           &lt; 
                           2 
                         
                       
                     
                   
                   
                     
                       
                         
                           
                             
                               a 
                               
                                 
                                   k 
                                   ′ 
                                 
                                  
                                 t 
                               
                             
                             × 
                             
                               S 
                               
                                 
                                   i 
                                   ′ 
                                 
                                 + 
                                 1 
                               
                             
                             × 
                             
                               
                                 ∏ 
                                 
                                   q 
                                   = 
                                   0 
                                 
                                 
                                   
                                     i 
                                     ′ 
                                   
                                   - 
                                   1 
                                 
                               
                                
                               
                                 S 
                                 q 
                               
                             
                           
                           + 
                           
                             a 
                             
                               
                                 k 
                                 ′ 
                               
                                
                               l 
                             
                           
                         
                         , 
                         
                           
                             j 
                             ′ 
                           
                           = 
                           2 
                         
                       
                     
                   
                 
                 } 
               
             
             , 
           
         
       
     
         [0000]    wherein k′=0, 1, 2, . . . , m−1; l=0, 1, 2, . . . , m. 
         [0072]    The first processor updates the coefficient matrix B according to the calculated results of the normalization calculation and the elimination calculation, stores the updated coefficient matrix B into the memory, subtracts 1 from the value of the counter, and sets the present counter as j; then a next iteration is implemented until the counter reaches 0; 
         [0073]    if the counter is 0, then outputting the last column of the updated coefficient matrix B into the port as the solved results. 
         [0074]    Working procedure of the parallel device of the present invention is now further described taking the example of n=8 and m=12 (i.e. solving a 12×13 coefficient matrix over finite field GF(2 8 )). 
         [0075]    (1) The first processor receives a 12×13 coefficient matrix B and the irreducible polynomial p(x) which is selected over the field GF(2 n ) and outputted from the input port. 
         [0076]    Elements a(x) and p(x) of the coefficient matrix B have the following forms, respectively: 
         [0000]        a ( x )= a   7   x   7   +a   6   x   6   + . . . +a   0 ; 
         [0000]        p ( x )= x   8   +p   7   x   7   +p   6   x   6   + . . . +p   1   x+ 1; 
         [0000]    wherein a 7 , a 6 , . . . , a 0  and p 7 , p 6 , . . . , p 1  are elements in the field GF(2 n ). 
         [0077]    (2) The first processor stores the coefficient matrix B into the memory, determines the size of the coefficient matrix B to be 12×13, and sets the counter of the built-in calculator as 12; the processor sends the column (the 1 st  column) containing the pivot to be found in the first iteration process to the pivot finding component, and waits for feedback of the pivot finding component. 
         [0078]    (3) The pivot finding component firstly determines whether the element in first row and first column is a non-zero element. If yes, then this element is determined as the pivot, and the pivot finding component does not send any feedback to the first processor. If not, then finding none-zero elements one by one among the elements with row numbers of 2 to 12, and column number of 1; the finding process stops until a first non-zero element is found. The first none-zero element found is determined as the pivot, and the pivot finding component sends the row number of the pivot as feedback to the processor. 
         [0079]    (4) The first processor sends the pivot found in step (3) to the partial inversion component, and sends the row where the pivot is in to the normalization component, and sends the other lines of the input coefficient matrix B to the elimination component. 
         [0080]    (5) The partial inversion component implements partial inversion calculation: 
         [0081]    The pivot β is an element in the field GF(2 8 ); p(x) is an irreducible polynomial inputted to a selected field GF(2 n ); for i=1, . . . , 7, calculating β 2′ ; 
         [0082]    and for k=0, 1, calculating S k =MUL3(β 2     3k+1   , β 2     3k+2   , β 2     3k+3   ); MUL3 is the multiplication of three operands defined over GF(2 n ); 
         [0083]    finally, for k=0,1, outputting β 2     7    and S k  to the normalization component and the elimination component. 
         [0084]    (6) The normalization component receives the row where the pivot is in and data from the partial inversion component, and sends them to each normalization calculation unit. 
         [0085]    The l th  normalization calculation unit calculates 
         [0000]    
       
         
           
             
               
                 S 
                 
                   i 
                   ′ 
                 
               
               = 
               
                 
                   a 
                   tl 
                 
                 × 
                 
                   
                     ∏ 
                     
                       k 
                       = 
                       
                         
                           3 
                            
                           
                             i 
                             ′ 
                           
                         
                         + 
                         1 
                       
                     
                     
                       n 
                       - 
                       1 
                     
                   
                    
                   
                     β 
                     
                       2 
                       k 
                     
                   
                 
               
             
             , 
           
         
       
     
         [0000]    wherein i′=(n−1)+3, and a t1  is the element in row t and column l of the coefficient matrix B; t=1 is the row number of the pivot of the present iteration; then calculating 
         [0000]    
       
         
           
             
               
                 a 
                 tl 
               
               = 
               
                 
                   ∏ 
                   
                     q 
                     = 
                     0 
                   
                   
                     i 
                     ′ 
                   
                 
                  
                 
                   S 
                   q 
                 
               
             
             , 
           
         
       
     
         [0000]    wherein l=0, 1, . . . , 12; 
         [0086]    the normalization results are outputted to the first processor. 
         [0087]    (7) The elimination component receives data from the processor and data from the partial inversion component, and sends them to each elimination calculation unit. 
         [0088]    The elimination calculation unit numbered (k′, l) calculates S t′+1 =a k′t ×a t1 ×β 2     7   , wherein t=1; 
         [0089]    and calculates 
         [0000]    
       
         
           
             
               
                 a 
                 
                   
                     k 
                     ′ 
                   
                    
                   l 
                 
               
               = 
               
                 
                   
                     S 
                     3 
                   
                   × 
                   
                     
                       ∏ 
                       
                         q 
                         = 
                         0 
                       
                       1 
                     
                      
                     
                       S 
                       q 
                     
                   
                 
                 + 
                 
                   a 
                   
                     
                       k 
                       ′ 
                     
                      
                     l 
                   
                 
               
             
             , 
           
         
       
     
         [0000]    wherein l=0, 1, . . . , 12, k′=0, 1, 2, . . . , 11; 
         [0090]    the normalization results are outputted to the first processor, and the column number of the pivot of the next iteration (the 2 nd  column) is outputted to the pivot finding component; 
         [0091]    the pivot finding component receives elements in the column of the pivot of the next iteration, implements the operation of finding the pivot, and notifies the first processor to receive data from the normalization component and the elimination component. 
         [0092]    Upon receiving feedback from the pivot finding component, the processor subtracts 1 from the counter, receives data from the elimination component and the normalization component, stores the data into the memory, and performs a next iteration. After performing 12 times iteration operation in total and the counter returning to 0, the whole solving process is done. The solved results (last column of the coefficient matrix) are then outputted to the output port. 
         [0093]    The above embodiments are preferred embodiments of the present invention, which, however, is not intended to limit the implementation of the present invention. All of the variations, modifications, alternatives, combinations, simplifications that are not apart from the spirit of the invention shall be deemed as equivalences to those skilled in the art, and are within the protection scope of the present invention.