source: issm/oecreview/Archive/16133-16554/ISSM-16190-16191.diff@ 16556

../trunk-jpl/src/m/solvers/issmsolver.py

@@ -1 @@
-from pairoptions import *
-def asmoptions(*args):
-        #ISSMSOLVE - return issm solver options
-        #
-        #   Usage:
-        #      options=issmsolver;
-        
-        #retrieve options provided in varargin
-        arguments=pairoptions(*args) 
-        
-        options=[['toolkit','issm'],['mat_type','mpidense'],['vec_type','mpi'],['solver_type','mumps']];
-
-        #now, go through our arguments, and write over default options.
-        for i in range(len(arguments.list)):
-                arg1=arguments.list[i][0]
-                arg2=arguments.list[i][1]
-                found=0;
-                for j in range(len(options)):
-                        joption=options[j][0]
-                        if joption==arg1:
-                                joption[1]=arg2;
-                                options[j]=joption;
-                                found=1;
-                                break
-                if not found:
-                        #this option did not exist, add it: 
-                        options.append([arg1,arg2])
-        return options

../trunk-jpl/src/m/solvers/issmsolver.m

@@ -1 @@
-function issmoptions=issmsolver(varargin)
-%ISSMSOLVER - 
-%
-%   Usage:
-%      options=issmsolver;
-
-%retrieve options provided in varargin
-options=pairoptions(varargin{:});
-issmoptions=struct();
-
-%default issmoptions options
-issmoptions.toolkit='issm';
-issmoptions.mat_type=getfieldvalue(options,'mat_type','mpidense');
-issmoptions.vec_type=getfieldvalue(options,'vec_type','mpi');
-issmoptions.solver_type=getfieldvalue(options,'solver_type','mumps');

../trunk-jpl/src/m/solvers/issmmumpssolver.py

@@ +1 @@
+from pairoptions import *
+def issmmumpssolver(*args):
+        #ISSMSOLVE - return issm solver options
+        #
+        #   Usage:
+        #      options=issmsolver;
+        
+        #retrieve options provided in varargin
+        arguments=pairoptions(*args) 
+        
+        options=[['toolkit','issm'],['mat_type','mpidense'],['vec_type','mpi'],['solver_type','mumps']];
+
+        #now, go through our arguments, and write over default options.
+        for i in range(len(arguments.list)):
+                arg1=arguments.list[i][0]
+                arg2=arguments.list[i][1]
+                found=0;
+                for j in range(len(options)):
+                        joption=options[j][0]
+                        if joption==arg1:
+                                joption[1]=arg2;
+                                options[j]=joption;
+                                found=1;
+                                break
+                if not found:
+                        #this option did not exist, add it: 
+                        options.append([arg1,arg2])
+        return options

../trunk-jpl/src/m/solvers/issmmumpssolver.m

@@ +1 @@
+function issmoptions=issmmumpssolver(varargin)
+%ISSMSOLVER - 
+%
+%   Usage:
+%      options=issmsolver;
+
+%retrieve options provided in varargin
+options=pairoptions(varargin{:});
+issmoptions=struct();
+
+%default issmoptions options
+issmoptions.toolkit='issm';
+issmoptions.mat_type=getfieldvalue(options,'mat_type','mpidense');
+issmoptions.vec_type=getfieldvalue(options,'vec_type','mpi');
+issmoptions.solver_type=getfieldvalue(options,'solver_type','mumps');

../trunk-jpl/src/m/solvers/issmgslsolver.py

@@ +1 @@
+from pairoptions import *
+def issmgslsolver(*args):
+        #ISSMSOLVE - return issm solver options
+        #
+        #   Usage:
+        #      options=issmsolver;
+        
+        #retrieve options provided in varargin
+        arguments=pairoptions(*args) 
+        
+        options=[['toolkit','issm'],['mat_type','dense'],['vec_type','seq'],['solver_type','gsl']];
+
+        #now, go through our arguments, and write over default options.
+        for i in range(len(arguments.list)):
+                arg1=arguments.list[i][0]
+                arg2=arguments.list[i][1]
+                found=0;
+                for j in range(len(options)):
+                        joption=options[j][0]
+                        if joption==arg1:
+                                joption[1]=arg2;
+                                options[j]=joption;
+                                found=1;
+                                break
+                if not found:
+                        #this option did not exist, add it: 
+                        options.append([arg1,arg2])
+        return options

../trunk-jpl/src/m/solvers/issmgslsolver.m

@@ +1 @@
+function issmoptions=issmgslsolver(varargin)
+%ISSMSOLVER - 
+%
+%   Usage:
+%      options=issmsolver;
+
+%retrieve options provided in varargin
+options=pairoptions(varargin{:});
+issmoptions=struct();
+
+%default issmoptions options
+issmoptions.toolkit='issm';
+issmoptions.mat_type=getfieldvalue(options,'mat_type','dense');
+issmoptions.vec_type=getfieldvalue(options,'vec_type','seq');
+issmoptions.solver_type=getfieldvalue(options,'solver_type','gsl');