strnlspbc_solve() changes all x values to nan.

sgaumer — Thu, 01 Dec 2022 20:42:55 GMT

During the RCI loop, strnlspbc_solve() can suddenly change all the x parameter values to nan.

For example, below is a sequence of what happens with the (13) x values. The solver doesn't recover from this.

Is this a feature? NaN is certainly outside the bounds. Is this the solver's way of saying that it cannot converge? Is it guaranteed that all the values will be nan? Are there other pathological values that the solver can set? What's the best way to check for this condition and handle it?

Thank you.

P.S. This situation is way better than the unbounded version, strnlsp_solve(), which hangs internally and never returns.

-----------------------------------

1330.6,1142.44,10.9383,9.94877,1.77521,1.56813,0,3159.36,6.06169,5.05123,1.78558,1.70316,0
1330.6,1142.44,10.9383,9.94877,1.77521,1.56813,0,3159.36,6.06169,5.05123,1.78558,1.70316,0
1331.94,1142.44,10.9383,9.94877,1.77521,1.56813,0,3159.36,6.06169,5.05123,1.78558,1.70316,0
1329.27,1142.44,10.9383,9.94877,1.77521,1.56813,0,3159.36,6.06169,5.05123,1.78558,1.70316,0
1330.6,1143.58,10.9383,9.94877,1.77521,1.56813,0,3159.36,6.06169,5.05123,1.78558,1.70316,0
1330.6,1141.3,10.9383,9.94877,1.77521,1.56813,0,3159.36,6.06169,5.05123,1.78558,1.70316,0
1330.6,1142.44,10.9493,9.94877,1.77521,1.56813,0,3159.36,6.06169,5.05123,1.78558,1.70316,0
1330.6,1142.44,10.9274,9.94877,1.77521,1.56813,0,3159.36,6.06169,5.05123,1.78558,1.70316,0
1330.6,1142.44,10.9383,9.95872,1.77521,1.56813,0,3159.36,6.06169,5.05123,1.78558,1.70316,0
1330.6,1142.44,10.9383,9.93883,1.77521,1.56813,0,3159.36,6.06169,5.05123,1.78558,1.70316,0
1330.6,1142.44,10.9383,9.94877,1.77698,1.56813,0,3159.36,6.06169,5.05123,1.78558,1.70316,0
1330.6,1142.44,10.9383,9.94877,1.77343,1.56813,0,3159.36,6.06169,5.05123,1.78558,1.70316,0
1330.6,1142.44,10.9383,9.94877,1.77521,1.5697,0,3159.36,6.06169,5.05123,1.78558,1.70316,0
1330.6,1142.44,10.9383,9.94877,1.77521,1.56656,0,3159.36,6.06169,5.05123,1.78558,1.70316,0
1330.6,1142.44,10.9383,9.94877,1.77521,1.56813,0.001,3159.36,6.06169,5.05123,1.78558,1.70316,0
1330.6,1142.44,10.9383,9.94877,1.77521,1.56813,-0.001,3159.36,6.06169,5.05123,1.78558,1.70316,0
1330.6,1142.44,10.9383,9.94877,1.77521,1.56813,0,3162.52,6.06169,5.05123,1.78558,1.70316,0
1330.6,1142.44,10.9383,9.94877,1.77521,1.56813,0,3156.2,6.06169,5.05123,1.78558,1.70316,0
1330.6,1142.44,10.9383,9.94877,1.77521,1.56813,0,3159.36,6.06775,5.05123,1.78558,1.70316,0
1330.6,1142.44,10.9383,9.94877,1.77521,1.56813,0,3159.36,6.05563,5.05123,1.78558,1.70316,0
1330.6,1142.44,10.9383,9.94877,1.77521,1.56813,0,3159.36,6.06169,5.05628,1.78558,1.70316,0
1330.6,1142.44,10.9383,9.94877,1.77521,1.56813,0,3159.36,6.06169,5.04618,1.78558,1.70316,0
1330.6,1142.44,10.9383,9.94877,1.77521,1.56813,0,3159.36,6.06169,5.05123,1.78737,1.70316,0
1330.6,1142.44,10.9383,9.94877,1.77521,1.56813,0,3159.36,6.06169,5.05123,1.78379,1.70316,0
1330.6,1142.44,10.9383,9.94877,1.77521,1.56813,0,3159.36,6.06169,5.05123,1.78558,1.70486,0
1330.6,1142.44,10.9383,9.94877,1.77521,1.56813,0,3159.36,6.06169,5.05123,1.78558,1.70146,0
1330.6,1142.44,10.9383,9.94877,1.77521,1.56813,0,3159.36,6.06169,5.05123,1.78558,1.70316,0.001
1330.6,1142.44,10.9383,9.94877,1.77521,1.56813,0,3159.36,6.06169,5.05123,1.78558,1.70316,-0.001
1070.41,324.814,9.72218,9.39311,2.02915,1.51927,0.227844,1904.56,6.4404,4.79253,2.37283,2.31391,0.0743298
1071.48,324.814,9.72218,9.39311,2.02915,1.51927,0.227844,1904.56,6.4404,4.79253,2.37283,2.31391,0.0743298
1069.34,324.814,9.72218,9.39311,2.02915,1.51927,0.227844,1904.56,6.4404,4.79253,2.37283,2.31391,0.0743298
1070.41,325.139,9.72218,9.39311,2.02915,1.51927,0.227844,1904.56,6.4404,4.79253,2.37283,2.31391,0.0743298
1070.41,324.489,9.72218,9.39311,2.02915,1.51927,0.227844,1904.56,6.4404,4.79253,2.37283,2.31391,0.0743298
1070.41,324.814,9.7319,9.39311,2.02915,1.51927,0.227844,1904.56,6.4404,4.79253,2.37283,2.31391,0.0743298
1070.41,324.814,9.71246,9.39311,2.02915,1.51927,0.227844,1904.56,6.4404,4.79253,2.37283,2.31391,0.0743298
1070.41,324.814,9.72218,9.4025,2.02915,1.51927,0.227844,1904.56,6.4404,4.79253,2.37283,2.31391,0.0743298
1070.41,324.814,9.72218,9.38372,2.02915,1.51927,0.227844,1904.56,6.4404,4.79253,2.37283,2.31391,0.0743298
1070.41,324.814,9.72218,9.39311,2.03118,1.51927,0.227844,1904.56,6.4404,4.79253,2.37283,2.31391,0.0743298
1070.41,324.814,9.72218,9.39311,2.02712,1.51927,0.227844,1904.56,6.4404,4.79253,2.37283,2.31391,0.0743298
1070.41,324.814,9.72218,9.39311,2.02915,1.52079,0.227844,1904.56,6.4404,4.79253,2.37283,2.31391,0.0743298
1070.41,324.814,9.72218,9.39311,2.02915,1.51775,0.227844,1904.56,6.4404,4.79253,2.37283,2.31391,0.0743298
1070.41,324.814,9.72218,9.39311,2.02915,1.51927,0.228072,1904.56,6.4404,4.79253,2.37283,2.31391,0.0743298
1070.41,324.814,9.72218,9.39311,2.02915,1.51927,0.227616,1904.56,6.4404,4.79253,2.37283,2.31391,0.0743298
1070.41,324.814,9.72218,9.39311,2.02915,1.51927,0.227844,1906.46,6.4404,4.79253,2.37283,2.31391,0.0743298
1070.41,324.814,9.72218,9.39311,2.02915,1.51927,0.227844,1902.65,6.4404,4.79253,2.37283,2.31391,0.0743298
1070.41,324.814,9.72218,9.39311,2.02915,1.51927,0.227844,1904.56,6.44685,4.79253,2.37283,2.31391,0.0743298
1070.41,324.814,9.72218,9.39311,2.02915,1.51927,0.227844,1904.56,6.43396,4.79253,2.37283,2.31391,0.0743298
1070.41,324.814,9.72218,9.39311,2.02915,1.51927,0.227844,1904.56,6.4404,4.79733,2.37283,2.31391,0.0743298
1070.41,324.814,9.72218,9.39311,2.02915,1.51927,0.227844,1904.56,6.4404,4.78774,2.37283,2.31391,0.0743298
1070.41,324.814,9.72218,9.39311,2.02915,1.51927,0.227844,1904.56,6.4404,4.79253,2.3752,2.31391,0.0743298
1070.41,324.814,9.72218,9.39311,2.02915,1.51927,0.227844,1904.56,6.4404,4.79253,2.37045,2.31391,0.0743298
1070.41,324.814,9.72218,9.39311,2.02915,1.51927,0.227844,1904.56,6.4404,4.79253,2.37283,2.31622,0.0743298
1070.41,324.814,9.72218,9.39311,2.02915,1.51927,0.227844,1904.56,6.4404,4.79253,2.37283,2.31159,0.0743298
1070.41,324.814,9.72218,9.39311,2.02915,1.51927,0.227844,1904.56,6.4404,4.79253,2.37283,2.31391,0.0744041
1070.41,324.814,9.72218,9.39311,2.02915,1.51927,0.227844,1904.56,6.4404,4.79253,2.37283,2.31391,0.0742555
986.475,376.448,11.7741,10.8572,1.06561,0.941168,-0.999386,2246.47,6.43891,5.11562,2.5202,3.40577,0.269949
987.461,376.448,11.7741,10.8572,1.06561,0.941168,-0.999386,2246.47,6.43891,5.11562,2.5202,3.40577,0.269949
985.488,376.448,11.7741,10.8572,1.06561,0.941168,-0.999386,2246.47,6.43891,5.11562,2.5202,3.40577,0.269949
986.475,376.824,11.7741,10.8572,1.06561,0.941168,-0.999386,2246.47,6.43891,5.11562,2.5202,3.40577,0.269949
986.475,376.072,11.7741,10.8572,1.06561,0.941168,-0.999386,2246.47,6.43891,5.11562,2.5202,3.40577,0.269949
986.475,376.448,11.7859,10.8572,1.06561,0.941168,-0.999386,2246.47,6.43891,5.11562,2.5202,3.40577,0.269949
986.475,376.448,11.7623,10.8572,1.06561,0.941168,-0.999386,2246.47,6.43891,5.11562,2.5202,3.40577,0.269949
986.475,376.448,11.7741,10.8681,1.06561,0.941168,-0.999386,2246.47,6.43891,5.11562,2.5202,3.40577,0.269949
986.475,376.448,11.7741,10.8463,1.06561,0.941168,-0.999386,2246.47,6.43891,5.11562,2.5202,3.40577,0.269949
986.475,376.448,11.7741,10.8572,1.06667,0.941168,-0.999386,2246.47,6.43891,5.11562,2.5202,3.40577,0.269949
986.475,376.448,11.7741,10.8572,1.06454,0.941168,-0.999386,2246.47,6.43891,5.11562,2.5202,3.40577,0.269949
986.475,376.448,11.7741,10.8572,1.06561,0.942109,-0.999386,2246.47,6.43891,5.11562,2.5202,3.40577,0.269949
986.475,376.448,11.7741,10.8572,1.06561,0.940226,-0.999386,2246.47,6.43891,5.11562,2.5202,3.40577,0.269949
986.475,376.448,11.7741,10.8572,1.06561,0.941168,-1.00039,2246.47,6.43891,5.11562,2.5202,3.40577,0.269949
986.475,376.448,11.7741,10.8572,1.06561,0.941168,-0.998387,2246.47,6.43891,5.11562,2.5202,3.40577,0.269949
986.475,376.448,11.7741,10.8572,1.06561,0.941168,-0.999386,2248.72,6.43891,5.11562,2.5202,3.40577,0.269949
986.475,376.448,11.7741,10.8572,1.06561,0.941168,-0.999386,2244.22,6.43891,5.11562,2.5202,3.40577,0.269949
986.475,376.448,11.7741,10.8572,1.06561,0.941168,-0.999386,2246.47,6.44535,5.11562,2.5202,3.40577,0.269949
986.475,376.448,11.7741,10.8572,1.06561,0.941168,-0.999386,2246.47,6.43247,5.11562,2.5202,3.40577,0.269949
986.475,376.448,11.7741,10.8572,1.06561,0.941168,-0.999386,2246.47,6.43891,5.12074,2.5202,3.40577,0.269949
986.475,376.448,11.7741,10.8572,1.06561,0.941168,-0.999386,2246.47,6.43891,5.11051,2.5202,3.40577,0.269949
986.475,376.448,11.7741,10.8572,1.06561,0.941168,-0.999386,2246.47,6.43891,5.11562,2.52272,3.40577,0.269949
986.475,376.448,11.7741,10.8572,1.06561,0.941168,-0.999386,2246.47,6.43891,5.11562,2.51768,3.40577,0.269949
986.475,376.448,11.7741,10.8572,1.06561,0.941168,-0.999386,2246.47,6.43891,5.11562,2.5202,3.40918,0.269949
986.475,376.448,11.7741,10.8572,1.06561,0.941168,-0.999386,2246.47,6.43891,5.11562,2.5202,3.40237,0.269949
986.475,376.448,11.7741,10.8572,1.06561,0.941168,-0.999386,2246.47,6.43891,5.11562,2.5202,3.40577,0.270219
986.475,376.448,11.7741,10.8572,1.06561,0.941168,-0.999386,2246.47,6.43891,5.11562,2.5202,3.40577,0.269679
nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan,nan

