<div dir="ltr"><div><div><div><div>Hi,<br><br></div>I am trying to use PCBDDC for my substructured problem. I have each subdomains stiffneesses and load vectors and also mapping of numbering to global dofs. I copied ex59 and replaced proper functions for my problem. Everything seems working but, it does not produce correct results! I checked kspconvergedreaseon and it is 3 and iteration count is 0 !. I started debugging and know I am suspicious about the RHS vector. <br><br> // assemble global matrix <br> ierr = ComputeMatrix(dd,&K);CHKERRQ(ierr);<br> // assemble BDDC rhs <br> ierr = MatGetVecs(K,&F,NULL);CHKERRQ(ierr);<br> ierr = VecZeroEntries(F);CHKERRQ(ierr);<br> ierr = VecSetValues(F,Mapping.rows(),Mapping.data(),b.data(),ADD_VALUES);CHKERRQ(ierr);<br><br> ierr = VecAssemblyBegin(F);CHKERRQ(ierr);<br> ierr = VecAssemblyEnd(F);CHKERRQ(ierr);<br><br><br></div>when I use matview and vecview, I see that K is matis and proc0 has 18x18 part of K, proc 1 has 12 part of K. 6 dofs are overlapping (dirichlet boundaries) so K is 24x24. However F seems that it has 12x1 at each processor. I guess sizes of K and F at each proc should be same and MatGetVecs should ensure same mapping for both K and F but it seems it does not ! <br><br></div>Same problem is solved in other ksp solvers, I wonder how to prepare RHS for PCBDDC. <br><br>thank you<br><br></div>PS: I am at the latest   master + stefano_zampini/pcbddc-primalfixes<br></div>