← Back to team overview

dolfin team mailing list archive

Re: Function file io

 

On 29 September 2011 09:33, Anders Logg <logg@xxxxxxxxx> wrote:
> On Thu, Sep 29, 2011 at 12:34:57AM +0100, Garth N. Wells wrote:
>> After a long absence, Function i/o to file has made a come back. It
>> was particularly necessary to support running code in parallel that
>> reads in data for Functions (such as the advection-diffusion demo). We
>> have more demos running in parallel now, and it's now possible to
>> restart parallel computations reliably (and even on a different number
>> of processes).
>>
>> When sending a Function to an XML file, there is an entry for each
>> vector (dof) entry, e.g.
>>
>>     <dof global_index="2" value="0" cell_index="3" local_dof_index="0" />
>>
>> 'global_index' is the vector index, 'value' is the entry value,
>> 'cell_index' is a cell that the dof belongs to and 'local_dof_index'
>> is the local dof number relative to the given cell. With this
>> information, it's possible to read in and use function data
>> independently of the new dof map. If anyone has suggestions or
>> objections to this format, speak up now.
>
> Is the cell index global to the mesh?
>

Yes.

> For consistency, I'd suggest the following names:
>
>  <dof index="2" value="0" cell_index="3" cell_dof_index="0" />
>

OK.

Garth

> --
> Anders
>
>
>> When reading in a Function, the Function must have a FunctionSpace,
>> and the vector must already have been sized (since the partitioning of
>> the vector in parallel will depend on the mesh partition).
>>
>> Garth
>>
>> _______________________________________________
>> Mailing list: https://launchpad.net/~dolfin
>> Post to     : dolfin@xxxxxxxxxxxxxxxxxxx
>> Unsubscribe : https://launchpad.net/~dolfin
>> More help   : https://help.launchpad.net/ListHelp
>


References