Changeset 0b59f5 in git for Singular/LIB/linalg.lib


Ignore:
Timestamp:
Dec 13, 1999, 4:33:50 PM (24 years ago)
Author:
Olaf Bachmann <obachman@…>
Branches:
(u'spielwiese', 'fe61d9c35bf7c61f2b6cbf1b56e25e2f08d536cc')
Children:
82ed244cff3158b7a79a6eed4e0548d485960cfe
Parents:
925cab8e04ecb2f3c2e451b913c11bcac28b2374
Message:
* GMG's changes


git-svn-id: file:///usr/local/Singular/svn/trunk@4006 2c84dea3-7e68-4137-9b89-c4e89433aadc
File:
1 edited

Legend:

Unmodified
Added
Removed

  • Singular/LIB/linalg.lib

    r925cab r0b59f5  
    33//////////////////////////////////////////////////////////////////////////////
    44
    5 version="$Id: linalg.lib,v 1.3 1999-10-25 12:45:16 obachman Exp $";
     5version="$Id: linalg.lib,v 1.4 1999-12-13 15:33:48 obachman Exp $";
    66info="
    77LIBRARY:  linalg.lib    PROCEDURES FOR ALGORITHMIC LINEAR ALGEBRA
    … …  
    2121 gaussred(A);       gaussian reduction: P*A=U*S, S is a row reduced form of A
    2222 gaussred_pivot(A); gaussian reduction: P*A=U*S, S is a row reduced form of A
    23  gnf(A);            gauss normal form of A
    24  rk(A);             rank of A
    25  UDO(A);            P*A=U*D*O 
     23 gauss_nf(A);            gauss normal form of A
     24 mat_rk(A);             rank of A
     25 U_D_O(A);            P*A=U*D*O 
    2626 pos_def(A,i);      test of positive symmetric matrizes.       
    2727";
    … …  
    824824
    825825//////////////////////////////////////////////////////////////////////////////
    826 proc gnf(matrix A)
    827 "USAGE:    gnf(A); any constant matrix A
     826proc gauss_nf(matrix A)
     827"USAGE:    gauss_nf(A); any constant matrix A
    828828RETURN:   matrix; gauss normal form of A
    829 EXAMPLE:  example gnf; shows an example"
     829EXAMPLE:  example gauss_nf; shows an example"
    830830{
    831831  list Z;
    … …  
    841841ring r=7,(x),dp;
    842842matrix A[4][4]= 1,4,4,7,2,5,5,4,4,1,1,3,0,2,2,7;
    843 print(gnf(A));
    844 }
    845 
    846 //////////////////////////////////////////////////////////////////////////////
    847 proc rk(matrix A)
    848 "USAGE:    rk(A); any constant matrix A
     843print(gauss_nf(A));
     844}
     845
     846//////////////////////////////////////////////////////////////////////////////
     847proc mat_rk(matrix A)
     848"USAGE:    mat_rk(A); any constant matrix A
    849849RETURN:   int; rank of A
    850 EXAMPLE:  example rk; shows an example"
     850EXAMPLE:  example mat_rk; shows an example"
    851851{
    852852  list Z;
    … …  
    862862ring r=7,(x),dp;
    863863matrix A[4][4]= 1,4,4,7,2,5,5,4,4,1,1,3,0,2,2,7;
    864 rk(A);
    865 }
    866 
    867 //////////////////////////////////////////////////////////////////////////////
    868 proc UDO(matrix A)
     864mat_rk(A);
     865}
     866
     867//////////////////////////////////////////////////////////////////////////////
     868proc U_D_O(matrix A)
    869869"USAGE:     gaussred(A);   constant invertible matrix A
    870870RETURN:    list Z:  Z[1]=P , Z[2]=U , Z[3]=D , Z[4]=O
    … …  
    922922               -10,10,-9,5,-1,-6,5;
    923923print(A);
    924 list Z=UDO(A);
     924list Z=U_D_O(A);
    925925"P=";print(Z[1]);
    926926"U=";print(Z[2]);
    … …  
    956956  }
    957957 
    958   Z=UDO(A);
     958  Z=U_D_O(A);
    959959
    960960  if(Z[1]==-1){return(0);} //A not regular , therefor not pos.definit
    … …  
    10231023  }
    10241024 
    1025   r=rk(A);
     1025  r=mat_rk(A);
    10261026       
    10271027  //1. b constant vector
Note: See TracChangeset for help on using the changeset viewer.