dolfin team mailing list archive

["Alessio Quaglino" <alessio@xxxxxx>]

>
>
> Alessio Quaglino wrote:
>>> Alessio Quaglino wrote:
>>>>> Alessio Quaglino wrote:
>>>>>> With "the mesh remains constant" is meant that the vertices cannot
>>>>>> move
>>>>>> or
>>>>>> that there should not be any connectivity change?
>>>>>>
>>>>> No connectivity changes.
>>>> I see that in
>>>>
>>>> void Assembler::assemble(GenericTensor& A, const Form& form,
>>>>                          bool reset_tensor)
>>>> {
>>>>   assemble(A, form.form(), form.coefficients(), 0, 0, 0);
>>>> }
>>>>
>>>> reset_tensor is not used. Is this wanted? The initGlobalTensor takes
>>>> 90%
>>>> of my runtime (I think it wasn't like this before the last release).
>>>>
>>>
>>> That's exactly why it's there. If the matrix has already been
>>> initialised, then there is not need to re-initialise. This will
>>> typically be useful for nonlinear problems.
>>
>> Sorry my question was probably ill-posed. I mean that the matrix is
>> re-initialized even if I set the flag reset_tensor to false because it
>> is
>> not passed to the inner assemble().
>>
>
> Try it out now and let me know if it works.

I guess it needs to be inserted also in the call:

void Assembler::assemble(GenericTensor& A, const Form& form, bool
reset_tensor)

Thanks,
Alessio


> Garth
>
>> Alessio
>>
>>
>>> If initialising the matrix is requiring 90% of your runtime, you're
>>> probably doing something strange. Initialising the matrix does now take
>>> more time, with the trade-off that assembly is now much faster and less
>>> memory is required.
>>>
>>> Garth
>>>
>>>
>>>> Alessio
>>>>
>>>>
>>>>> Garth
>>>>>
>>>>>> Alessio
>>>>>>
>>>>>>
>>>>>>> One or more new changesets pushed to the primary dolfin repository.
>>>>>>> A short summary of the last three changesets is included below.
>>>>>>>
>>>>>>> changeset:   3358:fe232f61370532bd51b0f77df3fdd6aeaf5aa82c
>>>>>>> tag:         tip
>>>>>>> user:        "Anders Logg <logg at simula.no>"
>>>>>>> date:        Tue Aug 28 14:50:21 2007 +0200
>>>>>>> files:       src/kernel/fem/Assembler.cpp
>>>>>>> src/kernel/fem/assemble.cpp
>>>>>>> src/kernel/fem/dolfin/Assembler.h src/kernel/pde/LinearPDE.cpp
>>>>>>> src/kernel/pde/NonlinearPDE.cpp
>>>>>>> description:
>>>>>>> Add mesh to constructor of Assembler. Reuse of dof maps requires
>>>>>>> that the mesh remains constant.
>>>>>>>
>>>>>>>
>>>>>>> changeset:   3357:4ad8f4efd1d6c87e9d6dc82d85e2d8f3e4facff0
>>>>>>> user:        "Anders Logg <logg at simula.no>"
>>>>>>> date:        Tue Aug 28 14:26:43 2007 +0200
>>>>>>> files:       src/kernel/function/UserFunction.cpp
>>>>>>> src/kernel/function/dolfin/Function.h
>>>>>>> src/kernel/function/dolfin/UserFunction.h
>>>>>>> description:
>>>>>>> Make rank() and dim() virtual (overloadable by user).
>>>>>>>
>>>>>>>
>>>>>>> changeset:   3356:24c98268945163d11ab3b45f08f6015c25b4546c
>>>>>>> user:        "Anders Logg <logg at simula.no>"
>>>>>>> date:        Tue Aug 28 13:38:22 2007 +0200
>>>>>>> files:       src/demo/pde/poisson/python/demo.py
>>>>>>> description:
>>>>>>> test
>>>>>>>
>>>>>>> ----------------------------------------------------------------------
>>>>>>> For more details, visit http://www.fenics.org/hg/dolfin
>>>>>>>
>>>>>>>
>>>>>> _______________________________________________
>>>>>> DOLFIN-dev mailing list
>>>>>> DOLFIN-dev at fenics.org
>>>>>> http://www.fenics.org/mailman/listinfo/dolfin-dev
>>>>>>
>>>>>
>>>>
>>>> _______________________________________________
>>>> DOLFIN-dev mailing list
>>>> DOLFIN-dev@xxxxxxxxxx
>>>> http://www.fenics.org/mailman/listinfo/dolfin-dev
>>>>
>>>
>>>
>>
>>
>>
>> _______________________________________________
>> DOLFIN-dev mailing list
>> DOLFIN-dev@xxxxxxxxxx
>> http://www.fenics.org/mailman/listinfo/dolfin-dev
>>
>
>
>