dolfin team mailing list archive

["Garth N. Wells" <gnw20@xxxxxxxxx>]

Anders Logg wrote:
> On Tue, Feb 23, 2010 at 11:44:44PM +0000, Garth N. Wells wrote:
>>
>> Anders Logg wrote:
>>> On Tue, Feb 23, 2010 at 10:38:38PM +0000, Garth N. Wells wrote:
>>>> Marie Rognes wrote:
>>>>> Am I the only one getting the following error with a fresh dolfin?
>>>>>
>>>> I've commented out the offending function. I don't see why it was added.
>>>>
>>>> Garth
>>> It's in the commit log. :-)
>>>
>> Where's the implementation? I don't see why it was added without an
>> implementation.
> 
> The .cpp file was missing. It's there now.
> 
>> I also think that this isn't a good idea from the la standpoint,
>> especially in parallel.
> 
> It works through ident() which works fine in parallel. 

Fine with PETSc, but not in general.

> The problem is
> getrow() that doesn't work in parallel. There might be other ways to
> figure out which rows are zero.
> 

get() functions, especially for matrices, should generally be avoided.

Garth

> The function can be removed when we get restrictions to work.
> 
> --
> Anders
> 
> 
>> Garth
>>
>>> The function ident_zeros inserts 1 in the diagonal in all rows where
>>> all entries are zero. This can be used to constrain inactive dofs.
>>>
>>> The best option for this would be to implement support for restricted
>>> function spaces (where dofs in some regions are inactive) but it's
>>> currently broken.
>>>
>>>>> meg@localhost ~ $ python
>>>>> Python 2.6.4 (r264:75706, Dec  7 2009, 23:01:27)
>>>>> [GCC 4.3.2] on linux2
>>>>> Type "help", "copyright", "credits" or "license" for more information.
>>>>>>>> import dolfin
>>>>> Traceback (most recent call last):
>>>>>  File "<stdin>", line 1, in <module>
>>>>>  File
>>>>> "/home/meg/local/lib64/python2.6/site-packages/dolfin/__init__.py", line
>>>>> 13, in <module>
>>>>>    import dolfin.cpp as cpp
>>>>>  File "/home/meg/local/lib64/python2.6/site-packages/dolfin/cpp.py",
>>>>> line 6, in <module>
>>>>>    import _cpp
>>>>> ImportError: /home/meg/local/lib/libdolfin.so.0: undefined symbol:
>>>>> _ZN6dolfin13GenericMatrix11ident_zerosEv
>>>> _______________________________________________
>>>> Mailing list: https://launchpad.net/~dolfin
>>>> Post to     : dolfin@xxxxxxxxxxxxxxxxxxx
>>>> Unsubscribe : https://launchpad.net/~dolfin
>>>> More help   : https://help.launchpad.net/ListHelp