Thank you very much, Prof.Bangerth. I have solved this problem using the first way you mentioned. 在 2019年8月18日星期日 UTC+8上午5:33:26，Wolfgang Bangerth写道： > > On 8/16/19 9:41 PM, hua zhang wrote: > > I am studying goal oriented adaptivity for FSI problems and > code it > > in deal.ii. The system consists of velocity, pressure, > displacement(u,p,w), > > so I use block sparsity matrix for system matrix. When I solve the dual > > equation, I find what I need is a transpose matrix of system matrix, so > how > > should I do to get the transpose of block sparsity matrix? > > I have tried to get transpose of every block using > > transpose_operator, but I don't know how to combined these terms into a > new > > block sparsity matrix. > > In essence, you only need an object whose vmult function calls the Tvmult > function of your block sparse matrix. That can be achieved with only a few > lines of code. But I think it can be done in an even easier way: create a > LinearOperator object from your block matrix, and then apply transpose() > to it > to get the transpose linear operator. > > Best > W. > > -- > ------------------------------------------------------------------------ > Wolfgang Bangerth email: [email protected] > <javascript:> > www: http://www.math.colostate.edu/~bangerth/ > >
-- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to [email protected]. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/6fcded48-f0c7-4614-a832-a0ec080f8899%40googlegroups.com.