Re:strnlspbc_solve() changes all x values to nan.

VidyalathaB_Intel — Fri, 02 Dec 2022 12:24:23 GMT

Hi Scott,

Thanks for reaching out to us.

Could you please provide us with a sample reproducer code along with the steps to reproduce the issue so that we can test it from our end?

Please do let us know the OS environment details and MKL version you are using in this case.

Regards,

Vidya.

Re: Re:strnlspbc_solve() changes all x values to nan.

sgaumer — Wed, 07 Dec 2022 01:45:40 GMT

It's Windows 10 and 2022.2.1.
Is there a way I can share the code privately?

Re:strnlspbc_solve() changes all x values to nan.

VidyalathaB_Intel — Wed, 07 Dec 2022 06:12:08 GMT

Hi Scott,

>>Is there a way I can share the code privately?

Thanks for getting back to us.

Yeah sure, we will contact you privately and you can share your code there.

Regards,

Vidya.

Re:strnlspbc_solve() changes all x values to nan.

VidyalathaB_Intel — Wed, 14 Dec 2022 12:18:41 GMT

Hi Scott,

We have already contacted you privately but we don't see any response there. As a reminder note could you please provide us with an update regarding the issue?

Regards,

Vidya.

Re:strnlspbc_solve() changes all x values to nan.

