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

mongodb - How to keep track of users making Stripe Payments -

i'm trying figure out way mean stack application keep track of users have paid grant them access portion of webpage. i've considered several options: stripe customer id, mongodb record, , html attribute can update. my mean stack keeps track of users jwt, , appears stripe assigns them own customer id isn't ideal. can done jwt opposed forced cutomer id? mongodb record. i'm thinking might best option. when new user has been created, i'll give attribute of haspaid = no. update record of customer when payment submitted. guess run script set unpaid each day? html element/attribute. don't know if possible; cool create key carried during html session after payment received. if person closers browser session closed? i'm looking guidance on 3 options determine if they're best solution. also, if has suggestions alternatives, i'm ears! thanks in advance. speaking generally, common approach second one: use attribute in data model indicates
Read more