<div dir="ltr"><div dir="ltr"><div>Hi,</div><div><br></div><div>I am trying to solve a Ax=b complex system. the vector b and "matrix" A are both binary and NOT created by PETSc. So I keep getting error messages that they are not correct format when I read the files with PetscViewBinaryOpen, after some digging it seems that one cannot just read a binary file that was created by another software. <br></div><div><br></div><div>How would I go on to solve this problem? <br></div><div><br></div><div>More info and trials:</div><div><br></div><div>"matrix" A consists of two files, one that contains row column index numbers and one that contains the non-zero values. So what I would have to do is multiply the last term in a+b with PETSC_i to get a real + imaginary vector A.</div><div><br></div><div>vector b is in binary, so what I have done so far (not sure if it works) is:<br></div><div><br></div><div>std::ifstream input("Vector_b.bin", std::ios::binary );<br>while (input.read(reinterpret_cast<char*>(&v), sizeof(float)))<br>     ierr    = VecSetValues(u,1,&iglobal,&v,INSERT_VALUES);CHKERRQ(ierr);</div><div><br></div><div>where v is a PetscScalar.</div><div><br></div><div>Once I am able to read both matrices I think I can figure out the solvers to solve the system.</div><div><br></div><div>All the best,</div><div>S<br></div></div></div>