Some strange results running in PETSc.

2009-04-29 Thread (Rebecca) Xuefei YUAN
Hi, I am running some codes and stores the solution in the text file. However, I found that some results are wired in the sense that some processors are "eating" my (i,j) index and the corresponding solution. For example, the solution at time step =235 on processor 6 is right, but at time s

Convex optimization with linear constraint?

2009-04-29 Thread David Fuentes
Liu, Unless I'm missing something, I don't think you will directly find what you are looking for. You will prob have to solve A * vec_x = vec_b directly in your FormObjective function then use an adjoint method or something to compute your gradient directly in your FormGradient routine.