#!/usr/bin/env python

import vtk

sphere1 = vtk.vtkSphereSource()
sphere1.SetCenter(1, 0, 0)
sphere1.SetRadius(0.5)
sphere2 = vtk.vtkSphereSource()
sphere2.SetRadius(0.5)
## sphere1.Update()
## sphere2.Update()

## multiblock = vtk.vtkMultiBlockDataSet()
## multiblock.SetNumberOfBlocks(2)
## multiblock.SetBlock(0, sphere1.GetOutput())
## multiblock.SetBlock(1, sphere2.GetOutput())

## sphereMapper = vtk.vtkCompositePolyDataMapper2()
## sphereMapper.SetInputConnection(multiblock.GetProducerPort())

multiblock = vtk.vtkAppendPolyData()
multiblock.AddInput(sphere1.GetOutput())
multiblock.AddInput(sphere2.GetOutput())

sphereMapper = vtk.vtkPolyDataMapper()
sphereMapper.SetInputConnection(multiblock.GetOutputPort())

## sphereActor = vtk.vtkActor()
sphereActor = vtk.vtkQuadricLODActor()
sphereActor.SetMapper(sphereMapper)

# following line not needed with patch:
# http://public.kitware.com/Bug/view.php?id=8659
sphereActor.GetProperty()

ren = vtk.vtkRenderer()
ren.AddActor(sphereActor)
ren.ResetCamera()

# Create the RenderWindow and RenderWindowInteractor
renWin = vtk.vtkRenderWindow()
renWin.AddRenderer(ren)
iren = vtk.vtkRenderWindowInteractor()
iren.SetRenderWindow(renWin)

iren.Start()