To check out this repository please hg clone the following URL, or open the URL using EasyMercurial or your preferred Mercurial client.

Statistics Download as Zip
| Branch: | Revision:

root / _FullBNT / BNT / general / partition_matrix_vec_3.m @ 8:b5b38998ef3b

History | View | Annotate | Download (547 Bytes)

1 
function [A1, A2, B1, B2, C11, C12, C21, C22] = partition_matrix_vec_3(A, B, C, n1, n2, bs)
2 

    		
  

  

    
      3
    
    
      
dom = myunion(n1, n2);


    		
  

  

    
      4
    
    
      
n1i = block(find_equiv_posns(n1, dom), bs(dom));


    		
  

  

    
      5
    
    
      
n2i = block(find_equiv_posns(n2, dom), bs(dom));


    		
  

  

    
      6
    
    
      

    		
  

  

    
      7
    
    
      

    		
  

  

    
      8
    
    
      
    A1 = A(n1i);


    		
  

  

    
      9
    
    
      
    A2 = A(n2i);


    		
  

  

    
      10
    
    
      
    if isempty(B)


    		
  

  

    
      11
    
    
      
        B1 = zeros(size(n1i, 2),size(B, 2));


    		
  

  

    
      12
    
    
      
        B2 = zeros(size(n2i, 2),size(B, 2));


    		
  

  

    
      13
    
    
      
    else


    		
  

  

    
      14
    
    
      
        B1 = B(n1i, :);


    		
  

  

    
      15
    
    
      
        B2 = B(n2i, :);


    		
  

  

    
      16
    
    
      
    end


    		
  

  

    
      17
    
    
      
    


    		
  

  

    
      18
    
    
      
    


    		
  

  

    
      19
    
    
      
    C11 = C(n1i, n1i);


    		
  

  

    
      20
    
    
      
    C12 = C(n1i, n2i);


    		
  

  

    
      21
    
    
      
    C21 = C(n2i, n1i);


    		
  

  

    
      22
    
    
      
    C22 = C(n2i, n2i);