graphics_toolkit("fltk");

GoCs      = load ("GoCcoordinates.sorted.dat" );
GoCadends = load ("GoCadendcoordinates.sorted.dat" );
#GoCbdends = load ("GoCbdendcoordinates.sorted.dat" );
apoints = size(GoCadends)(2)/3
#bpoints = size(GoCbdends)(2)/3

GCs            = load ("GCcoordinates.sorted.dat" );
PFs            = load ("PFcoordinates.dat" );
PFtoGoCsources = load ("PFtoGoCsources.dat" );
PFtoGoCtargets = load ("PFtoGoCtargets.dat" );

myGoCind = 0; #randi(size(GoCs)(1),1);
myGoCcoords = [GoCs(myGoCind+1,1) GoCs(myGoCind+1,2) GoCs(myGoCind+1,3)];
myGoCacoords = GoCadends(myGoCind+1,:);
#myGoCbcoords = GoCbdends(myGoCind+1,:);

myPFtoGoCinds    = find (PFtoGoCtargets == myGoCind);
myPFtoGoCsources = PFtoGoCsources(myPFtoGoCinds);

plot3 (GoCs(myGoCind+1,1),GoCs(myGoCind+1,2),GoCs(myGoCind+1,3),'r.','markersize',5);
hold on;
for i=1:apoints
  x = myGoCacoords(:,(i-1)*3+1);
  y = myGoCacoords(:,(i-1)*3+2);
  z = myGoCacoords(:,(i-1)*3+3);
  plot3(x,y,z,'g.');
endfor       

#for i=1:bpoints
#  x = myGoCbcoords(:,(i-1)*3+1);
#  y = myGoCbcoords(:,(i-1)*3+2);
#  z = myGoCbcoords(:,(i-1)*3+3);
#  plot3(x,y,z,'b.');
#endfor       
for i = 1:length(myPFtoGoCsources)
    id = myPFtoGoCsources (i)
    myPFs = find (PFs(:,1) == id);
    plot3 (PFs(myPFs,2),PFs(myPFs,3),PFs(myPFs,4),'r');
    plot3 (GCs(id+1,1),GCs(id+1,2),GCs(id+1,3),'b.');
endfor

PFs = load ("PFcoordinates.dat" );

nGC = 1000;
myGCinds = randi (size(GCs)(1),nGC,1);

for i=1:nGC
  i
  myPFinds = find (PFs(:,1) == myGCinds(i));
  plot3 (PFs(myPFinds,2),PFs(myPFinds,3),PFs(myPFinds,4),'b');
endfor

title ("Parallel Fiber and Golgi Cell connections",'fontsize',24);