Changeset 2473c9 in git
 Timestamp:
 Apr 3, 2019, 12:39:17 PM (4 years ago)
 Branches:
 (u'spielwiese', '0d6b7fcd9813a1ca1ed4220cfa2b104b97a0a003')
 Children:
 56ffd981992fd934df3847d95789268d842efbcc
 Parents:
 a74fa8837632268cfdb724eabc42506f5e99818e
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

Singular/LIB/realclassify.lib
ra74fa8 r2473c9 1333 1333 ring R = 0,x,dp; 1334 1334 poly k = fetch(br,k); 1335 int n = nrroots(k);1335 int nn = nrroots(k); 1336 1336 kill R; 1337 1337 setring br; 1338 if(n ==1)1338 if(nn==1) 1339 1339 { 1340 1340 typeofsing1 = "J[10]+"; … … 1364 1364 } 1365 1365 if(deg(MP)==1) 1366 1367 1368 1369 1370 1371 1372 1373 1374 1375 1376 1377 1378 1379 1380 1381 1382 1383 1384 1385 1386 1387 1388 1389 1390 1391 1392 return(F1)1366 { 1367 number a = number(subst(MP,var(1),0)); 1368 setring br1; 1369 number a = fetch(br,a); 1370 } 1371 else 1372 { 1373 string mord = ringlist(br1)[3][1][1]; 1374 ring R1 = (0,a),(var(1),var(2)),ds; 1375 list L1 = ringlist(R1); 1376 L1[3][1][1]=mord; 1377 def R2 = ring(L1); 1378 setring R2; 1379 kill R1; 1380 poly ff = fetch(br,MP); 1381 minpoly = number(subst(ff,var(1),a)); 1382 } 1383 Poly nfe1 = var(1)^3+a*var(1)^2*var(2)^2+var(1)*var(2)^4; 1384 number n1 = fetch(br,n1); 1385 number n2 = fetch(br,n2); 1386 interval1 = list(n1,n2); 1387 NormalFormEquation F1; 1388 F1.normalFormEquation = nfe1; 1389 F1.parameters = list(list(a*var(1)^2*var(2)^2,interval1)); 1390 F1.singularityType = typeofsing1; 1391 setring br; 1392 return(F1); 1393 1393 } 1394 1394 else … … 1462 1462 kill r; 1463 1463 kill j; 1464 int k = size(r1);1464 int kk = size(r1); 1465 1465 setring R1; 1466 1466 list r2; 1467 for(int i =1; i<=k;i++)1467 for(int ii=1; ii<=kk;ii++) 1468 1468 { 1469 1469 number m; 1470 m=number(r1[i ])/(10^(epsilon));1470 m=number(r1[ii])/(10^(epsilon)); 1471 1471 r2 = insert(r2,m); 1472 1472 kill m; 1473 1473 } 1474 kill i ;1474 kill ii; 1475 1475 list L1 = ringlist(br1); 1476 1476 number t = (number(epsilon))^(1); … … 1504 1504 kill r; 1505 1505 kill j; 1506 k = size(r1);1507 1506 setring R1; 1507 kk = size(r1); 1508 1508 list r2; 1509 for(int i =1; i<=k;i++)1509 for(int ii=1; ii<=kk;ii++) 1510 1510 { 1511 1511 number m; 1512 m=number(r1[i ])/(10^(epsilon));1512 m=number(r1[ii])/(10^(epsilon)); 1513 1513 r2 = insert(r2,m); 1514 1514 kill m; 1515 1515 } 1516 kill i ;1516 kill ii; 1517 1517 number t = (number(epsilon))^(1); 1518 1518 number n1 = sturm(pplus,r2[1]t,r2[1]+t); … … 1524 1524 kill R1; 1525 1525 list r2; 1526 for(int i=1; i<=k;i++) 1526 kk = size(r1); 1527 for(int ii=1; ii<=kk;ii++) 1527 1528 { 1528 1529 number m; 1529 m=number(r1[i ])/(10^(epsilon));1530 m=number(r1[ii])/(10^(epsilon)); 1530 1531 r2 = insert(r2,m); 1531 1532 kill m; 1532 1533 } 1533 kill i ;1534 kill ii; 1534 1535 number t = (number(epsilon))^(1); 1535 1536 typeofsing1 = "J[10]";
Note: See TracChangeset
for help on using the changeset viewer.