/*  K3 surface: Quartic CP3 Polynomial in 16D Veronese: Include file
  NOTE: Bug in NDviewer makes XI steps 65 -> 21, TH steps 33-> 11,
    so it divides by 3?  
 */


#define XIMAX    6.0    		// range of exponential
#define XISTEPS  65    		//  exponential steps must be ODD
#define XISCALE  0.05 		//  scale of sinh
#define THSTEPS  33    		//  theta steps must be ODD

/* 

a:region 1 xi   phi=pi/2, th=0:pi/2
b:region 1 rho  phi=pi/2, th=0:pi/2
c:region 2 xi   phi=0:pi/2, th=0
d:region 2 rho  phi=0:pi/2, th=0
e:region 3 xi   phi=0:pi/2, th=pi/2
f:region 3 rho  phi=0:pi/2, th=pi/2
g: polar region  xi and rho  phi=0:pi/2, th=pi/2

./makeplugins K3-000.cpp   full 4 manifold

   two-manifold boundaries, I think!

./makeplugins K3-000a.cpp
./makeplugins K3-000b.cpp
./makeplugins K3-000c.cpp
./makeplugins K3-000d.cpp
./makeplugins K3-000e.cpp
./makeplugins K3-000f.cpp
./makeplugins K3-000g.cpp
./makeplugins K3-000h.cpp
./makeplugins K3-000i.cpp
./makeplugins K3-000j.cpp
./makeplugins K3-000k.cpp
./makeplugins K3-000l.cpp
./makeplugins K3-000m.cpp




 */