Singular

D.5.9.6 prepEmbDiv

Procedure from library `reszeta.lib` (see reszeta_lib).

Usage:
prepEmbDiv(L[,a]);
L = list of rings
a = integer

Assume:
L is output of resolution of singularities

Compute:
if a is not present: exceptional divisors including components of the strict transform
otherwise: only exceptional divisors

Return:
list of Q-irreducible exceptional divisors (embedded case)

Example:
 ```LIB "reszeta.lib"; ring R=0,(x,y,z),dp; ideal I=x2+y2+z11; list L=resolve(I); prepEmbDiv(L); ==> [1]: ==> [1]: ==> 2,1 ==> [2]: ==> 3,1 ==> [3]: ==> 4,1 ==> [2]: ==> [1]: ==> 4,2 ==> [2]: ==> 5,2 ==> [3]: ==> 6,2 ==> [3]: ==> [1]: ==> 6,3 ==> [2]: ==> 7,3 ==> [3]: ==> 8,3 ==> [4]: ==> [1]: ==> 8,4 ==> [2]: ==> 9,4 ==> [3]: ==> 10,4 ==> [5]: ==> [1]: ==> 10,5 ==> [2]: ==> 11,5 ==> [3]: ==> 12,5 ==> [4]: ==> 13,5 ==> [5]: ==> 15,5 ==> [6]: ==> 17,5 ==> [6]: ==> [1]: ==> 12,6 ==> [2]: ==> 13,6 ==> [3]: ==> 14,6 ==> [4]: ==> 16,6 ==> [7]: ==> [1]: ==> 14,7 ==> [2]: ==> 15,7 ==> [3]: ==> 16,7 ==> [4]: ==> 17,7 ==> [8]: ==> [1]: ==> 3,2 ==> [2]: ==> 4,3 ==> [3]: ==> 6,4 ==> [4]: ==> 8,5 ==> [5]: ==> 10,6 ==> [6]: ==> 14,8 ==> [7]: ==> 15,8 ==> [8]: ==> 16,8 ==> [9]: ==> 17,8 ```