<div dir="ltr">the problem is with M. what about M?</div><div class="gmail_extra"><br><div class="gmail_quote">2016-12-16 17:33 GMT+03:00 Eric Chamberland <span dir="ltr"><<a href="mailto:Eric.Chamberland@giref.ulaval.ca" target="_blank">Eric.Chamberland@giref.ulaval.ca</a>></span>:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Ok, the offending call made is:<br>
<br>
 ierr = (*conv)(mat,newtype,reuse,M);C<wbr>HKERRQ(ierr);<br>
<br>
with:<br>
print newtype<br>
$1 = (MatType) 0x7fdc7bf8414a "seqdense"<br>
print reuse<br>
$2 = MAT_REUSE_MATRIX<br>
<br>
<br>
and mat is:<br>
<br>
(gdb) print *mat<br>
<br>
$4 = {hdr = {classid = 1211214, bops = {{getcomm = 0x7fdc6adda05b <PetscObjectGetComm_Petsc>, view = 0x7fdc6b7205e8 <MatView>, destroy = 0x7fdc6b7234b7 <MatDestroy>, compose = 0x7fdc6adda98d <PetscObjectCompose_Petsc>, query = 0x7fdc6addadbd <PetscObjectQuery_Petsc>, composefunction = 0x7fdc6addb28e <PetscObjectComposeFunction_Pe<wbr>tsc>, queryfunction = 0x7fdc6addb75f <PetscObjectQueryFunction_Pets<wbr>c>}}, comm = 0x5054930, type = 0, flops = 0, time = 0, mem = 16584, memchildren = 0, id = 41, refct = 1, tag = 2147483601, qlist = 0x4c0fc60, olist = 0x0, class_name = 0x7fdc6cc195e0 "Mat", description = 0x7fdc6cc195e4 "Matrix", mansec = 0x7fdc6cc195e0 "Mat", type_name = 0x4c0f600 "seqaij", parent = 0x0, parentid = 0, name = 0x0, prefix = 0x0, tablevel = 0, cpp = 0x0, state = 4, int_idmax = 0, intstar_idmax = 0, intcomposedstate = 0x0, intstarcomposedstate = 0x0, intcomposeddata = 0x0, intstarcomposeddata = 0x0, real_idmax = 0, realstar_idmax = 0, realcomposedstate = 0x0, realstarcomposedstate = 0x0, realcomposeddata = 0x0, realstarcomposeddata = 0x0, scalar_idmax = 0, scalarstar_idmax = 0, scalarcomposedstate = 0x0, scalarstarcomposedstate = 0x0, scalarcomposeddata = 0x0, scalarstarcomposeddata = 0x0, fortran_func_pointers = 0x0, num_fortran_func_pointers = 0, fortrancallback = {0x0, 0x0}, num_fortrancallback = {0, 0}, python_context = 0x0, python_destroy = 0x0, noptionhandler = 0, optionhandler = {0x0, 0x0, 0x0, 0x0, 0x0}, optiondestroy = {0x0, 0x0, 0x0, 0x0, 0x0}, optionctx = {0x0, 0x0, 0x0, 0x0, 0x0}, precision = PETSC_PRECISION_DOUBLE, optionsprinted = PETSC_FALSE, options = 0x0}, ops = {{setvalues = 0x7fdc6b53198a <MatSetValues_SeqAIJ>, getrow = 0x7fdc6b544343 <MatGetRow_SeqAIJ>, restorerow = 0x7fdc6b54476f <MatRestoreRow_SeqAIJ>, mult = 0x7fdc6b5bc57c <MatMult_SeqAIJ_Inode>, multadd = 0x7fdc6b5bdbda <MatMultAdd_SeqAIJ_Inode>, multtranspose = 0x7fdc6b53cd04 <MatMultTranspose_SeqAIJ>, multtransposeadd = 0x7fdc6b53c574 <MatMultTransposeAdd_SeqAIJ>, solve = 0x0, solveadd = 0x0, solvetranspose = 0x0, solvetransposeadd = 0x0, lufactor = 0x7fdc6b56e475 <MatLUFactor_SeqAIJ>, choleskyfactor = 0x0, sor = 0x7fdc6b5d3913 <MatSOR_SeqAIJ_Inode>, transpose = 0x7fdc6b545be6 <MatTranspose_SeqAIJ>, getinfo = 0x7fdc6b5426c4 <MatGetInfo_SeqAIJ>, equal = 0x7fdc6b55e32e <MatEqual_SeqAIJ>, getdiagonal = 0x7fdc6b53bcc9 <MatGetDiagonal_SeqAIJ>, diagonalscale = 0x7fdc6b548097 <MatDiagonalScale_SeqAIJ>, norm = 0x7fdc6b544a61 <MatNorm_SeqAIJ>, assemblybegin = 0x0, assemblyend = 0x7fdc6b538f0a <MatAssemblyEnd_SeqAIJ>, setoption = 0x7fdc6b53b7b6 <MatSetOption_SeqAIJ>, zeroentries = 0x7fdc6b53a23a <MatZeroEntries_SeqAIJ>, zerorows = 0x7fdc6b542abf <MatZeroRows_SeqAIJ>, lufactorsymbolic = 0x0, lufactornumeric = 0x0, choleskyfactorsymbolic = 0x0, choleskyfactornumeric = 0x0, setup = 0x7fdc6b54e099 <MatSetUp_SeqAIJ>, ilufactorsymbolic = 0x0, iccfactorsymbolic = 0x0, getdiagonalblock = 0x0, placeholder_33 = 0x0, duplicate = 0x7fdc6b55cc4f <MatDuplicate_SeqAIJ>, forwardsolve = 0x0, backwardsolve = 0x0, ilufactor = 0x7fdc6b54adf8 <MatILUFactor_SeqAIJ>, iccfactor = 0x0, axpy = 0x7fdc6b54f265 <MatAXPY_SeqAIJ>, getsubmatrices = 0x7fdc6b54bfd8 <MatGetSubMatrices_SeqAIJ>, increaseoverlap = 0x7fdc6b54c432 <MatIncreaseOverlap_SeqAIJ>, getvalues = 0x7fdc6b532b7e <MatGetValues_SeqAIJ>, copy = 0x7fdc6b54dbd1 <MatCopy_SeqAIJ>, getrowmax = 0x7fdc6b5508bd <MatGetRowMax_SeqAIJ>, scale = 0x7fdc6b54b8d1 <MatScale_SeqAIJ>, shift = 0x7fdc6b5541b6 <MatShift_SeqAIJ>, diagonalset = 0x7fdc6b52e223 <MatDiagonalSet_SeqAIJ>, zerorowscolumns = 0x7fdc6b54373d <MatZeroRowsColumns_SeqAIJ>, setrandom = 0x7fdc6b553bdb <MatSetRandom_SeqAIJ>, getrowij = 0x7fdc6b52e980 <MatGetRowIJ_SeqAIJ>, restorerowij = 0x7fdc6b52f0c0 <MatRestoreRowIJ_SeqAIJ>, getcolumnij = 0x7fdc6b52f642 <MatGetColumnIJ_SeqAIJ>, restorecolumnij = 0x7fdc6b52ffd6 <MatRestoreColumnIJ_SeqAIJ>, fdcoloringcreate = 0x7fdc6b58f28d <MatFDColoringCreate_SeqXAIJ>, coloringpatch = 0x0, setunfactored = 0x0, permute = 0x7fdc6b54cdc6 <MatPermute_SeqAIJ>, setvaluesblocked = 0x0, getsubmatrix = 0x0, destroy = 0x7fdc6b53a623 <MatDestroy_SeqAIJ>, view = 0x7fdc6b538922 <MatView_SeqAIJ>, convertfrom = 0x0, matmatmult = 0x7fdc6b5e8922 <MatMatMatMult_SeqAIJ_SeqAIJ_S<wbr>eqAIJ>, matmatmultsymbolic = 0x7fdc6b5e917a <MatMatMatMultSymbolic_SeqAIJ_<wbr>SeqAIJ_SeqAIJ>, matmatmultnumeric = 0x7fdc6b5e9a2f <MatMatMatMultNumeric_SeqAIJ_S<wbr>eqAIJ_SeqAIJ>, setlocaltoglobalmapping = 0x0, setvalueslocal = 0x0, zerorowslocal = 0x0, getrowmaxabs = 0x7fdc6b550244 <MatGetRowMaxAbs_SeqAIJ>, getrowminabs = 0x7fdc6b550fca <MatGetRowMinAbs_SeqAIJ>, convert = 0x0, dummy72 = 0x0, placeholder_73 = 0x0, setvaluesadifor = 0x0, fdcoloringapply = 0x7fdc6b6a0408 <MatFDColoringApply_AIJ>, setfromoptions = 0x0, multconstrained = 0x0, multtransposeconstrained = 0x0, findzerodiagonals = 0x7fdc6b52d703 <MatFindZeroDiagonals_SeqAIJ>, mults = 0x0, solves = 0x0, getinertia = 0x0, load = 0x7fdc6b55d1e7 <MatLoad_SeqAIJ>, issymmetric = 0x7fdc6b5479d3 <MatIsSymmetric_SeqAIJ>, ishermitian = 0x7fdc6b547d35 <MatIsHermitian_SeqAIJ>, isstructurallysymmetric = 0x0, setvaluesblockedlocal = 0x0, getvecs = 0x0, matmult = 0x7fdc6b595db4 <MatMatMult_SeqAIJ_SeqAIJ>, matmultsymbolic = 0x7fdc6b59e054 <MatMatMultSymbolic_SeqAIJ_Seq<wbr>AIJ>, matmultnumeric = 0x7fdc6b597e13 <MatMatMultNumeric_SeqAIJ_SeqA<wbr>IJ>, ptap = 0x7fdc6b5ab519 <MatPtAP_SeqAIJ_SeqAIJ>, ptapsymbolic = 0x7fdc6b5aeda4 <MatPtAPSymbolic_SeqAIJ_SeqAIJ<wbr>_DenseAxpy>, ptapnumeric = 0x7fdc6b5af661 <MatPtAPNumeric_SeqAIJ_SeqAIJ><wbr>, mattransposemult = 0x7fdc6b59f105 <MatMatTransposeMult_SeqAIJ_Se<wbr>qAIJ>, mattransposemultsymbolic = 0x7fdc6b59fe65 <MatMatTransposeMultSymbolic_S<wbr>eqAIJ_SeqAIJ>, mattransposemultnumeric = 0x7fdc6b5a0e80 <MatMatTransposeMultNumeric_Se<wbr>qAIJ_SeqAIJ>, placeholder_98 = 0x0, placeholder_99 = 0x0, placeholder_100 = 0x0, placeholder_101 = 0x0, conjugate = 0x7fdc6b54ff61 <MatConjugate_SeqAIJ>, placeholder_103 = 0x0, setvaluesrow = 0x7fdc6b531318 <MatSetValuesRow_SeqAIJ>, realpart = 0x7fdc6b539aee <MatRealPart_SeqAIJ>, imaginarypart = 0x7fdc6b539e9d <MatImaginaryPart_SeqAIJ>, getrowuppertriangular = 0x0, restorerowuppertriangular = 0x0, matsolve = 0x7fdc6b5704d4 <MatMatSolve_SeqAIJ>, placeholder_110 = 0x0, getrowmin = 0x7fdc6b5517ae <MatGetRowMin_SeqAIJ>, getcolumnvector = 0x0, missingdiagonal = 0x7fdc6b53f8ea <MatMissingDiagonal_SeqAIJ>, getseqnonzerostructure = 0x0, create = 0x0, getghosts = 0x0, getlocalsubmatrix = 0x0, restorelocalsubmatrix = 0x0, multdiagonalblock = 0x7fdc6b5e1597 <MatMultDiagonalBlock_SeqAIJ_I<wbr>node>, hermitiantranspose = 0x0, multhermitiantranspose = 0x0, multhermitiantransposeadd = 0x0, getmultiprocblock = 0x7fdc6b54a83c <MatGetMultiProcBlock_SeqAIJ>, findnonzerorows = 0x7fdc6b52daec <MatFindNonzeroRows_SeqAIJ>, getcolumnnorms = 0x7fdc6b52c3d1 <MatGetColumnNorms_SeqAIJ>, invertblockdiagonal = 0x7fdc6b551ebb <MatInvertBlockDiagonal_SeqAIJ<wbr>>, placeholder_127 = 0x0, getsubmatricesmpi = 0x0, setvaluesbatch = 0x0, transposematmult = 0x7fdc6b5a1c5a <MatTransposeMatMult_SeqAIJ_Se<wbr>qAIJ>, transposematmultsymbolic = 0x7fdc6b5a2498 <MatTransposeMatMultSymbolic_S<wbr>eqAIJ_SeqAIJ>, transposematmultnumeric = 0x7fdc6b5a2a65 <MatTransposeMatMultNumeric_Se<wbr>qAIJ_SeqAIJ>, transposecoloringcreate = 0x7fdc6b5a6ee6 <MatTransposeColoringCreate_Se<wbr>qAIJ>, transcoloringapplysptoden = 0x7fdc6b5a6094 <MatTransColoringApplySpToDen_<wbr>SeqAIJ>, transcoloringapplydentosp = 0x7fdc6b5a663b <MatTransColoringApplyDenToSp_<wbr>SeqAIJ>, rart = 0x7fdc6b5b4b9e <MatRARt_SeqAIJ_SeqAIJ>, rartsymbolic = 0x7fdc6b5b4258 <MatRARtSymbolic_SeqAIJ_SeqAIJ<wbr>>, rartnumeric = 0x7fdc6b5b47b2 <MatRARtNumeric_SeqAIJ_SeqAIJ><wbr>, setblocksizes = 0x0, aypx = 0x0, residual = 0x0, fdcoloringsetup = 0x7fdc6b5907a3 <MatFDColoringSetUp_SeqXAIJ>, findoffblockdiagonalentries = 0x7fdc6b52cb21 <MatFindOffBlockDiagonalEntrie<wbr>s_SeqAIJ>, creatempimatconcatenateseqmat = 0x7fdc6b5602f6 <MatCreateMPIMatConcatenateSeq<wbr>Mat_SeqAIJ>}}, rmap = 0x4c0e120, cmap = 0x4c0e7a0, data = 0x4c0ee20, factortype = MAT_FACTOR_NONE, assembled = PETSC_TRUE, was_assembled = PETSC_TRUE, num_ass = 2, nonzerostate = 1, info = {block_size = 0, nz_allocated = 0, nz_used = 0, nz_unneeded = 0, memory = 0, assemblies = 0, mallocs = 0, fill_ratio_given = 0, fill_ratio_needed = 0, factor_mallocs = 0}, insertmode = NOT_SET_VALUES, stash = {nmax = 0, umax = 0, oldnmax = 0, n = 0, bs = 0, reallocs = 0, space_head = 0x0, space = 0x0, ScatterBegin = 0x0, ScatterGetMesg = 0x0, ScatterEnd = 0x0, ScatterDestroy = 0x0, comm = 0x0, size = 0, rank = 0, tag1 = 0, tag2 = 0, send_waits = 0x0, recv_waits = 0x0, send_status = 0x0, nsends = 0, nrecvs = 0, svalues = 0x0, sindices = 0x0, rvalues = 0x0, rindices = 0x0, nprocessed = 0, flg_v = 0x0, reproduce = PETSC_FALSE, reproduce_count = 0, subset_off_proc = PETSC_FALSE, use_status = PETSC_FALSE, nsendranks = 0, nrecvranks = 0, sendranks = 0x0, recvranks = 0x0, sendhdr = 0x0, recvhdr = 0x0, sendframes = 0x0, recvframes = 0x0, recvframe_active = 0x0, recvframe_i = 0, recvframe_count = 0, recvcount = 0, some_indices = 0x0, some_statuses = 0x0, some_count = 0, some_i = 0, sendreqs = 0x0, recvreqs = 0x0, segsendblocks = 0x0, segrecvframe = 0x0, segrecvblocks = 0x0, blocktype = 0x0, blocktype_size = 0, insertmode = 0x0}, bstash = {nmax = 0, umax = 0, oldnmax = 0, n = 0, bs = 0, reallocs = 0, space_head = 0x0, space = 0x0, ScatterBegin = 0x0, ScatterGetMesg = 0x0, ScatterEnd = 0x0, ScatterDestroy = 0x0, comm = 0x0, size = 0, rank = 0, tag1 = 0, tag2 = 0, send_waits = 0x0, recv_waits = 0x0, send_status = 0x0, nsends = 0, nrecvs = 0, svalues = 0x0, sindices = 0x0, rvalues = 0x0, rindices = 0x0, nprocessed = 0, flg_v = 0x0, reproduce = PETSC_FALSE, reproduce_count = 0, subset_off_proc = PETSC_FALSE, use_status = PETSC_FALSE, nsendranks = 0, nrecvranks = 0, sendranks = 0x0, recvranks = 0x0, sendhdr = 0x0, recvhdr = 0x0, sendframes = 0x0, recvframes = 0x0, recvframe_active = 0x0, recvframe_i = 0, recvframe_count = 0, recvcount = 0, some_indices = 0x0, some_statuses = 0x0, some_count = 0, some_i = 0, sendreqs = 0x0, recvreqs = 0x0, segsendblocks = 0x0, segrecvframe = 0x0, segrecvblocks = 0x0, blocktype = 0x0, blocktype_size = 0, insertmode = 0x0}, nullsp = 0x0, transnullsp = 0x0, nearnullsp = 0x0, congruentlayouts = -1, preallocated = PETSC_TRUE, stencil = {dim = 0, dims = {0, 0, 0, 0}, starts = {0, 0, 0, 0}, noc = PETSC_FALSE}, symmetric = PETSC_FALSE, hermitian = PETSC_FALSE, structurally_symmetric = PETSC_FALSE, spd = PETSC_FALSE, symmetric_set = PETSC_FALSE, hermitian_set = PETSC_FALSE, structurally_symmetric_set = PETSC_FALSE, spd_set = PETSC_FALSE, symmetric_eternal = PETSC_FALSE, nooffprocentries = PETSC_FALSE, nooffproczerorows = PETSC_TRUE, subsetoffprocentries = PETSC_FALSE, submat_singleis = PETSC_FALSE, spptr = 0x0, solvertype = 0x0, checksymmetryonassembly = PETSC_FALSE, checknullspaceonassembly = PETSC_FALSE, checksymmetrytol = 0, redundant = 0x0, erroriffailure = PETSC_FALSE, factorerrortype = MAT_FACTOR_NOERROR, factorerror_zeropivot_value = 0, factorerror_zeropivot_row = 0}<span class="HOEnZb"><font color="#888888"><br>
<br>
Eric<br>
<br>
</font></span></blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature">Stefano</div>
</div>