Changeset 1628b14 in git
- Timestamp:
- Mar 24, 2011, 10:25:14 AM (12 years ago)
- Branches:
- (u'jengelh-datetime', 'ceac47cbc86fe4a15902392bdbb9bd2ae0ea02c6')(u'spielwiese', 'a800fe4b3e9d37a38c5a10cc0ae9dfa0c15a4ee6')
- Children:
- 2ab830a563f54a42a41eb941f715cec3422755af
- Parents:
- a87b34163475a42f2c9e60652769c83ac2b67427
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
Singular/LIB/paraplanecurves.lib
ra87b34 r1628b14 759 759 " 760 760 { 761 "#############################################"; 762 int totalTime = timer; 761 763 poly dhf = subst(f,var(3),1); 762 764 def Roriginal = basering; … … 780 782 def QF = ring(rl); // make var(1) transcendental 781 783 list LIntB; 784 int t = timer; 782 785 if(DATA[1] <= 4) // use normalization algorithm 783 786 { … … 796 799 setring Rdummy; 797 800 } 801 t = timer - t; 802 "### adjointIdealIB, IB-time:", t, "msec"; 798 803 int i,j,k,l; 799 804 ideal IB = LIntB[1]; … … 809 814 poly d = imap(Rdummy,d); 810 815 M=1/d*M; 816 t = timer; 811 817 list LUM = ludecomp(M); 818 t = timer - t; 819 "### adjointIdealIB, time for 1st ludecomp:", t, "msec"; 812 820 list LS; 813 821 matrix dummyvector[sL][1]; … … 818 826 list Iij, empty; 819 827 matrix Gij[sL][sL]; 828 t = timer; 820 829 for(i = 1; i <= sL; i++) 821 830 { … … 848 857 } 849 858 } 859 t = timer - t; 860 "### adjointIdealIB, time for linear solving:", t, "msec"; 861 t = timer; 850 862 LUM = ludecomp(Tr); 863 t = timer - t; 864 "### adjointIdealIB, time for 2nd ludecomp:", t, "msec"; 851 865 setring Rred; 852 866 poly d2f = imap(Rdummy,d2f); … … 860 874 setring Roriginal; 861 875 ideal AI = imap(QF,LL); 876 totalTime = timer - totalTime; 877 "### adjointIdealIB, total time:", totalTime, "msec"; 862 878 return(AI); 863 879 } … … 866 882 static proc adjointIdealIQ(poly f, list choices) 867 883 { 884 "#############################################"; 885 int totalTime = timer; 868 886 poly dhf = subst(f,var(3),1); 869 887 def Roriginal = basering; … … 894 912 } 895 913 t = timer - t; 896 "#############################################"; 897 "### IB-time:", t, "msec"; 914 "### adjointIdealIQ, IB-time:", t, "msec"; 898 915 int i,j,k,l; 899 916 ideal IB = LIntB[1]; … … 904 921 ideal AI = quotient(fd, IBf); 905 922 t = timer - t; 906 "### quot-time:", t, "msec";923 "### adjointIdealIQ, quot-time:", t, "msec"; 907 924 //"#### IB:"; IB; 908 925 //"#### d:", d; … … 910 927 ideal AI = imap(Rdummy,AI); 911 928 //"#### AI:"; AI; 929 totalTime = timer - totalTime; 930 "### adjointIdealIQ, total time:", totalTime, "msec"; 912 931 return(AI); 913 932 }
Note: See TracChangeset
for help on using the changeset viewer.