# generate four hyperstreamlines

# get the supporting scripts
source vtkInt.tcl
source vtkInclude.tcl

# create tensor ellipsoids
vtkRenderMaster rm;

# Now create the RenderWindow, Renderer and interactive renderer
#
set renWin [rm MakeRenderWindow];
set ren1 [$renWin MakeRenderer];
set iren [$renWin MakeRenderWindowInteractor];

#
# Create tensor ellipsoids
#
# generate tensors
vtkPointLoad ptLoad;
    ptLoad SetLoadValue 100.0;
    ptLoad SetSampleDimensions 30 30 30;
    ptLoad ComputeEffectiveStressOn;
    ptLoad SetModelBounds -10 10 -10 10 -10 10;

# Generate hyperstreamlines
vtkHyperStreamline s1;
    s1 SetInput [ptLoad GetOutput];
    s1 SetStartPosition 9 9 -9;
    s1 IntegrateMinorEigenvector;
    s1 SetMaximumPropagationDistance 18.0;
    s1 SetIntegrationStepLength 0.1;
    s1 SetStepLength 0.01;
    s1 SetRadius 0.25;
    s1 SetNumberOfSides 18;
    s1 SetIntegrationDirection $VTK_INTEGRATE_BOTH_DIRECTIONS;
    s1 Update;
# Map hyperstreamlines
vtkLogLookupTable lut;
    lut SetHueRange .6667 0.0;
vtkPolyMapper s1Mapper;
    s1Mapper SetInput [s1 GetOutput];
    s1Mapper SetLookupTable lut;
    ptLoad Update;#force update for scalar range
    eval s1Mapper SetScalarRange [[ptLoad GetOutput] GetScalarRange];
vtkActor s1Actor;
    s1Actor SetMapper s1Mapper;

vtkHyperStreamline s2;
    s2 SetInput [ptLoad GetOutput];
    s2 SetStartPosition -9 -9 -9;
    s2 IntegrateMinorEigenvector;
    s2 SetMaximumPropagationDistance 18.0;
    s2 SetIntegrationStepLength 0.1;
    s2 SetStepLength 0.01;
    s2 SetRadius 0.25;
    s2 SetNumberOfSides 18;
    s2 SetIntegrationDirection $VTK_INTEGRATE_BOTH_DIRECTIONS;
    s2 Update;
vtkPolyMapper s2Mapper;
    s2Mapper SetInput [s2 GetOutput];
    s2Mapper SetLookupTable lut;
    ptLoad Update;#force update for scalar range
    eval s2Mapper SetScalarRange [[ptLoad GetOutput] GetScalarRange];
vtkActor s2Actor;
    s2Actor SetMapper s2Mapper;

vtkHyperStreamline s3;
    s3 SetInput [ptLoad GetOutput];
    s3 SetStartPosition 9 -9 -9;
    s3 IntegrateMinorEigenvector;
    s3 SetMaximumPropagationDistance 18.0;
    s3 SetIntegrationStepLength 0.1;
    s3 SetStepLength 0.01;
    s3 SetRadius 0.25;
    s3 SetNumberOfSides 18;
    s3 SetIntegrationDirection $VTK_INTEGRATE_BOTH_DIRECTIONS;
    s3 Update;
vtkPolyMapper s3Mapper;
    s3Mapper SetInput [s3 GetOutput];
    s3Mapper SetLookupTable lut;
    ptLoad Update;#force update for scalar range
    eval s3Mapper SetScalarRange [[ptLoad GetOutput] GetScalarRange];
vtkActor s3Actor;
    s3Actor SetMapper s3Mapper;

vtkHyperStreamline s4;
    s4 SetInput [ptLoad GetOutput];
    s4 SetStartPosition -9 9 -9;
    s4 IntegrateMinorEigenvector;
    s4 SetMaximumPropagationDistance 18.0;
    s4 SetIntegrationStepLength 0.1;
    s4 SetStepLength 0.01;
    s4 SetRadius 0.25;
    s4 SetNumberOfSides 18;
    s4 SetIntegrationDirection $VTK_INTEGRATE_BOTH_DIRECTIONS;
    s4 Update;
vtkPolyMapper s4Mapper;
    s4Mapper SetInput [s4 GetOutput];
    s4Mapper SetLookupTable lut;
    ptLoad Update;#force update for scalar range
    eval s4Mapper SetScalarRange [[ptLoad GetOutput] GetScalarRange];
vtkActor s4Actor;
    s4Actor SetMapper s4Mapper;

#
# plane for context
#
vtkStructuredPointsGeometryFilter g;
    g SetInput [ptLoad GetOutput];
    g SetExtent 0 100 0 100 0 0;
    g Update;#for scalar range
vtkPolyMapper gm;
    gm SetInput [g GetOutput];
    eval gm SetScalarRange [[g GetOutput] GetScalarRange];
vtkActor ga;
    ga SetMapper gm;

#
# Create outline around data
#
vtkOutlineFilter outline;
    outline SetInput [ptLoad GetOutput];
vtkPolyMapper outlineMapper;
    outlineMapper SetInput [outline GetOutput];
vtkActor outlineActor;
    outlineActor SetMapper outlineMapper;
    eval [outlineActor GetProperty] SetColor 0 0 0;
#
# Create cone indicating application of load
#
vtkConeSource coneSrc;
    coneSrc  SetRadius .5;
    coneSrc  SetHeight 2;
vtkPolyMapper coneMap
    coneMap SetInput [coneSrc GetOutput];
vtkActor coneActor;
    coneActor SetMapper coneMap;    
    coneActor SetPosition 0 0 11;
    coneActor RotateY 90;
    eval [coneActor GetProperty] SetColor 1 0 0;

vtkCamera camera;
    camera SetFocalPoint 0.113766 -1.13665 -1.01919;
    camera SetPosition -29.4886 -63.1488 26.5807;
    camera CalcViewPlaneNormal;
    camera SetViewAngle 24.4617;
    camera SetViewUp 0.17138 0.331163 0.927879;

$ren1 AddActors s1Actor;
$ren1 AddActors s2Actor;
$ren1 AddActors s3Actor;
$ren1 AddActors s4Actor;
$ren1 AddActors outlineActor;
$ren1 AddActors coneActor;
$ren1 AddActors ga;
$ren1 SetBackground 1.0 1.0 1.0;
$ren1 SetActiveCamera camera;

$renWin SetSize 500 500;
$renWin Render;
$iren SetUserMethod {wm deiconify .vtkInteract};

#$renWin SetFileName Hyper.tcl.ppm;
#$renWin SaveImageAsPPM;

# prevent the tk window from showing up then start the event loop
wm withdraw .