VidyalathaB_Intel — Wed, 21 Dec 2022 06:39:26 GMT

Hi Scott,

As we haven't heard back from you, we are closing this thread. Please post a new question if you need any additional assistance from Intel as this thread will no longer be monitored.

Regards,

Vidya.

Re: Re:strnlspbc_solve() changes all x values to nan.

sgaumer — Thu, 12 Jan 2023 17:19:46 GMT

Vidya,

I replied to the email and they said "Please note that the reproducer sample was asked by our representative from the
engineering team of oneMKL. Since we do not handle the oneAPI Math Kernel Library
(oneMKL) issues, we suggest you to post your query directly to the relevant team by
clicking here: Intel® oneAPI Math Kernel Library...".

Can you provide a way I can share the reproducible code examples with you? Otherwise, perhaps you could look at 05710106 and 05711432.

Thank you.

Re: Re:strnlspbc_solve() changes all x values to nan.

VidyalathaB_Intel — Mon, 16 Jan 2023 09:55:38 GMT

Hi Scott,

Thanks for getting back to us, and I apologize that we didn't reply to you earlier.

Since this specific issue is archived, I would recommend you create a new thread (here in the MKL forum) so that we can resume working on this issue.

Regards,

Vidya.

topic strnlspbc_solve() changes all x values to nan. in Intel® oneAPI Math Kernel Library

strnlspbc_solve() changes all x values to nan.

Re:strnlspbc_solve() changes all x values to nan.

Re: Re:strnlspbc_solve() changes all x values to nan.

Re:strnlspbc_solve() changes all x values to nan.

Re:strnlspbc_solve() changes all x values to nan.

Re:strnlspbc_solve() changes all x values to nan.

Re: Re:strnlspbc_solve() changes all x values to nan.

Re: Re:strnlspbc_solve() changes all x values to nan.