How to construct petsc matrices? -


i'm using petsc4py , face problems. have number of small petsc dense matrices mij, , want construct them big matrix m this:

     [  m11  m12  m13  ] m =  |  m21  m22  m23  |   ,      [  m31  m32  m33  ]

a mcve code showing below, , i'm using python wrap of petsc, however, grammars similar. 

import numpy np petsc4py import petsc  msizes = (5, 8, 6) mij = []  # create sub-matrices mij in range(len(msizes)):     j in range(len(msizes)):         temp_m = petsc.mat().create(comm=petsc.comm_world)         temp_m.setsizes(((none, msizes[i]), (none, msizes[j])))         temp_m.settype('mpidense')         temp_m.setfromoptions()         temp_m.setup()         temp_m[:, :] = np.random.random_sample((msizes[i], msizes[j]))         temp_m.assemble()         mij.append(temp_m)  # have 4 sub-matrices.  # construct them big matrix m. m = petsc.mat().create(comm=petsc.comm_world) m.setsizes(((none, np.sum(msizes)), (none, np.sum(msizes)))) m.settype('mpidense') m.setfromoptions() m.setup() mlocations = np.insert(np.cumsum(msizes), 0, 0)    # mlocations = [0, msizes] in range(len(msizes)):     j in range(len(msizes)):         m[mlocations[i]:mlocations[i+1], mlocations[j]:mlocations[j+1]] = \              mij[i*len(msizes)+j][:, :] m.assemble()

which report such error messages:

traceback (most recent call last):   file "trymatconstuct.py", line 29, in <module>     mij[i*len(msizes)+j][:, :]   file "petsc/mat.pyx", line 227, in petsc4py.petsc.mat.__getitem__ (src/petsc4py.petsc.c:110477)   file "petsc/petscmat.pxi", line 997, in petsc4py.petsc.mat_getitem (src/petsc4py.petsc.c:30481)   file "petsc/petscmat.pxi", line 917, in petsc4py.petsc.matgetvalues (src/petsc4py.petsc.c:29242) petsc4py.petsc.error: error code 56 [1] matgetvalues() line 1818 in /home/zhangji/pycharmprojects/petsc-petsc-31a1859eaff6/src/mat/interface/matrix.c [1] matgetvalues_mpidense() line 154 in /home/zhangji/pycharmprojects/petsc-petsc-31a1859eaff6/src/mat/impls/dense/mpi/mpidense.c [1] no support operation object type [1] local values supported

i have use 2d array instead of petsc matrix short sub-matrix mij.

in petsc4py function petsc.mat().getdensearray(), i.e.

temp_m = mij[i*len(msizes)+j].getdensearray()

-

Comments

Post a Comment

Popular posts from this blog

javascript - Thinglink image not visible until browser resize -

Image
i trying repurpose pre-existing theme, includes drop down menu displays image. trying embed thinglink drop down menu's content section, not visible until resize browser window. here menu looks like: then when clicked, can see empty space image should show: then when browser resized, image appears: this code content of menu: <div class="ajax_accordion_content" style="display: none;"> <div class="report-detail"> <table border="0" cellpadding="0" cellspacing="0" width="100%" class="table_basic"> <tbody> <!-- thinglink image embedded --> <img style="max-width: 100%;" src="http://cdn.thinglink.me/api/image/832705713204625409/1024/10/scaletowidth#tl-832705713204625409;1043138249'" class="alwaysthinglink"/> <script async charset="utf-8" src="http://cdn.thinglink.me/jse/embed.j...
Read more

firebird - Error "invalid transaction handle (expecting explicit transaction start)" executing script from Delphi -

i'm developing update app delphi 10. i'm running windows 7 64x, firebird 2.5.1.26351 32x. execution order: checks current version(select statement firebird database) downloads update(via ftp) apply it(.exe files, , execute .sql) error code -901 pops out when try execute .sql files. know there incompatibility issue w7 64x , fdb 32x, but, can everything: backup , restore(gbak), db commands(all). tried fdb 64x instead, , not working. code: //components fdwaitcursor: tfdguixwaitcursor; fddriverlink: tfdphysfbdriverlink; fdconnection:tfdconnection; fdscript: tfdscript; function tthr_script.createcomponents:boolean; begin try fddriverlink := tfdphysfbdriverlink.create(fddriverlink); fdwaitcursor := tfdguixwaitcursor.create(fdwaitcursor); fdconnection := tfdconnection.create(fdconnection); fdconnection begin drivername := 'fb'; params.username := thrbanco.bduser; params.password := thrbanco.bdpass; params.database := ...
Read more

Sound is not coming out while implementing Text-to-speech in Android activity -

i using texttospeak in 1 of app. while adding speech string texttospeech object instance, doesn't produce sound. though not getting error in log still sound not coming out. want device speak value of textview gets set resultmsg.settext() method. public class resultactivity extends activity { @override protected void oncreate(bundle savedinstancestate){ ............ resultmsg.settext(resultmsg_str); texttospeech = new texttospeech(getapplicationcontext(),newtexttospeech.oninitlistener(){ @override public void oninit(int status) { system.out.println("txt speech status = "+status); system.out.println("txt speech error status = "+texttospeech.error); system.out.println("txt speech success status = "+texttospeech.success); if(status == texttospeech.success){ texttospeech.setlanguage(locale.english); ...
Read more