MatSolveTranspose

[Hong Zhang]

Laslo,

MatSolveTranspose() should work for sequential BAIJ format.
The example ~petsc/src/mat/examples/tests/ex48.c tests it.
Would you please run this example and see if it works.

You may simplify your code and send it to us.
Then I'll test it to see where is the problem.

Hong

On Mon, 7 Aug 2006, Laslo Tibor Diosady wrote:

> Hi,
>
> I am trying to use MatSolveTranspose for a sequential BAIJ format, however
> whenever I call MatSolveTranspose I get a segmentation fault.
>
> The sequence of calls which I make is:
> MatILUFactorSymbolic
> MatLUFactorNumeric
>
> And I use the resulting matrix in calls to:
> MatSolve and MatSolveTranspose
>
> When I use MatSolve this works great. However, with MatSolveTranspose it
> appears petsc tries to allocate memory in a never ending loop until my
> machine runs out of memory and I get a seg fault.
>
> If I do a call to MatHasOperation with this matrix the result is
> Petsc_True, so in theory I think the call to MatSolveTranspose should
> work.
>
> Am I doing something wrong or is this a problem with MatSolveTranspose for
> seqBAIJ matrices?
>
> Any help would be greatly appreciated.
>
> Thanks,
>
> Laslo
>
>