dolfin team mailing list archive

[Patrick Riesen <question147820@xxxxxxxxxxxxxxxxxxxxx>]

Question #147820 on DOLFIN changed:
https://answers.launchpad.net/dolfin/+question/147820

Patrick Riesen proposed the following answer:
i'm sorry, the last line should be:
solution = ah.vtk2array(vtk_values)

patrick

On 03/07/2011 09:16 AM, Patrick Riesen wrote:
> in python, I use this to read in a vtu file to an array:
> this doesn't need the GetArray("name")
>
> # -------------------------------------------------
> import vtk as v
> from enthought.tvtk import array_handler as ah
>
> datareader = v.vtkXMLUnstructuredGridReader()
> datareader.SetFileName(filename)
> datareader.Update()
>
> data = datareader.GetOutput()
> # if you have a vector (e.g. velocity):
> vtk_values = data.GetPointData().GetVectors()
> # if you have a scalar:
> # vtk_values = data.GetPointData().GetScalars()
>
> # this is convenient to convert vtkarrays to numpy arrays
> solution = ah.vtk2array(data)
> # ------------------------------------------------------
>
>
> regards,
> patrick
>
>
> On 03/04/2011 09:30 PM, JJ Jiang wrote:
>> Question #147820 on DOLFIN changed:
>> https://answers.launchpad.net/dolfin/+question/147820
>>
>> Status: Answered => Open
>>
>> JJ Jiang is still having a problem:
>> Hi Johan and Kevin,
>> Thanks for the notes. I think you can dump simulated temperature or
>> mass to raw files. That can be manipulated when the mesh data are also
>> given either in a vtu format or in a raw format.
>> By the way, I used Paraview to see the simulated "temperature" scalar
>> has been attached to a vtu. Surprisingly, I couldn't read it out by
>> its name. I don't know why.
>> I guess Johan's suggestions will give users more flexibilities to
>> handle the data.
>> Thanks again to both of you.
>> Best regards,
>> --JJ
>>
>> On 03/04/11, Johan Hake wrote:
>>> Your question #147820 on DOLFIN changed:
>>> https://answers.launchpad.net/dolfin/+question/147820
>>>
>>> Johan Hake proposed the following answer:
>>> Maybe Viper could provide a high level interface to read a vtu file
>>> into a
>>> NumPy array, and a mesh into a coordinate and cell arrays?
>>>
>>> Then we eventually can use this to pass data to Viper for plotting in
>>> the C++
>>> interface. This might fix the parallel ploting issue we have.
>>>
>>> It would also be cool to implement a helper function in the Python
>>> interface
>>> of DOLFIN that uses the potential added functionality in Viper to
>>> generate a
>>> Function, with or without a given mesh. A FunctionSpace could be
>>> automatically
>>> generated based on the dimensions of the data (is this true?) If a
>>> mesh is not
>>> given we generate it based on the stored mesh data from the file.
>>>
>>> Johan
>>>
>>>
>>> On Friday March 4 2011 11:30:52 Kent-Andre Mardal wrote:
>>>> Question #147820 on DOLFIN changed:
>>>> https://answers.launchpad.net/dolfin/+question/147820
>>>>
>>>> Status: Open => Answered
>>>>
>>>> Kent-Andre Mardal proposed the following answer:
>>>> You should be able to use GetPointData().GetArray(). Have you looked
>>>> at the
>>>> name of the variable in the vtu file? It is in ASCII format.
>>>>
>>>> I attach something I did earlier.
>>>>
>>>> Kent
>>>>
>>>> On 4 March 2011 19:34, JJ Jiang
>>>>
>>>> <question147820@xxxxxxxxxxxxxxxxxxxxx>wrote:
>>>>> New question #147820 on DOLFIN:
>>>>> https://answers.launchpad.net/dolfin/+question/147820
>>>>>
>>>>> Hi all,
>>>>>
>>>>> I am running into a problem. I run the advection-diffusion demo.
>>>>> All the
>>>>>
>>>>> simulated mass transport quantities u were save to a vtk PVD file
>>>>> with an
>>>>> unstructured mesh.
>>>>>
>>>>> My question is how to read out the u from the associated individual
>>>>> vtu
>>>>> files. I thought the u is attached to either PointData or CellData so
>>>>> that I can read out the scalar u using one of the following:
>>>>> grid->GetPointData()->GetArray('name');
>>>>> or
>>>>> grid->GetCellData()->GetArray('name');
>>>>>
>>>>> But my attempts have failed. Please provide your advice. Thanks so
>>>>> much
>>>>> for your help!
>>>>>
>>>>> Best regards,
>>>>> --JJ
>>>>>
>>>>>
>>>>> --
>>>>> You received this question notification because you are a member of
>>>>> DOLFIN Team, which is an answer contact for DOLFIN.
>>>>>
>>>>> _______________________________________________
>>>>> Mailing list: https://launchpad.net/~dolfin
>>>>> Post to : dolfin@xxxxxxxxxxxxxxxxxxx
>>>>> Unsubscribe : https://launchpad.net/~dolfin
>>>>> More help : https://help.launchpad.net/ListHelp
>>>
>>> --
>>> If this answers your question, please go to the following page to let us
>>> know that it is solved:
>>> https://answers.launchpad.net/dolfin/+question/147820/+confirm?answer_id=1
>>>
>>>
>>> If you still need help, you can reply to this email or go to the
>>> following page to enter your feedback:
>>> https://answers.launchpad.net/dolfin/+question/147820
>>>
>>> You received this question notification because you are a direct
>>> subscriber of the question.
>>
>> --
>> +++++++++++++++++++++++++++++++++++
>> Jingfeng Jiang (JJ)
>> University of Wisconsin
>> Medical Physics Department
>> 1530 Medical Sciences Center
>> 1300 University Avenue
>> Madison, WI 53706
>> Phone: (608) 262-4197
>> FAX: (608) 262-2413
>> +++++++++++++++++++++++++++++++++++
>>
>> You received this question notification because you are a member of
>> DOLFIN Team, which is an answer contact for DOLFIN.
>>
>> _______________________________________________
>> Mailing list: https://launchpad.net/~dolfin
>> Post to : dolfin@xxxxxxxxxxxxxxxxxxx
>> Unsubscribe : https://launchpad.net/~dolfin
>> More help : https://help.launchpad.net/ListHelp
>
>
> _______________________________________________
> Mailing list: https://launchpad.net/~dolfin
> Post to : dolfin@xxxxxxxxxxxxxxxxxxx
> Unsubscribe : https://launchpad.net/~dolfin
> More help : https://help.launchpad.net/ListHelp

-- 
You received this question notification because you are a member of
DOLFIN Team, which is an answer contact for DOLFIN.