> I have tried compiling the latest svn trunk (r19474) which includes > Martin's changes to > lac/include/lac/constraint_matrix.templates.h > The errors reported in my previous post do not get resolved. > > I have also tried making the change suggested by Wolfgang to r19474, > vec(constraint_line->entries[q].first) > += (vec(constraint_line->line) * > static_cast<typename VectorType::value_type> > (constraint_line->entries[q].second)); > This too does not resolve the error.
I have resolved the places in lac/ where we couldn't compile with complex PETSc installations (and these changes are now in svn), but there are many more places in deal.II. A lot of those don't compile because that just wouldn't make sense; examples are - GridRefinement::refine_and_coarsen* can be called with a PETSc vector, but since they expect refinement indicators it clearly doesn't make any sense if these indicators are complex-valued - the DoF iterators have functions get_dof_values that can take all kinds of vectors as inputs and are instantiated for real and double output vectors, allowing for example to extract values from a double vector into a float vector. But this fails if we try, for example, to extract double values from a complex PETSc vector. In some of these cases it's not quite clear what to do, or it will be a lot of work. Any suggestions? W. ------------------------------------------------------------------------- Wolfgang Bangerth email: [email protected] www: http://www.math.tamu.edu/~bangerth/ _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
