
<html>


<head>


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


</head>


<body>


<div dir="ltr">


<div class="gmail_quote">


<div dir="ltr">Markus :<br>


</div>


<div>Can you provide a concrete example of your application?</div>


<div>Is your model built on a network? Is so, give us a simple example of such network and the math systems over it.</div>


<div>This would help us understand your request.</div>


<div>Hong</div>


<div><br>


</div>


<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">


Dear PETSc users and developers,<br>


<br>


in particular those experienced with the relatively new DMNetwork object,<br>


I would like to get some advice on wether it makes sense for my application to be built using this PETSc feature<br>


or if I am equally well served if I use plain Vec and Mat objects.<br>


<br>


In the latter case, you might nevertheless have some good advice for a novice PETSc user<br>


or you might know about something that helps me to come up with a well architected solution.<br>


<br>


So the application context is closely linked to LTI state-space models (and in particular two-port / n-port network models and their interconnections):<br>


x,t = A x + B u<br>


y   = C x + D u<br>


<br>


More specifically, input ‘u' and output ‘y' are vectors (of same length for all components).<br>


Different components will have different dimension of state ‘x' (and hence also ‘A', ‘B', ‘C').<br>


<br>


These component models then have to be interconnected according to a given topology:<br>


Some pair of outputs of model 1  feeds into  the corresponding pair of inputs of model 2  (and also the other way round by symmetry), etc.<br>


<br>


After most (or even all) of the original inputs 'u_i' / outputs 'y_i' have been eliminated (based on the given interconnection structure amongst components),<br>


it will be necessary to use an iterative eigenvalue solver to obtain eigenvectors for some interesting part of the spectrum.<br>


<br>


The models will probably not be "very large" in the foreseeable future but this still doesn’t make e.g. MATLAB’s control toolbox an option.<br>


<br>


I have seen the presentation by Hong Zhang (1) at this year’s user meeting and I have looked at the paper “Scalable Multiphysics Network Simulation Using PETSc DMNetwork”.<br>


The use cases and concept of network presented therein were slightly different from this one of interconnected multi-ports in state-space formulation.<br>


<br>


Thanks you very much for your advice and reading this post.<br>


<br>


Best,<br>


Markus</blockquote>


</div>


</div>


</body>


</html>