<div dir="ltr"><div dir="ltr"><br></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sat, 2 Feb 2019 at 02:12, Sajid Ali <<a href="mailto:sajidsyed2021@u.northwestern.edu">sajidsyed2021@u.northwestern.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr">The vector is essentially snapshots in time of a data array. I should probably store this as a 2D dense matrix of dimensions (dim_x*dim_y) * dim_z. Now I can pick one column at a time and use it for my TS Jacobian. Apologies for being a little unclear. <br clear="all"><div><div><div><div><br></div></div></div></div></div></blockquote><div><br></div><div>m, n = dim_x*dim_y, dim_z</div><div>A = PETSc.Mat.createDense([m,n], array=your_array)</div><div><br></div></div>-- <br><div dir="ltr" class="gmail_signature">Lisandro Dalcin<br>============<br>Research Scientist<br>Computer, Electrical and Mathematical Sciences & Engineering (CEMSE)<br>Extreme Computing Research Center (ECRC)<br>King Abdullah University of Science and Technology (KAUST)<br><a href="http://ecrc.kaust.edu.sa/" target="_blank">http://ecrc.kaust.edu.sa/</a><br><br>4700 King Abdullah University of Science and Technology<br>al-Khawarizmi Bldg (Bldg 1), Office # 0109<br>Thuwal 23955-6900, Kingdom of Saudi Arabia<br><a href="http://www.kaust.edu.sa" target="_blank">http://www.kaust.edu.sa</a><br><br>Office Phone: +966 12 808-0459</div></div>