graphics_toolkit("fltk");

GoCs       = load ("GoCcoordinates.sorted.dat" );
GoCadends  = load ("GoCadendcoordinates.sorted.dat" );
#GoCbdends = load ("GoCbdendcoordinates.sorted.dat" );
GoCaxons   = load ("GoCaxoncoordinates.sorted.dat" );

nGoC = 1;

apoints    = size(GoCadends)(2)/3
#bpoints   = size(GoCbdends)(2)/3
axonpoints = size(GoCaxons)(2)/3

myGoCinds = randi(size(GoCs)(1),nGoC,1);

myGoCcoords = [GoCs(myGoCinds,1) GoCs(myGoCinds,2) GoCs(myGoCinds,3)];

myGoCacoords    = GoCadends(myGoCinds,:);
#myGoCbcoords   = GoCbdends(myGoCinds,:);
myGoCaxoncoords = GoCaxons(myGoCinds,:);

plot3 (GoCs(myGoCinds,1),GoCs(myGoCinds,2),GoCs(myGoCinds,3),'r.');
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:axonpoints
  x = myGoCaxoncoords(:,(i-1)*3+1);
  y = myGoCaxoncoords(:,(i-1)*3+2);
  z = myGoCaxoncoords(:,(i-1)*3+3);
  plot3(x,y,z,'k.');
endfor       

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

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

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