
<html>

<head>

<meta http-equiv="Content-Type" content="text/html; charset=us-ascii">

<style type="text/css" style="display:none;"><!-- P {margin-top:0;margin-bottom:0;} --></style>

</head>

<body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;">

<div id="divtagdefaultwrapper" style="font-size: 12pt; color: rgb(0, 0, 0); background-color: rgb(255, 255, 255); font-family: Calibri, Arial, Helvetica, sans-serif; position: static; z-index: auto;">

Barry, Mark,<br>

<br>

Apologies for taking a while to respond.<br>

<div>

<div><br>

</div>

<blockquote type="cite">Are you saying it works for a while but fails when the problem is large, or that it never works with fieldsplit_1? And how many processors are you using?</blockquote>

</div>

<br>

<div>All of these are serial runs. The preconditioner always fails for large problems; but I also found an example where it fails for a fairly small matrix. fieldsplit_1 is the one that fails, although fieldpslit_0<span style="font-size: 12pt;"> </span><span style="font-size: 12pt;">always</span><span style="font-size: 12pt;"> works. </span></div>

<div><br>

</div>

<div>

<div>

<blockquote type="cite">This should be fixed in the master branch </blockquote>

</div>

<br>

</div>

<div>I tried the developmental version (from today), but the same bug persists.</div>

<div><br>

</div>

<div><br>

</div>

<div><br>

</div>

<div>I have some code that reproduces the example. The index sets that I use are non-overlapping; I also checked with matlab that each fieldsplit submatrix is non-singular.</div>

<div>The issue is with fieldsplit_TA: if I do not use -fieldsplit_TA-pc_type_gamg, then the solver converges. Otherwise, it crashes with the same message.</div>

<div><br>

</div>

<div><span style="font-size: 12pt;">Artur</span></div>

<div><br>

</div>

<div><br>

</div>

</div>

</body>

</html>