Re: [OPM] Porsol: IncompFlowSolverHybrid - decrypting the procedure

2015-03-16 Thread Jørgen Kvalsvik
On 03/16/2015 07:07 PM, Bård Skaflestad wrote: > On 16/03/15 18:43, Jørgen Kvalsvik wrote: >> Right, but isn't this also, I suppose implicitly, reforming the >> unchanged sparsity structure? > > "It depends". > > Specifically, it depends specifically on what you mean by the term > "sparsity struc

Re: [OPM] Porsol: IncompFlowSolverHybrid - decrypting the procedure

2015-03-16 Thread Bård Skaflestad
On 16/03/15 18:43, Jørgen Kvalsvik wrote: On 03/16/2015 06:34 PM, Bård Skaflestad wrote: On 16/03/15 18:28, Jørgen Kvalsvik wrote: On 03/16/2015 05:07 PM, Bård Skaflestad wrote: We currently implement the second approach. The matrix sparsity (connection) structure does not change between call

Re: [OPM] Porsol: IncompFlowSolverHybrid - decrypting the procedure

2015-03-16 Thread Jørgen Kvalsvik
On 03/16/2015 06:34 PM, Bård Skaflestad wrote: > On 16/03/15 18:28, Jørgen Kvalsvik wrote: >> On 03/16/2015 05:07 PM, Bård Skaflestad wrote: >>> We currently implement the second approach. The matrix sparsity >>> (connection) structure does not change between calls to solve(), >>> although the num

Re: [OPM] Porsol: IncompFlowSolverHybrid - decrypting the procedure

2015-03-16 Thread Bård Skaflestad
On 16/03/15 16:51, Jørgen Kvalsvik wrote: On 03/16/2015 04:43 PM, Bård Skaflestad wrote: As for the overall construction approach, it basically boils down to which is more expensive: * Having two copies of the matrix in memory and re-forming the (unchanging) sparsity structure on ea

Re: [OPM] Porsol: IncompFlowSolverHybrid - decrypting the procedure

2015-03-16 Thread Jørgen Kvalsvik
On 03/16/2015 04:43 PM, Bård Skaflestad wrote: > On 16/03/15 16:11, Jørgen Kvalsvik wrote: > in the periodic case, it uses plain assignment S_[i][j] = x, >> whereas in all the other cases it uses +=. > > Do you mean *Dirichlet* case rather than periodic in the above statement? > > I ported it to

[OPM] Fwd: Re: Porsol: IncompFlowSolverHybrid - decrypting the procedure

2015-03-16 Thread Bård Skaflestad
Sorry, I hit "reply" rather than "reply all". Bård Forwarded Message Subject: Re: [OPM] Porsol: IncompFlowSolverHybrid - decrypting the procedure Date: Mon, 16 Mar 2015 16:43:23 +0100 From: Bård Skaflestad To: Jørgen Kvalsvik On 16/03/15 16:11, Jørgen Kvalsvik wrote: in t

Re: [OPM] Porsol: IncompFlowSolverHybrid - decrypting the procedure

2015-03-16 Thread Jørgen Kvalsvik
On 03/16/2015 04:06 PM, Bård Skaflestad wrote: > On 16/03/15 15:51, Jørgen Kvalsvik wrote: >> On 03/16/2015 02:49 PM, Bård Skaflestad wrote: >>> On 26/01/15 13:00, Jørgen Kvalsvik wrote: In the periodic case, will some indices be encountered twice? >>> >>> I'm not sure I understand the questio

Re: [OPM] Porsol: IncompFlowSolverHybrid - decrypting the procedure

2015-03-16 Thread Bård Skaflestad
On 16/03/15 15:51, Jørgen Kvalsvik wrote: On 03/16/2015 02:49 PM, Bård Skaflestad wrote: On 26/01/15 13:00, Jørgen Kvalsvik wrote: In the periodic case, will some indices be encountered twice? I'm not sure I understand the question. All indices may, and typically will, be encountered multipl

Re: [OPM] Porsol: IncompFlowSolverHybrid - decrypting the procedure

2015-03-16 Thread Jørgen Kvalsvik
On 03/16/2015 02:49 PM, Bård Skaflestad wrote: > On 26/01/15 13:00, Jørgen Kvalsvik wrote: >> Hi, > > Hi, > > Sorry about the delay. I was just going through my old e-mail and > noticed that I didn't see an answer to this question. I'll provide one > for the mailing list record, even it much ti

Re: [OPM] Porsol: IncompFlowSolverHybrid - decrypting the procedure

2015-03-16 Thread Bård Skaflestad
On 26/01/15 13:00, Jørgen Kvalsvik wrote: Hi, Hi, Sorry about the delay. I was just going through my old e-mail and noticed that I didn't see an answer to this question. I'll provide one for the mailing list record, even it much time has passed since the question was posed. Can the Dir