dgesv failing in windows when computing Ax=b

[arunkumar palanisamy]

When computing numerical computation for data Reconciliation, the dgesv solver fails to compute Ax=b in windows, where as for the same problem the solution is successfully computed in linux. An example to reproduce the problem is available in the testsuite in the following location

testsuite/openmodelica/dataReconciliation/TSP_FourFlows11.mos
testsuite/openmodelica/dataReconciliation/TSP_FourFlows10.mos

use the old FrontEnd to run the tests

omc -d=-newInst TSPFourFlows11.mos

I am attaching the numerical computation log for TSP_FourFlows11.mos ran from linux and windows, look for the heading in the log file

Calculations of Matrix (F*Sx*Ft) f* = c(x,y) 
============================================

From the log file we can see that the calculation made for (F*Sx*Ft) and c(x,y) is the same for both in windows and linux, but when computing f*, dgesv cannot compute the solution in windows.

[arunkumar palanisamy]

debug log from windows

[arunkumar palanisamy]

debug log from linux

[Andreas Heuermann]

I tried testsuite/openmodelica/dataReconciliation/TSP_FourFlows11.mos with todays master and for me this fails on Linux and on Windows.

record SimulationResult
    resultFile = "",
    simulationOptions = "startTime = 0.0, stopTime = 1.0, numberOfIntervals = 500, tolerance = 1e-06, method = 'dassl', fileNamePrefix = 'NewDataReconciliationSimpleTests.TSP_FourFlows11', options = '', outputFormat = 'mat', variableFilter = '.*', cflags = '', simflags = '-reconcile -sx=./NewDataReconciliationSimpleTests/resources/NewDataReconciliationSimpleTests.TSP_FourFlows9_Inputs.csv -eps=0.0023 -lv=LOG_JAC'",
    messages = "Simulation execution failed for model: NewDataReconciliationSimpleTests.TSP_FourFlows11
LOG_SUCCESS       | info    | The initialization finished successfully without homotopy method.
LOG_SUCCESS       | info    | The simulation finished successfully.
stdout            | info    | DataReconciliation Starting!
stdout            | info    | NewDataReconciliationSimpleTests.TSP_FourFlows11
stdout            | error   | solveSystemFstar() Failed !, The solution could not be computed, The info satus is 13
",
    timeFrontend = 0.6608404,
    timeBackend = 0.1515216,
    timeSimCode = 0.011052,
    timeTemplates = 0.09209050000000001,
    timeCompile = 0.900345,
    timeSimulation = 0.06415620000000001,
    timeTotal = 1.8801648
end SimulationResult;

The error message is identical as well and suggest that U(13,13) is 0, so the matrix A is singular.

So I would say it is not a problem depending on Windows or Linux.

[arunkumar palanisamy]

@Andreas, then the test (TSP_FourFlows11.mos) should fail in jenkins as well, but it is working there and also I tested in dev.openmodelica.org server which is linux Linux 4.15.0-51-generic x86_64, The baseline of TSP_FourFlows11.mos is updated using dev server

[Francesco Casella]

Retargeted to 1.18.0 because of 1.17.0 timed release.

[Francesco Casella]

Ticket retargeted after milestone closed