xref: /trunk/main/sccomp/source/solver/solver.cxx (revision cf6516809c57e1bb0a940545cca99cdad54d4ce2)
1e4f93722SAndrew Rist /**************************************************************
2cdf0e10cSrcweir  *
3e4f93722SAndrew Rist  * Licensed to the Apache Software Foundation (ASF) under one
4e4f93722SAndrew Rist  * or more contributor license agreements.  See the NOTICE file
5e4f93722SAndrew Rist  * distributed with this work for additional information
6e4f93722SAndrew Rist  * regarding copyright ownership.  The ASF licenses this file
7e4f93722SAndrew Rist  * to you under the Apache License, Version 2.0 (the
8e4f93722SAndrew Rist  * "License"); you may not use this file except in compliance
9e4f93722SAndrew Rist  * with the License.  You may obtain a copy of the License at
10cdf0e10cSrcweir  *
11e4f93722SAndrew Rist  *   http://www.apache.org/licenses/LICENSE-2.0
12cdf0e10cSrcweir  *
13e4f93722SAndrew Rist  * Unless required by applicable law or agreed to in writing,
14e4f93722SAndrew Rist  * software distributed under the License is distributed on an
15e4f93722SAndrew Rist  * "AS IS" BASIS, WITHOUT WARRANTIES OR CONDITIONS OF ANY
16e4f93722SAndrew Rist  * KIND, either express or implied.  See the License for the
17e4f93722SAndrew Rist  * specific language governing permissions and limitations
18e4f93722SAndrew Rist  * under the License.
19cdf0e10cSrcweir  *
20e4f93722SAndrew Rist  *************************************************************/
21e4f93722SAndrew Rist 
22e4f93722SAndrew Rist 
23*1f56509eSdamjan #include <CoinMP.h>
24cdf0e10cSrcweir 
25cdf0e10cSrcweir #include "solver.hxx"
26cdf0e10cSrcweir #include "solver.hrc"
27cdf0e10cSrcweir 
28cdf0e10cSrcweir #include <com/sun/star/beans/XPropertySet.hpp>
29cdf0e10cSrcweir #include <com/sun/star/container/XIndexAccess.hpp>
30cdf0e10cSrcweir #include <com/sun/star/frame/XModel.hpp>
31cdf0e10cSrcweir #include <com/sun/star/lang/XMultiServiceFactory.hpp>
32cdf0e10cSrcweir #include <com/sun/star/sheet/XSpreadsheetDocument.hpp>
33cdf0e10cSrcweir #include <com/sun/star/sheet/XSpreadsheet.hpp>
34cdf0e10cSrcweir #include <com/sun/star/table/CellAddress.hpp>
35cdf0e10cSrcweir #include <com/sun/star/table/CellRangeAddress.hpp>
36cdf0e10cSrcweir #include <com/sun/star/text/XTextRange.hpp>
37cdf0e10cSrcweir 
38cdf0e10cSrcweir #include <rtl/math.hxx>
39cdf0e10cSrcweir #include <rtl/ustrbuf.hxx>
40cdf0e10cSrcweir #include <cppuhelper/factory.hxx>
41cdf0e10cSrcweir #include <vector>
42cdf0e10cSrcweir #include <hash_map>
43cdf0e10cSrcweir 
44cdf0e10cSrcweir #include <tools/resmgr.hxx>
45cdf0e10cSrcweir 
46cdf0e10cSrcweir using namespace com::sun::star;
47cdf0e10cSrcweir 
48cdf0e10cSrcweir using ::rtl::OUString;
49cdf0e10cSrcweir 
50cdf0e10cSrcweir #define C2U(constAsciiStr) (::rtl::OUString( RTL_CONSTASCII_USTRINGPARAM( constAsciiStr ) ))
51cdf0e10cSrcweir 
52cdf0e10cSrcweir #define STR_NONNEGATIVE   "NonNegative"
53cdf0e10cSrcweir #define STR_INTEGER       "Integer"
54cdf0e10cSrcweir #define STR_TIMEOUT       "Timeout"
55cdf0e10cSrcweir #define STR_EPSILONLEVEL  "EpsilonLevel"
56cdf0e10cSrcweir #define STR_LIMITBBDEPTH  "LimitBBDepth"
57b1f1da72SPedro Giffuni #define STR_NONLINEARTEST "NonLinearTest"
58cdf0e10cSrcweir 
59cdf0e10cSrcweir // -----------------------------------------------------------------------
60cdf0e10cSrcweir //  Resources from tools are used for translated strings
61cdf0e10cSrcweir 
62cdf0e10cSrcweir static ResMgr* pSolverResMgr = NULL;
63cdf0e10cSrcweir 
lcl_GetResourceString(sal_uInt32 nId)64cdf0e10cSrcweir OUString lcl_GetResourceString( sal_uInt32 nId )
65cdf0e10cSrcweir {
66cdf0e10cSrcweir     if (!pSolverResMgr)
67cdf0e10cSrcweir         pSolverResMgr = CREATEVERSIONRESMGR( solver );
68cdf0e10cSrcweir 
69cdf0e10cSrcweir     return String( ResId( nId, *pSolverResMgr ) );
70cdf0e10cSrcweir }
71cdf0e10cSrcweir 
72cdf0e10cSrcweir // -----------------------------------------------------------------------
73cdf0e10cSrcweir 
74cdf0e10cSrcweir namespace
75cdf0e10cSrcweir {
76cdf0e10cSrcweir     enum
77cdf0e10cSrcweir     {
78cdf0e10cSrcweir         PROP_NONNEGATIVE,
79cdf0e10cSrcweir         PROP_INTEGER,
80cdf0e10cSrcweir         PROP_TIMEOUT,
81cdf0e10cSrcweir         PROP_EPSILONLEVEL,
82b1f1da72SPedro Giffuni         PROP_LIMITBBDEPTH,
83b1f1da72SPedro Giffuni         PROP_NONLINEARTEST
84cdf0e10cSrcweir     };
85cdf0e10cSrcweir }
86cdf0e10cSrcweir 
87cdf0e10cSrcweir // -----------------------------------------------------------------------
88cdf0e10cSrcweir 
89cdf0e10cSrcweir // hash map for the coefficients of a dependent cell (objective or constraint)
90cdf0e10cSrcweir // The size of each vector is the number of columns (variable cells) plus one, first entry is initial value.
91cdf0e10cSrcweir 
92cdf0e10cSrcweir struct ScSolverCellHash
93cdf0e10cSrcweir {
operator ()ScSolverCellHash94cdf0e10cSrcweir     size_t operator()( const table::CellAddress& rAddress ) const
95cdf0e10cSrcweir     {
96cdf0e10cSrcweir         return ( rAddress.Sheet << 24 ) | ( rAddress.Column << 16 ) | rAddress.Row;
97cdf0e10cSrcweir     }
98cdf0e10cSrcweir };
99cdf0e10cSrcweir 
AddressEqual(const table::CellAddress & rAddr1,const table::CellAddress & rAddr2)100cdf0e10cSrcweir inline bool AddressEqual( const table::CellAddress& rAddr1, const table::CellAddress& rAddr2 )
101cdf0e10cSrcweir {
102cdf0e10cSrcweir     return rAddr1.Sheet == rAddr2.Sheet && rAddr1.Column == rAddr2.Column && rAddr1.Row == rAddr2.Row;
103cdf0e10cSrcweir }
104cdf0e10cSrcweir 
105cdf0e10cSrcweir struct ScSolverCellEqual
106cdf0e10cSrcweir {
operator ()ScSolverCellEqual107cdf0e10cSrcweir     bool operator()( const table::CellAddress& rAddr1, const table::CellAddress& rAddr2 ) const
108cdf0e10cSrcweir     {
109cdf0e10cSrcweir         return AddressEqual( rAddr1, rAddr2 );
110cdf0e10cSrcweir     }
111cdf0e10cSrcweir };
112cdf0e10cSrcweir 
113cdf0e10cSrcweir typedef std::hash_map< table::CellAddress, std::vector<double>, ScSolverCellHash, ScSolverCellEqual > ScSolverCellHashMap;
114cdf0e10cSrcweir 
115cdf0e10cSrcweir // -----------------------------------------------------------------------
116cdf0e10cSrcweir 
lcl_GetCell(const uno::Reference<sheet::XSpreadsheetDocument> & xDoc,const table::CellAddress & rPos)117cdf0e10cSrcweir uno::Reference<table::XCell> lcl_GetCell( const uno::Reference<sheet::XSpreadsheetDocument>& xDoc,
118cdf0e10cSrcweir                                           const table::CellAddress& rPos )
119cdf0e10cSrcweir {
120cdf0e10cSrcweir     uno::Reference<container::XIndexAccess> xSheets( xDoc->getSheets(), uno::UNO_QUERY );
121cdf0e10cSrcweir     uno::Reference<sheet::XSpreadsheet> xSheet( xSheets->getByIndex( rPos.Sheet ), uno::UNO_QUERY );
122cdf0e10cSrcweir     return xSheet->getCellByPosition( rPos.Column, rPos.Row );
123cdf0e10cSrcweir }
124cdf0e10cSrcweir 
lcl_SetValue(const uno::Reference<sheet::XSpreadsheetDocument> & xDoc,const table::CellAddress & rPos,double fValue)125cdf0e10cSrcweir void lcl_SetValue( const uno::Reference<sheet::XSpreadsheetDocument>& xDoc,
126cdf0e10cSrcweir                    const table::CellAddress& rPos, double fValue )
127cdf0e10cSrcweir {
128cdf0e10cSrcweir     lcl_GetCell( xDoc, rPos )->setValue( fValue );
129cdf0e10cSrcweir }
130cdf0e10cSrcweir 
lcl_GetValue(const uno::Reference<sheet::XSpreadsheetDocument> & xDoc,const table::CellAddress & rPos)131cdf0e10cSrcweir double lcl_GetValue( const uno::Reference<sheet::XSpreadsheetDocument>& xDoc,
132cdf0e10cSrcweir                      const table::CellAddress& rPos )
133cdf0e10cSrcweir {
134cdf0e10cSrcweir     return lcl_GetCell( xDoc, rPos )->getValue();
135cdf0e10cSrcweir }
136cdf0e10cSrcweir 
137cdf0e10cSrcweir // -------------------------------------------------------------------------
138cdf0e10cSrcweir 
SolverComponent(const uno::Reference<uno::XComponentContext> &)139cdf0e10cSrcweir SolverComponent::SolverComponent( const uno::Reference<uno::XComponentContext>& /* rSMgr */ ) :
140cdf0e10cSrcweir     OPropertyContainer( GetBroadcastHelper() ),
141cdf0e10cSrcweir     mbMaximize( sal_True ),
142cdf0e10cSrcweir     mbNonNegative( sal_False ),
143cdf0e10cSrcweir     mbInteger( sal_False ),
144e74d7e62SPedro Giffuni     mnTimeout( 120 ),
145cdf0e10cSrcweir     mnEpsilonLevel( 0 ),
146cdf0e10cSrcweir     mbLimitBBDepth( sal_True ),
147b1f1da72SPedro Giffuni     mbNonLinearTest( sal_True ),
148cdf0e10cSrcweir     mbSuccess( sal_False ),
149cdf0e10cSrcweir     mfResultValue( 0.0 )
150cdf0e10cSrcweir {
151cdf0e10cSrcweir     // for XPropertySet implementation:
152cdf0e10cSrcweir     registerProperty( C2U(STR_NONNEGATIVE),  PROP_NONNEGATIVE,  0, &mbNonNegative,  getCppuType( &mbNonNegative )  );
153cdf0e10cSrcweir     registerProperty( C2U(STR_INTEGER),      PROP_INTEGER,      0, &mbInteger,      getCppuType( &mbInteger )      );
154cdf0e10cSrcweir     registerProperty( C2U(STR_TIMEOUT),      PROP_TIMEOUT,      0, &mnTimeout,      getCppuType( &mnTimeout )      );
155cdf0e10cSrcweir     registerProperty( C2U(STR_EPSILONLEVEL), PROP_EPSILONLEVEL, 0, &mnEpsilonLevel, getCppuType( &mnEpsilonLevel ) );
156cdf0e10cSrcweir     registerProperty( C2U(STR_LIMITBBDEPTH), PROP_LIMITBBDEPTH, 0, &mbLimitBBDepth, getCppuType( &mbLimitBBDepth ) );
157b1f1da72SPedro Giffuni     registerProperty( C2U(STR_NONLINEARTEST), PROP_NONLINEARTEST, 0, &mbNonLinearTest, getCppuType( &mbNonLinearTest ) );
158cdf0e10cSrcweir }
159cdf0e10cSrcweir 
~SolverComponent()160cdf0e10cSrcweir SolverComponent::~SolverComponent()
161cdf0e10cSrcweir {
162cdf0e10cSrcweir }
163cdf0e10cSrcweir 
IMPLEMENT_FORWARD_XINTERFACE2(SolverComponent,SolverComponent_Base,OPropertyContainer)164cdf0e10cSrcweir IMPLEMENT_FORWARD_XINTERFACE2( SolverComponent, SolverComponent_Base, OPropertyContainer )
165cdf0e10cSrcweir IMPLEMENT_FORWARD_XTYPEPROVIDER2( SolverComponent, SolverComponent_Base, OPropertyContainer )
166cdf0e10cSrcweir 
167cdf0e10cSrcweir cppu::IPropertyArrayHelper* SolverComponent::createArrayHelper() const
168cdf0e10cSrcweir {
169cdf0e10cSrcweir     uno::Sequence<beans::Property> aProps;
170cdf0e10cSrcweir     describeProperties( aProps );
171cdf0e10cSrcweir     return new cppu::OPropertyArrayHelper( aProps );
172cdf0e10cSrcweir }
173cdf0e10cSrcweir 
getInfoHelper()174cdf0e10cSrcweir cppu::IPropertyArrayHelper& SAL_CALL SolverComponent::getInfoHelper()
175cdf0e10cSrcweir {
176cdf0e10cSrcweir     return *getArrayHelper();
177cdf0e10cSrcweir }
178cdf0e10cSrcweir 
getPropertySetInfo()179cdf0e10cSrcweir uno::Reference<beans::XPropertySetInfo> SAL_CALL SolverComponent::getPropertySetInfo() throw(uno::RuntimeException)
180cdf0e10cSrcweir {
181cdf0e10cSrcweir     return createPropertySetInfo( getInfoHelper() );
182cdf0e10cSrcweir }
183cdf0e10cSrcweir 
184cdf0e10cSrcweir // XSolverDescription
185cdf0e10cSrcweir 
getComponentDescription()186cdf0e10cSrcweir OUString SAL_CALL SolverComponent::getComponentDescription() throw (uno::RuntimeException)
187cdf0e10cSrcweir {
188cdf0e10cSrcweir     return lcl_GetResourceString( RID_SOLVER_COMPONENT );
189cdf0e10cSrcweir }
190cdf0e10cSrcweir 
getStatusDescription()191cdf0e10cSrcweir OUString SAL_CALL SolverComponent::getStatusDescription() throw (uno::RuntimeException)
192cdf0e10cSrcweir {
193cdf0e10cSrcweir     return maStatus;
194cdf0e10cSrcweir }
195cdf0e10cSrcweir 
getPropertyDescription(const OUString & rPropertyName)196cdf0e10cSrcweir OUString SAL_CALL SolverComponent::getPropertyDescription( const OUString& rPropertyName ) throw (uno::RuntimeException)
197cdf0e10cSrcweir {
198cdf0e10cSrcweir     sal_uInt32 nResId = 0;
199cdf0e10cSrcweir     sal_Int32 nHandle = getInfoHelper().getHandleByName( rPropertyName );
200cdf0e10cSrcweir     switch (nHandle)
201cdf0e10cSrcweir     {
202cdf0e10cSrcweir         case PROP_NONNEGATIVE:
203cdf0e10cSrcweir             nResId = RID_PROPERTY_NONNEGATIVE;
204cdf0e10cSrcweir             break;
205cdf0e10cSrcweir         case PROP_INTEGER:
206cdf0e10cSrcweir             nResId = RID_PROPERTY_INTEGER;
207cdf0e10cSrcweir             break;
208cdf0e10cSrcweir         case PROP_TIMEOUT:
209cdf0e10cSrcweir             nResId = RID_PROPERTY_TIMEOUT;
210cdf0e10cSrcweir             break;
211cdf0e10cSrcweir         case PROP_EPSILONLEVEL:
212cdf0e10cSrcweir             nResId = RID_PROPERTY_EPSILONLEVEL;
213cdf0e10cSrcweir             break;
214cdf0e10cSrcweir         case PROP_LIMITBBDEPTH:
215cdf0e10cSrcweir             nResId = RID_PROPERTY_LIMITBBDEPTH;
216cdf0e10cSrcweir             break;
217b1f1da72SPedro Giffuni         case PROP_NONLINEARTEST:
218b1f1da72SPedro Giffuni             nResId = RID_PROPERTY_NONLINEARTEST;
219b1f1da72SPedro Giffuni             break;
220cdf0e10cSrcweir         default:
221cdf0e10cSrcweir             {
222cdf0e10cSrcweir                 // unknown - leave empty
223cdf0e10cSrcweir             }
224cdf0e10cSrcweir     }
225cdf0e10cSrcweir     OUString aRet;
226cdf0e10cSrcweir     if ( nResId )
227cdf0e10cSrcweir         aRet = lcl_GetResourceString( nResId );
228cdf0e10cSrcweir     return aRet;
229cdf0e10cSrcweir }
230cdf0e10cSrcweir 
231cdf0e10cSrcweir // XSolver: settings
232cdf0e10cSrcweir 
getDocument()233cdf0e10cSrcweir uno::Reference<sheet::XSpreadsheetDocument> SAL_CALL SolverComponent::getDocument() throw(uno::RuntimeException)
234cdf0e10cSrcweir {
235cdf0e10cSrcweir     return mxDoc;
236cdf0e10cSrcweir }
237cdf0e10cSrcweir 
setDocument(const uno::Reference<sheet::XSpreadsheetDocument> & _document)238cdf0e10cSrcweir void SAL_CALL SolverComponent::setDocument( const uno::Reference<sheet::XSpreadsheetDocument>& _document )
239cdf0e10cSrcweir                                 throw(uno::RuntimeException)
240cdf0e10cSrcweir {
241cdf0e10cSrcweir     mxDoc = _document;
242cdf0e10cSrcweir }
243cdf0e10cSrcweir 
getObjective()244cdf0e10cSrcweir table::CellAddress SAL_CALL SolverComponent::getObjective() throw(uno::RuntimeException)
245cdf0e10cSrcweir {
246cdf0e10cSrcweir     return maObjective;
247cdf0e10cSrcweir }
248cdf0e10cSrcweir 
setObjective(const table::CellAddress & _objective)249cdf0e10cSrcweir void SAL_CALL SolverComponent::setObjective( const table::CellAddress& _objective ) throw(uno::RuntimeException)
250cdf0e10cSrcweir {
251cdf0e10cSrcweir     maObjective = _objective;
252cdf0e10cSrcweir }
253cdf0e10cSrcweir 
getVariables()254cdf0e10cSrcweir uno::Sequence<table::CellAddress> SAL_CALL SolverComponent::getVariables() throw(uno::RuntimeException)
255cdf0e10cSrcweir {
256cdf0e10cSrcweir     return maVariables;
257cdf0e10cSrcweir }
258cdf0e10cSrcweir 
setVariables(const uno::Sequence<table::CellAddress> & _variables)259cdf0e10cSrcweir void SAL_CALL SolverComponent::setVariables( const uno::Sequence<table::CellAddress>& _variables )
260cdf0e10cSrcweir                                 throw(uno::RuntimeException)
261cdf0e10cSrcweir {
262cdf0e10cSrcweir     maVariables = _variables;
263cdf0e10cSrcweir }
264cdf0e10cSrcweir 
getConstraints()265cdf0e10cSrcweir uno::Sequence<sheet::SolverConstraint> SAL_CALL SolverComponent::getConstraints() throw(uno::RuntimeException)
266cdf0e10cSrcweir {
267cdf0e10cSrcweir     return maConstraints;
268cdf0e10cSrcweir }
269cdf0e10cSrcweir 
setConstraints(const uno::Sequence<sheet::SolverConstraint> & _constraints)270cdf0e10cSrcweir void SAL_CALL SolverComponent::setConstraints( const uno::Sequence<sheet::SolverConstraint>& _constraints )
271cdf0e10cSrcweir                                 throw(uno::RuntimeException)
272cdf0e10cSrcweir {
273cdf0e10cSrcweir     maConstraints = _constraints;
274cdf0e10cSrcweir }
275cdf0e10cSrcweir 
getMaximize()276cdf0e10cSrcweir sal_Bool SAL_CALL SolverComponent::getMaximize() throw(uno::RuntimeException)
277cdf0e10cSrcweir {
278cdf0e10cSrcweir     return mbMaximize;
279cdf0e10cSrcweir }
280cdf0e10cSrcweir 
setMaximize(sal_Bool _maximize)281cdf0e10cSrcweir void SAL_CALL SolverComponent::setMaximize( sal_Bool _maximize ) throw(uno::RuntimeException)
282cdf0e10cSrcweir {
283cdf0e10cSrcweir     mbMaximize = _maximize;
284cdf0e10cSrcweir }
285cdf0e10cSrcweir 
286cdf0e10cSrcweir // XSolver: get results
287cdf0e10cSrcweir 
getSuccess()288cdf0e10cSrcweir sal_Bool SAL_CALL SolverComponent::getSuccess() throw(uno::RuntimeException)
289cdf0e10cSrcweir {
290cdf0e10cSrcweir     return mbSuccess;
291cdf0e10cSrcweir }
292cdf0e10cSrcweir 
getResultValue()293cdf0e10cSrcweir double SAL_CALL SolverComponent::getResultValue() throw(uno::RuntimeException)
294cdf0e10cSrcweir {
295cdf0e10cSrcweir     return mfResultValue;
296cdf0e10cSrcweir }
297cdf0e10cSrcweir 
getSolution()298cdf0e10cSrcweir uno::Sequence<double> SAL_CALL SolverComponent::getSolution() throw(uno::RuntimeException)
299cdf0e10cSrcweir {
300cdf0e10cSrcweir     return maSolution;
301cdf0e10cSrcweir }
302cdf0e10cSrcweir 
303cdf0e10cSrcweir // -------------------------------------------------------------------------
304cdf0e10cSrcweir 
solve()305cdf0e10cSrcweir void SAL_CALL SolverComponent::solve() throw(uno::RuntimeException)
306cdf0e10cSrcweir {
307cdf0e10cSrcweir     uno::Reference<frame::XModel> xModel( mxDoc, uno::UNO_QUERY );
308cdf0e10cSrcweir     if ( !xModel.is() )
309cdf0e10cSrcweir         throw uno::RuntimeException();
310cdf0e10cSrcweir 
311cdf0e10cSrcweir     maStatus = OUString();
312cdf0e10cSrcweir     mbSuccess = false;
313cdf0e10cSrcweir 
314cdf0e10cSrcweir     xModel->lockControllers();
315cdf0e10cSrcweir 
316cdf0e10cSrcweir     // collect variables in vector (?)
317cdf0e10cSrcweir 
318cdf0e10cSrcweir     std::vector<table::CellAddress> aVariableCells;
319cdf0e10cSrcweir     for (sal_Int32 nPos=0; nPos<maVariables.getLength(); nPos++)
320cdf0e10cSrcweir         aVariableCells.push_back( maVariables[nPos] );
321cdf0e10cSrcweir     size_t nVariables = aVariableCells.size();
322cdf0e10cSrcweir     size_t nVar = 0;
323cdf0e10cSrcweir 
324cdf0e10cSrcweir     // collect all dependent cells
325cdf0e10cSrcweir 
326cdf0e10cSrcweir     ScSolverCellHashMap aCellsHash;
327cdf0e10cSrcweir     aCellsHash[maObjective].reserve( nVariables + 1 );                  // objective function
328cdf0e10cSrcweir 
329cdf0e10cSrcweir     for (sal_Int32 nConstrPos = 0; nConstrPos < maConstraints.getLength(); ++nConstrPos)
330cdf0e10cSrcweir     {
331cdf0e10cSrcweir         table::CellAddress aCellAddr = maConstraints[nConstrPos].Left;
332cdf0e10cSrcweir         aCellsHash[aCellAddr].reserve( nVariables + 1 );                // constraints: left hand side
333cdf0e10cSrcweir 
334cdf0e10cSrcweir         if ( maConstraints[nConstrPos].Right >>= aCellAddr )
335cdf0e10cSrcweir             aCellsHash[aCellAddr].reserve( nVariables + 1 );            // constraints: right hand side
336cdf0e10cSrcweir     }
337cdf0e10cSrcweir 
338cdf0e10cSrcweir     // set all variables to zero
339cdf0e10cSrcweir     //! store old values?
340cdf0e10cSrcweir     //! use old values as initial values?
341cdf0e10cSrcweir     std::vector<table::CellAddress>::const_iterator aVarIter;
342cdf0e10cSrcweir     for ( aVarIter = aVariableCells.begin(); aVarIter != aVariableCells.end(); ++aVarIter )
343cdf0e10cSrcweir     {
344cdf0e10cSrcweir         lcl_SetValue( mxDoc, *aVarIter, 0.0 );
345cdf0e10cSrcweir     }
346cdf0e10cSrcweir 
347cdf0e10cSrcweir     // read initial values from all dependent cells
348cdf0e10cSrcweir     ScSolverCellHashMap::iterator aCellsIter;
349cdf0e10cSrcweir     for ( aCellsIter = aCellsHash.begin(); aCellsIter != aCellsHash.end(); ++aCellsIter )
350cdf0e10cSrcweir     {
351cdf0e10cSrcweir         double fValue = lcl_GetValue( mxDoc, aCellsIter->first );
352cdf0e10cSrcweir         aCellsIter->second.push_back( fValue );                         // store as first element, as-is
353cdf0e10cSrcweir     }
354cdf0e10cSrcweir 
355cdf0e10cSrcweir     // loop through variables
356cdf0e10cSrcweir     for ( aVarIter = aVariableCells.begin(); aVarIter != aVariableCells.end(); ++aVarIter )
357cdf0e10cSrcweir     {
358cdf0e10cSrcweir         lcl_SetValue( mxDoc, *aVarIter, 1.0 );      // set to 1 to examine influence
359cdf0e10cSrcweir 
360cdf0e10cSrcweir         // read value change from all dependent cells
361cdf0e10cSrcweir         for ( aCellsIter = aCellsHash.begin(); aCellsIter != aCellsHash.end(); ++aCellsIter )
362cdf0e10cSrcweir         {
363cdf0e10cSrcweir             double fChanged = lcl_GetValue( mxDoc, aCellsIter->first );
364cdf0e10cSrcweir             double fInitial = aCellsIter->second.front();
365cdf0e10cSrcweir             aCellsIter->second.push_back( fChanged - fInitial );
366cdf0e10cSrcweir         }
367cdf0e10cSrcweir 
368cdf0e10cSrcweir         lcl_SetValue( mxDoc, *aVarIter, 2.0 );      // minimal test for linearity
369cdf0e10cSrcweir 
370cdf0e10cSrcweir         for ( aCellsIter = aCellsHash.begin(); aCellsIter != aCellsHash.end(); ++aCellsIter )
371cdf0e10cSrcweir         {
372cdf0e10cSrcweir             double fInitial = aCellsIter->second.front();
373cdf0e10cSrcweir             double fCoeff   = aCellsIter->second.back();       // last appended: coefficient for this variable
374cdf0e10cSrcweir             double fTwo     = lcl_GetValue( mxDoc, aCellsIter->first );
375b1f1da72SPedro Giffuni 
376b1f1da72SPedro Giffuni           if ( mbNonLinearTest )
377b1f1da72SPedro Giffuni           {
378b1f1da72SPedro Giffuni               bool bLinear ( sal_True );
379b1f1da72SPedro Giffuni               bLinear = rtl::math::approxEqual( fTwo, fInitial + 2.0 * fCoeff ) ||
380b1f1da72SPedro Giffuni               rtl::math::approxEqual( fInitial, fTwo - 2.0 * fCoeff );
381b1f1da72SPedro Giffuni             // second comparison is needed in case fTwo is zero
382b1f1da72SPedro Giffuni               if ( !bLinear )
383b1f1da72SPedro Giffuni                   maStatus = lcl_GetResourceString( RID_ERROR_NONLINEAR );
384b1f1da72SPedro Giffuni            }
385cdf0e10cSrcweir         }
386cdf0e10cSrcweir 
387cdf0e10cSrcweir         lcl_SetValue( mxDoc, *aVarIter, 0.0 );      // set back to zero for examining next variable
388cdf0e10cSrcweir     }
389cdf0e10cSrcweir 
390cdf0e10cSrcweir     xModel->unlockControllers();
391cdf0e10cSrcweir 
392cdf0e10cSrcweir     if ( maStatus.getLength() )
393cdf0e10cSrcweir         return;
394cdf0e10cSrcweir 
395cdf0e10cSrcweir     //
39670840c81SAndre Fischer     // build parameter arrays for CoinMP
397cdf0e10cSrcweir     //
398cdf0e10cSrcweir 
399cdf0e10cSrcweir     // set objective function
400cdf0e10cSrcweir 
401cdf0e10cSrcweir     const std::vector<double>& rObjCoeff = aCellsHash[maObjective];
40270840c81SAndre Fischer     double* pObjectCoeffs = new double[nVariables];
403cdf0e10cSrcweir     for (nVar=0; nVar<nVariables; nVar++)
40470840c81SAndre Fischer         pObjectCoeffs[nVar] = rObjCoeff[nVar+1];
40570840c81SAndre Fischer     double nObjectConst = rObjCoeff[0];             // constant term of objective
406cdf0e10cSrcweir 
407cdf0e10cSrcweir     // add rows
408cdf0e10cSrcweir 
40970840c81SAndre Fischer     size_t nRows = maConstraints.getLength();
41070840c81SAndre Fischer     size_t nCompSize = nVariables * nRows;
41170840c81SAndre Fischer     double* pCompMatrix = new double[nCompSize];    // first collect all coefficients, row-wise
41270840c81SAndre Fischer     for (size_t i=0; i<nCompSize; i++)
41370840c81SAndre Fischer         pCompMatrix[i] = 0.0;
41470840c81SAndre Fischer 
41570840c81SAndre Fischer     double* pRHS = new double[nRows];
41670840c81SAndre Fischer     char* pRowType = new char[nRows];
41770840c81SAndre Fischer     for (size_t i=0; i<nRows; i++)
41870840c81SAndre Fischer     {
41970840c81SAndre Fischer         pRHS[i] = 0.0;
42070840c81SAndre Fischer         pRowType[i] = 'N';
42170840c81SAndre Fischer     }
422cdf0e10cSrcweir 
423cdf0e10cSrcweir     for (sal_Int32 nConstrPos = 0; nConstrPos < maConstraints.getLength(); ++nConstrPos)
424cdf0e10cSrcweir     {
425cdf0e10cSrcweir         // integer constraints are set later
426cdf0e10cSrcweir         sheet::SolverConstraintOperator eOp = maConstraints[nConstrPos].Operator;
427cdf0e10cSrcweir         if ( eOp == sheet::SolverConstraintOperator_LESS_EQUAL ||
428cdf0e10cSrcweir              eOp == sheet::SolverConstraintOperator_GREATER_EQUAL ||
429cdf0e10cSrcweir              eOp == sheet::SolverConstraintOperator_EQUAL )
430cdf0e10cSrcweir         {
431cdf0e10cSrcweir             double fDirectValue = 0.0;
432cdf0e10cSrcweir             bool bRightCell = false;
433cdf0e10cSrcweir             table::CellAddress aRightAddr;
434cdf0e10cSrcweir             const uno::Any& rRightAny = maConstraints[nConstrPos].Right;
435cdf0e10cSrcweir             if ( rRightAny >>= aRightAddr )
436cdf0e10cSrcweir                 bRightCell = true;                  // cell specified as right-hand side
437cdf0e10cSrcweir             else
438cdf0e10cSrcweir                 rRightAny >>= fDirectValue;         // constant value
439cdf0e10cSrcweir 
440cdf0e10cSrcweir             table::CellAddress aLeftAddr = maConstraints[nConstrPos].Left;
441cdf0e10cSrcweir 
442cdf0e10cSrcweir             const std::vector<double>& rLeftCoeff = aCellsHash[aLeftAddr];
44370840c81SAndre Fischer             double* pValues = &pCompMatrix[nConstrPos * nVariables];
444cdf0e10cSrcweir             for (nVar=0; nVar<nVariables; nVar++)
44570840c81SAndre Fischer                 pValues[nVar] = rLeftCoeff[nVar+1];
446cdf0e10cSrcweir 
447cdf0e10cSrcweir             // if left hand cell has a constant term, put into rhs value
448cdf0e10cSrcweir             double fRightValue = -rLeftCoeff[0];
449cdf0e10cSrcweir 
450cdf0e10cSrcweir             if ( bRightCell )
451cdf0e10cSrcweir             {
452cdf0e10cSrcweir                 const std::vector<double>& rRightCoeff = aCellsHash[aRightAddr];
453cdf0e10cSrcweir                 // modify pValues with rhs coefficients
454cdf0e10cSrcweir                 for (nVar=0; nVar<nVariables; nVar++)
45570840c81SAndre Fischer                     pValues[nVar] -= rRightCoeff[nVar+1];
456cdf0e10cSrcweir 
457cdf0e10cSrcweir                 fRightValue += rRightCoeff[0];      // constant term
458cdf0e10cSrcweir             }
459cdf0e10cSrcweir             else
460cdf0e10cSrcweir                 fRightValue += fDirectValue;
461cdf0e10cSrcweir 
462cdf0e10cSrcweir             switch ( eOp )
463cdf0e10cSrcweir             {
46470840c81SAndre Fischer                 case sheet::SolverConstraintOperator_LESS_EQUAL:    pRowType[nConstrPos] = 'L'; break;
46570840c81SAndre Fischer                 case sheet::SolverConstraintOperator_GREATER_EQUAL: pRowType[nConstrPos] = 'G'; break;
46670840c81SAndre Fischer                 case sheet::SolverConstraintOperator_EQUAL:         pRowType[nConstrPos] = 'E'; break;
467cdf0e10cSrcweir                 default:
468cdf0e10cSrcweir                     OSL_ENSURE( false, "unexpected enum type" );
469cdf0e10cSrcweir             }
47070840c81SAndre Fischer             pRHS[nConstrPos] = fRightValue;
471cdf0e10cSrcweir         }
472cdf0e10cSrcweir     }
473cdf0e10cSrcweir 
47470840c81SAndre Fischer     // Find non-zero coefficients, column-wise
47570840c81SAndre Fischer 
47670840c81SAndre Fischer     int* pMatrixBegin = new int[nVariables+1];
47770840c81SAndre Fischer     int* pMatrixCount = new int[nVariables];
47870840c81SAndre Fischer     double* pMatrix = new double[nCompSize];    // not always completely used
47970840c81SAndre Fischer     int* pMatrixIndex = new int[nCompSize];
48070840c81SAndre Fischer     int nMatrixPos = 0;
48170840c81SAndre Fischer     for (nVar=0; nVar<nVariables; nVar++)
48270840c81SAndre Fischer     {
48370840c81SAndre Fischer         int nBegin = nMatrixPos;
48470840c81SAndre Fischer         for (size_t nRow=0; nRow<nRows; nRow++)
48570840c81SAndre Fischer         {
48670840c81SAndre Fischer             double fCoeff = pCompMatrix[ nRow * nVariables + nVar ];    // row-wise
48770840c81SAndre Fischer             if ( fCoeff != 0.0 )
48870840c81SAndre Fischer             {
48970840c81SAndre Fischer                 pMatrix[nMatrixPos] = fCoeff;
49070840c81SAndre Fischer                 pMatrixIndex[nMatrixPos] = nRow;
49170840c81SAndre Fischer                 ++nMatrixPos;
49270840c81SAndre Fischer             }
49370840c81SAndre Fischer         }
49470840c81SAndre Fischer         pMatrixBegin[nVar] = nBegin;
49570840c81SAndre Fischer         pMatrixCount[nVar] = nMatrixPos - nBegin;
49670840c81SAndre Fischer     }
49770840c81SAndre Fischer     pMatrixBegin[nVariables] = nMatrixPos;
49870840c81SAndre Fischer     delete[] pCompMatrix;
49970840c81SAndre Fischer     pCompMatrix = NULL;
500cdf0e10cSrcweir 
501cdf0e10cSrcweir     // apply settings to all variables
502cdf0e10cSrcweir 
50370840c81SAndre Fischer     double* pLowerBounds = new double[nVariables];
50470840c81SAndre Fischer     double* pUpperBounds = new double[nVariables];
505cdf0e10cSrcweir     for (nVar=0; nVar<nVariables; nVar++)
506cdf0e10cSrcweir     {
50770840c81SAndre Fischer         pLowerBounds[nVar] = mbNonNegative ? 0.0 : -DBL_MAX;
50870840c81SAndre Fischer         pUpperBounds[nVar] = DBL_MAX;
50970840c81SAndre Fischer 
51070840c81SAndre Fischer         // bounds could possibly be further restricted from single-cell constraints
511cdf0e10cSrcweir     }
512cdf0e10cSrcweir 
51370840c81SAndre Fischer     char* pColType = new char[nVariables];
51470840c81SAndre Fischer     for (nVar=0; nVar<nVariables; nVar++)
51570840c81SAndre Fischer         pColType[nVar] = mbInteger ? 'I' : 'C';
51670840c81SAndre Fischer 
517cdf0e10cSrcweir     // apply single-var integer constraints
518cdf0e10cSrcweir 
519cdf0e10cSrcweir     for (sal_Int32 nConstrPos = 0; nConstrPos < maConstraints.getLength(); ++nConstrPos)
520cdf0e10cSrcweir     {
521cdf0e10cSrcweir         sheet::SolverConstraintOperator eOp = maConstraints[nConstrPos].Operator;
522cdf0e10cSrcweir         if ( eOp == sheet::SolverConstraintOperator_INTEGER ||
523cdf0e10cSrcweir              eOp == sheet::SolverConstraintOperator_BINARY )
524cdf0e10cSrcweir         {
525cdf0e10cSrcweir             table::CellAddress aLeftAddr = maConstraints[nConstrPos].Left;
526cdf0e10cSrcweir             // find variable index for cell
527cdf0e10cSrcweir             for (nVar=0; nVar<nVariables; nVar++)
528cdf0e10cSrcweir                 if ( AddressEqual( aVariableCells[nVar], aLeftAddr ) )
529cdf0e10cSrcweir                 {
530cdf0e10cSrcweir                     if ( eOp == sheet::SolverConstraintOperator_INTEGER )
53170840c81SAndre Fischer                         pColType[nVar] = 'I';
532cdf0e10cSrcweir                     else
53370840c81SAndre Fischer                     {
53470840c81SAndre Fischer                         pColType[nVar] = 'B';
53570840c81SAndre Fischer                         pLowerBounds[nVar] = 0.0;
53670840c81SAndre Fischer                         pUpperBounds[nVar] = 1.0;
53770840c81SAndre Fischer                     }
538cdf0e10cSrcweir                 }
539cdf0e10cSrcweir         }
540cdf0e10cSrcweir     }
541cdf0e10cSrcweir 
54270840c81SAndre Fischer     int nObjectSense = mbMaximize ? SOLV_OBJSENS_MAX : SOLV_OBJSENS_MIN;
543cdf0e10cSrcweir 
54470840c81SAndre Fischer     HPROB hProb = CoinCreateProblem("");
54570840c81SAndre Fischer     int nResult = CoinLoadProblem( hProb, nVariables, nRows, nMatrixPos, 0,
54670840c81SAndre Fischer                     nObjectSense, nObjectConst, pObjectCoeffs,
54770840c81SAndre Fischer                     pLowerBounds, pUpperBounds, pRowType, pRHS, NULL,
54870840c81SAndre Fischer                     pMatrixBegin, pMatrixCount, pMatrixIndex, pMatrix,
54970840c81SAndre Fischer                     NULL, NULL, NULL );
55070840c81SAndre Fischer     nResult = CoinLoadInteger( hProb, pColType );
551cdf0e10cSrcweir 
55270840c81SAndre Fischer     delete[] pColType;
55370840c81SAndre Fischer     delete[] pMatrixIndex;
55470840c81SAndre Fischer     delete[] pMatrix;
55570840c81SAndre Fischer     delete[] pMatrixCount;
55670840c81SAndre Fischer     delete[] pMatrixBegin;
55770840c81SAndre Fischer     delete[] pUpperBounds;
55870840c81SAndre Fischer     delete[] pLowerBounds;
55970840c81SAndre Fischer     delete[] pRowType;
56070840c81SAndre Fischer     delete[] pRHS;
56170840c81SAndre Fischer     delete[] pObjectCoeffs;
56270840c81SAndre Fischer 
56370840c81SAndre Fischer     CoinSetRealOption( hProb, COIN_REAL_MAXSECONDS, mnTimeout );
56470840c81SAndre Fischer     CoinSetRealOption( hProb, COIN_REAL_MIPMAXSEC, mnTimeout );
56570840c81SAndre Fischer 
56670840c81SAndre Fischer     // TODO: handle (or remove) settings: epsilon, B&B depth
567cdf0e10cSrcweir 
568cdf0e10cSrcweir     // solve model
569cdf0e10cSrcweir 
57070840c81SAndre Fischer     nResult = CoinCheckProblem( hProb );
57170840c81SAndre Fischer     nResult = CoinOptimizeProblem( hProb, 0 );
572cdf0e10cSrcweir 
57370840c81SAndre Fischer     mbSuccess = ( nResult == SOLV_CALL_SUCCESS );
574cdf0e10cSrcweir     if ( mbSuccess )
575cdf0e10cSrcweir     {
576cdf0e10cSrcweir         // get solution
577cdf0e10cSrcweir 
578cdf0e10cSrcweir         maSolution.realloc( nVariables );
57970840c81SAndre Fischer         CoinGetSolutionValues( hProb, maSolution.getArray(), NULL, NULL, NULL );
58070840c81SAndre Fischer         mfResultValue = CoinGetObjectValue( hProb );
581cdf0e10cSrcweir     }
58270840c81SAndre Fischer     else
58370840c81SAndre Fischer     {
58470840c81SAndre Fischer         int nSolutionStatus = CoinGetSolutionStatus( hProb );
58570840c81SAndre Fischer         if ( nSolutionStatus == 1 )
586cdf0e10cSrcweir             maStatus = lcl_GetResourceString( RID_ERROR_INFEASIBLE );
58770840c81SAndre Fischer         else if ( nSolutionStatus == 2 )
588cdf0e10cSrcweir             maStatus = lcl_GetResourceString( RID_ERROR_UNBOUNDED );
58970840c81SAndre Fischer         // TODO: detect timeout condition and report as RID_ERROR_TIMEOUT
59070840c81SAndre Fischer         // (currently reported as infeasible)
59170840c81SAndre Fischer     }
592cdf0e10cSrcweir 
59370840c81SAndre Fischer     CoinUnloadProblem( hProb );
594cdf0e10cSrcweir }
595cdf0e10cSrcweir 
596cdf0e10cSrcweir // -------------------------------------------------------------------------
597cdf0e10cSrcweir 
598cdf0e10cSrcweir // XServiceInfo
599cdf0e10cSrcweir 
SolverComponent_getSupportedServiceNames()600cdf0e10cSrcweir uno::Sequence< OUString > SolverComponent_getSupportedServiceNames()
601cdf0e10cSrcweir {
602cdf0e10cSrcweir     uno::Sequence< OUString > aServiceNames( 1 );
603cdf0e10cSrcweir     aServiceNames[ 0 ] = OUString::createFromAscii( "com.sun.star.sheet.Solver" );
604cdf0e10cSrcweir     return aServiceNames;
605cdf0e10cSrcweir }
606cdf0e10cSrcweir 
SolverComponent_getImplementationName()607cdf0e10cSrcweir OUString SolverComponent_getImplementationName()
608cdf0e10cSrcweir {
609cdf0e10cSrcweir     return OUString::createFromAscii( "com.sun.star.comp.Calc.Solver" );
610cdf0e10cSrcweir }
611cdf0e10cSrcweir 
getImplementationName()612cdf0e10cSrcweir OUString SAL_CALL SolverComponent::getImplementationName() throw(uno::RuntimeException)
613cdf0e10cSrcweir {
614cdf0e10cSrcweir     return SolverComponent_getImplementationName();
615cdf0e10cSrcweir }
616cdf0e10cSrcweir 
supportsService(const OUString & rServiceName)617cdf0e10cSrcweir sal_Bool SAL_CALL SolverComponent::supportsService( const OUString& rServiceName ) throw(uno::RuntimeException)
618cdf0e10cSrcweir {
619cdf0e10cSrcweir     const uno::Sequence< OUString > aServices = SolverComponent_getSupportedServiceNames();
620cdf0e10cSrcweir     const OUString* pArray = aServices.getConstArray();
621cdf0e10cSrcweir     const OUString* pArrayEnd = pArray + aServices.getLength();
622cdf0e10cSrcweir     return ::std::find( pArray, pArrayEnd, rServiceName ) != pArrayEnd;
623cdf0e10cSrcweir }
624cdf0e10cSrcweir 
getSupportedServiceNames()625cdf0e10cSrcweir uno::Sequence<OUString> SAL_CALL SolverComponent::getSupportedServiceNames() throw(uno::RuntimeException)
626cdf0e10cSrcweir {
627cdf0e10cSrcweir     return SolverComponent_getSupportedServiceNames();
628cdf0e10cSrcweir }
629cdf0e10cSrcweir 
SolverComponent_createInstance(const uno::Reference<uno::XComponentContext> & rSMgr)630cdf0e10cSrcweir uno::Reference<uno::XInterface> SolverComponent_createInstance( const uno::Reference<uno::XComponentContext>& rSMgr )
631cdf0e10cSrcweir     throw(uno::Exception)
632cdf0e10cSrcweir {
633cdf0e10cSrcweir     return (cppu::OWeakObject*) new SolverComponent( rSMgr );
634cdf0e10cSrcweir }
635cdf0e10cSrcweir 
636cdf0e10cSrcweir // -------------------------------------------------------------------------
637cdf0e10cSrcweir 
638cdf0e10cSrcweir extern "C"
639cdf0e10cSrcweir {
component_getImplementationEnvironment(const sal_Char ** ppEnvTypeName,uno_Environment **)640cdf0e10cSrcweir     SAL_DLLPUBLIC_EXPORT void SAL_CALL component_getImplementationEnvironment(
641cdf0e10cSrcweir         const sal_Char ** ppEnvTypeName, uno_Environment ** )
642cdf0e10cSrcweir     {
643cdf0e10cSrcweir         *ppEnvTypeName = CPPU_CURRENT_LANGUAGE_BINDING_NAME;
644cdf0e10cSrcweir     }
645cdf0e10cSrcweir 
646cdf0e10cSrcweir     // -------------------------------------------------------------------------
647cdf0e10cSrcweir 
component_getFactory(const sal_Char * pImplName,void * pServiceManager,void *)648cdf0e10cSrcweir     SAL_DLLPUBLIC_EXPORT void* SAL_CALL component_getFactory( const sal_Char * pImplName, void * pServiceManager, void * /*pRegistryKey*/ )
649cdf0e10cSrcweir     {
650cdf0e10cSrcweir         OUString    aImplName( OUString::createFromAscii( pImplName ) );
651cdf0e10cSrcweir         void*       pRet = 0;
652cdf0e10cSrcweir 
653cdf0e10cSrcweir         if( pServiceManager )
654cdf0e10cSrcweir         {
655cdf0e10cSrcweir             uno::Reference< lang::XSingleComponentFactory > xFactory;
656cdf0e10cSrcweir             if( aImplName.equals( SolverComponent_getImplementationName() ) )
657cdf0e10cSrcweir                 xFactory = cppu::createSingleComponentFactory(
658cdf0e10cSrcweir                         SolverComponent_createInstance,
659cdf0e10cSrcweir                         OUString::createFromAscii( pImplName ),
660cdf0e10cSrcweir                         SolverComponent_getSupportedServiceNames() );
661cdf0e10cSrcweir 
662cdf0e10cSrcweir             if( xFactory.is() )
663cdf0e10cSrcweir             {
664cdf0e10cSrcweir                 xFactory->acquire();
665cdf0e10cSrcweir                 pRet = xFactory.get();
666cdf0e10cSrcweir             }
667cdf0e10cSrcweir         }
668cdf0e10cSrcweir         return pRet;
669cdf0e10cSrcweir     }
670cdf0e10cSrcweir }
671