Search found 357 matches

by Juha
06 Feb 2018, 11:15
Forum: Elmer/Ice
Topic: Changes to parstokes Solver?
Replies: 6
Views: 9237

Re: Changes to parstokes Solver?

Hi The patch below fixes the case. There seems to have been a misunderstaning of the role of the variable initialization statement in Fortran, which is very different from c/c++. Fortran initializes the variable only once (not at each re-entry) and so the initialization also implies saving of the va...
by Juha
05 Feb 2018, 14:08
Forum: Bug reports
Topic: Passive condition
Replies: 3
Views: 4548

Re: Passive condition

Hi
To my best understaning your fix is good, and I applied it to the "devel" -branch. Thanks for the
fix and sorry 'bout the breakage.
BR, juha
by Juha
29 Aug 2014, 11:53
Forum: Bug reports
Topic: SuperLU
Replies: 2
Views: 4301

Re: SuperLU

Hi,

made the changes as you suggested (in the SVN). Can't test it atm though, maybe you
could give it a try?

Thanks, Juha
by Juha
28 Jun 2013, 13:35
Forum: ElmerSolver
Topic: Integral COnstraint for single Solver
Replies: 12
Views: 9784

Re: Integral COnstraint for single Solver

Hi, there is also a second system if you have constraints such that you can attach the new rows to DOFs of the original matrix, namely A % ConstraintMatrix. An example of using this to constraint values on a boundary to given function would be: p => Solver % Variable % Perm CM => AllocateMatrix() CM...
by Juha
28 Jun 2013, 12:35
Forum: ElmerSolver
Topic: Integral COnstraint for single Solver
Replies: 12
Views: 9784

Re: Integral COnstraint for single Solver

Hi, The AddMatrix system adds new global row(s) to the matrix that are shared by all partitions. You can think of them having globally indices from 1 to max_partitions(addmatrix%numberofrows-a%numberofrows). So, no, the partition AddMatrices are not independent, but share the row numering. If you ne...
by Juha
13 Mar 2013, 12:19
Forum: ElmerGUI
Topic: Output Length Exceeded Error
Replies: 8
Views: 8849

Re: Output Length Exceeded Error

Hi,
unfortunately the gfortran compiler packaged with elmer windows
binaries doesn't support allocatable characters at all, so the current
svn sources are unusable with current windows binary installation.
I'll try to update the gcc in our packages, but that'll take some
time....
Regards,Juha
by Juha
13 Mar 2013, 11:46
Forum: ElmerGUI
Topic: Output Length Exceeded Error
Replies: 8
Views: 8849

Re: Output Length Exceeded Error

Hi, i changed the sources in SVN such that allocatable character strings are now used instead of fixed length character strings. This particular problem shoud not occur after fetching the sources and recompiling, regardless of the MAX_STRING_LEN setting. The allocatable characters are somewhat newis...
by Juha
12 Mar 2013, 08:19
Forum: Software development
Topic: dBasisdx over Edges
Replies: 1
Views: 3512

Re: dBasisdx over Edges

Hi,

the derivatives from ElementInfo are with respect to global coordinates.
In your example the edge might be oriented along y-coordinate, for
example? If you need tangential derivatives, take a dot product with
edge direction vector?

Regards, Juha
by Juha
17 Jan 2013, 12:40
Forum: Bug reports
Topic: parallel AVWhitney: DefaultDirichletBCs IterSolve Failed
Replies: 14
Views: 14080

Re: parallel AVWhitney: DefaultDirichletBCs IterSolve Failed

Hi Eelis,

impossible to say anything w/o having the whole model. Have you
tried with "Linear System Preconditioning=None"?

Regards,Juha
by Juha
05 Dec 2012, 13:14
Forum: ElmerSolver
Topic: Convergence of Whitney-AV Solver
Replies: 3
Views: 3451

Re: Convergence of Whitney-AV Solver

Hi,
If you haven't done so already, you could try using BiCGStabL instead of TFQMR, also
adding "BiCGStabL Polynomial Degree=4" might help(?) GMRES with the
"Linear System Gmres Restart=n" is also an option.

Using ILU might actually make things worse.
Regards, Juha