I am running out of memory when providing the Jacobian matrix of one of my components in my optimization problem (using SNOPT with pyOptSparse in OpenMDAO 0.13). Would it be possible to supply the Jacobian of one or more components in a sparse matrix format?

Thank you Pieter

asked 02 Sep '15, 16:07

pieterg's gravatar image


Without some modifications this isn't possible, as far as I know. However, there is code in there to zero pad any missing derivatives. You could probably adjust code to not do that. You can see this code here: https://github.com/OpenMDAO/OpenMDAO-Framework/blob/dev/openmdao.main/src/openmdao/main/linearsolver.py#L180 Specifically note lines 183 and 187. This is where the zero padding happens, and where you need to remove it to get more sparsity.

We did implement this sparsity in the alpha version of OpenMDAO though. So the features you need are available there. You could look at that wrapper code for hints as to how to make the older one work. The alpha code is here: https://github.com/OpenMDAO/OpenMDAO/blob/master/openmdao/core/problem.py#L801


answered 02 Sep '15, 16:19

justingray's gravatar image

justingray ♦♦

edited 02 Sep '15, 17:09

Follow this question

By Email:

Once you sign in you will be able to subscribe for any updates here



Answers and Comments

Markdown Basics

  • *italic* or _italic_
  • **bold** or __bold__
  • link:[text](http://url.com/ "title")
  • image?![alt text](/path/img.jpg "title")
  • numbered list: 1. Foo 2. Bar
  • to add a line break simply add two spaces to where you would like the new line to be.
  • basic HTML tags are also supported



Asked: 02 Sep '15, 16:07

Seen: 1,250 times

Last updated: 02 Sep '15, 17:09

powered by OSQA