<html xmlns="http://www.w3.org/1999/xhtml" xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office"><head><!--[if gte mso 9]><xml><o:OfficeDocumentSettings><o:AllowPNG/><o:PixelsPerInch>96</o:PixelsPerInch></o:OfficeDocumentSettings></xml><![endif]--></head><body><div style="font-family:Helvetica Neue, Helvetica, Arial, sans-serif;font-size:13px;"><div style="font-family:Helvetica Neue, Helvetica, Arial, sans-serif;font-size:13px;"><div><div>Hi</div><div><br></div><div><div style="color: rgb(0, 0, 0); font-family: "Helvetica Neue", Helvetica, Arial, sans-serif;">You are right  only the case dim=2 has duplicates if compute second order is set tot True.<br></div></div><div><div style="color: rgb(0, 0, 0); font-family: "Helvetica Neue", Helvetica, Arial, sans-serif;">I miss it, sorry!</div><div><br></div><div>An enhancement is to generate samples of size N * (dim + 1) in case dim=2 whatever second order is True or False</div></div><div><br></div><div>Thanks for the report.</div><div><br></div><div>Regards,</div><div>Sofiane</div></div>
            <div><br></div><div><br></div>
            
            </div><div id="yahoo_quoted_3323916293" class="yahoo_quoted">
                <div style="font-family:'Helvetica Neue', Helvetica, Arial, sans-serif;font-size:13px;color:#26282a;">
                    
                    <div>
                        Le samedi 2 décembre 2017 à 13:25:25 UTC+1, Pamphile ROY <roy@cerfacs.fr> a écrit :
                    </div>
                    <div><br></div>
                    <div><br></div>
                    <div><div id="yiv3179916410"><div><div style="font-family:arial, helvetica, sans-serif;font-size:12pt;color:#000000;"><div>Hi Sofiane,</div><div><br clear="none"></div><div>I got it now.</div><div><br clear="none"></div><div>But if in 2dim this behavior is to be expected, why not doing this internally?</div><div>The root of this was that I have an expensive numerical model. So having to compute twice a sample is not tractable.<br clear="none"></div><div><br clear="none"></div><div>Thanks for your support.</div><div><br clear="none"></div><div>Sincerely,</div><div><br clear="none"></div><div>Pamphile ROY</div><hr id="yiv3179916410zwchr"><div class="yiv3179916410yqt8012359966" id="yiv3179916410yqt96824"><div><b>De: </b>"HADDAD Sofiane" <sofiane_haddad@yahoo.fr><br clear="none"><b>À: </b>"users" <users@openturns.org>, "roy" <roy@cerfacs.fr><br clear="none"><b>Envoyé: </b>Vendredi 1 Décembre 2017 13:31:14<br clear="none"><b>Objet: </b>Re: [ot-users] duplicate with SobolIndicesExperiment<br clear="none"></div><br clear="none"><div><div style="font-family:Helvetica Neue, Helvetica, Arial, sans-serif;font-size:13px;"><div><div>Hi</div><br clear="none"><div>There is no problem here</div><br clear="none"><div>You can find here how the experiment is defined. As you set second order to True and your sub samples are of size 5, you have (2 * 2 + 2) blocks of size 5 Have a look at <a rel="nofollow" shape="rect" class="yiv3179916410enhancr_card_2701847989" target="_blank" href="http://openturns.github.io/openturns/latest/user_manual/_generated/openturns.SobolIndicesAlgorithm.html?highlight=sobolindices#openturns.SobolIndicesAlgorithm">SobolIndicesAlgorithm — OpenTURNS  documentation</a></div><br clear="none"><div class="yiv3179916410ydpe1883951yahoo-link-enhancr-card yiv3179916410ydpe1883951yahoo-link-enhancr-not-allow-cover yiv3179916410ydpe1883951ymail-preserve-class yiv3179916410ydpe1883951ymail-preserve-style" id="yiv3179916410ydpe1883951enhancr_card_2701847989" style="max-width:400px;font-family:Helvetica Neue, Helvetica, Segoe UI, Arial, sans-serif;"><a rel="nofollow" shape="rect" class="yiv3179916410ydpe1883951yahoo-enhancr-cardlink" target="_blank" href="http://openturns.github.io/openturns/latest/user_manual/_generated/openturns.SobolIndicesAlgorithm.html?highlight=sobolindices#openturns.SobolIndicesAlgorithm" style="text-decoration:none;color:#000;"><table class="yiv3179916410ydpe1883951card-wrapper yiv3179916410ydpe1883951yahoo-ignore-table" border="0" cellpadding="0" cellspacing="0" style="max-width:400px;"><tbody><tr><td colspan="1" rowspan="1" width="400"><table class="yiv3179916410ydpe1883951card yiv3179916410ydpe1883951yahoo-ignore-table" border="0" cellpadding="0" cellspacing="0" width="100%" style="max-width:400px;border-radius:2px;border:1px solid #e0e4e9;"><tbody><tr><td colspan="1" rowspan="1"><table class="yiv3179916410ydpe1883951card-info yiv3179916410ydpe1883951yahoo-ignore-table" border="0" cellpadding="0" cellspacing="0" style="background:#fff;position:relative;z-index:2;width:100%;max-width:400px;border-radius:0 0 2px 2px;border-top:1px solid #e0e4e9;"><tbody><tr><td colspan="1" rowspan="1" style="background-color:#ffffff;padding:16px 0 16px 12px;vertical-align:top;border-radius:0 0 0 2px;"></td><td colspan="1" rowspan="1" style="vertical-align:middle;padding:12px 24px 16px 12px;width:99%;font-family:Helvetica Neue, Helvetica, Segoe UI, Arial, sans-serif;border-radius:0 0 2px 0;"><h2 class="yiv3179916410ydpe1883951card-title" style="font-size:14px;line-height:19px;margin:0px 0px 6px;font-family:'Helvetica Neue', Helvetica, 'Segoe UI', Arial, sans-serif;color:#26282a;">SobolIndicesAlgorithm — OpenTURNS  documentation</h2><p class="yiv3179916410ydpe1883951card-description" style="font-size:12px;line-height:16px;margin:0px;color:#979ba7;"></p></td></tr></tbody></table></td></tr></tbody></table></td></tr></tbody></table></a></div><br clear="none"><br clear="none"><br clear="none"><div>Using this :</div><div><div>import openturns as ot</div><div>ot.RandomGenerator.SetSeed(0)</div><div>distribution = ot.ComposedDistribution([ot.Uniform(15, 60), ot.Normal(3000, 400)])</div><div>sample = np.array(ot.SobolIndicesAlgorithmImplementation.Generate(distribution, 5, True))</div>print(sample.reshape(-1,5,2))</div><br clear="none"><div>you will see all matrices defines in link.</div><div>Hope this helps</div><br clear="none"><div>Bien cordialement,<br clear="none"></div><div class="yiv3179916410ydpf8c753f7signature">Sofiane HADDAD</div></div>
            <br clear="none"><br clear="none">
            
            <div class="yiv3179916410yahoo_quoted" id="yiv3179916410yahoo_quoted_2487085360">
                <div style="font-family:'Helvetica Neue', Helvetica, Arial, sans-serif;font-size:13px;color:#26282a;">
                    
                    <div>
                        Le mercredi 29 novembre 2017 à 19:30:39 UTC+1, roy <roy@cerfacs.fr> a écrit :
                    </div>
                    <br clear="none">
                    <br clear="none">
                    <div><div id="yiv3179916410"><div>Hi,<div class="yiv3179916410"><br clear="none" class="yiv3179916410"></div><div class="yiv3179916410">I am using ot.SobolIndicesExperiment and if I set my input dimension to 2,</div><div class="yiv3179916410">I get some repeated points. From my understanding of the method, this should not be the case.</div><div class="yiv3179916410">Also going with a higher dimension does not do that.</div><div class="yiv3179916410"><br clear="none" class="yiv3179916410"></div><div class="yiv3179916410">Am I wrong?</div><div class="yiv3179916410"><br clear="none" class="yiv3179916410"></div><div class="yiv3179916410">Here is an example (same behaviour with the old method available in OT 1.9) where I highlighted some duplicate:</div><div class="yiv3179916410"><br clear="none" class="yiv3179916410"></div><div class="yiv3179916410">import openturns as ot</div><div class="yiv3179916410">distribution = ot.ComposedDistribution([ot.Uniform(15, 60), ot.Normal(3000, 400)])</div><div class="yiv3179916410">sample = np.array(ot.SobolIndicesAlgorithmImplementation.Generate(distribution, 5, True))</div><div class="yiv3179916410"><br clear="none" class="yiv3179916410"></div><div class="yiv3179916410"><div class="yiv3179916410">array([[   29.65970938,  2535.47991432],</div><div class="yiv3179916410">       [   33.01991727,  2559.28624639],</div><div class="yiv3179916410">       [   33.25474751,  2682.95080229],</div><div class="yiv3179916410">       [   32.95380182,  2419.44678937],</div><div class="yiv3179916410">       [   55.23575378,  3039.33121131],</div><div class="yiv3179916410">       [   26.05095231,  3271.18330661],</div><div class="yiv3179916410">       [   41.00594229,  3683.75154513],</div><div class="yiv3179916410">       [   54.81729255,  3428.24812578],</div><div class="yiv3179916410">       [   28.2423326 ,  2797.23010815],</div><div class="yiv3179916410">       [   52.36310769,  2335.65441869],</div><div class="yiv3179916410">       [   26.05095231,  2535.47991432],</div><div class="yiv3179916410">      <b class="yiv3179916410"> [   41.00594229,  2559.28624639],</b></div><div class="yiv3179916410">       [   54.81729255,  2682.95080229],</div><div class="yiv3179916410">       [   28.2423326 ,  2419.44678937],</div><div class="yiv3179916410">       [   52.36310769,  3039.33121131],</div><div class="yiv3179916410">       [   29.65970938,  3271.18330661],</div><div class="yiv3179916410">       <b class="yiv3179916410">[   33.01991727,  3683.75154513],</b></div><div class="yiv3179916410">       [   33.25474751,  3428.24812578],</div><div class="yiv3179916410">       [   32.95380182,  2797.23010815],</div><div class="yiv3179916410">       [   55.23575378,  2335.65441869],</div><div class="yiv3179916410">       [   29.65970938,  3271.18330661],</div><div class="yiv3179916410">      <b class="yiv3179916410"> [   33.01991727,  3683.75154513],</b></div><div class="yiv3179916410">       [   33.25474751,  3428.24812578],</div><div class="yiv3179916410">       [   32.95380182,  2797.23010815],</div><div class="yiv3179916410">       [   55.23575378,  2335.65441869],</div><div class="yiv3179916410">       [   26.05095231,  2535.47991432],</div><div class="yiv3179916410"><b class="yiv3179916410">       [   41.00594229,  2559.28624639],</b></div><div class="yiv3179916410">       [   54.81729255,  2682.95080229],</div><div class="yiv3179916410">       [   28.2423326 ,  2419.44678937],</div><div class="yiv3179916410">       [   52.36310769,  3039.33121131]])</div></div><div class="yiv3179916410"><br clear="none" class="yiv3179916410"></div><div class="yiv3179916410">Thanks for your support.</div><div class="yiv3179916410"><br clear="none" class="yiv3179916410"></div><div class="yiv3179916410">Sincerely,</div><div class="yiv3179916410"><br clear="none" class="yiv3179916410"></div><br clear="none" class="yiv3179916410"><div class="yiv3179916410">
<div class="yiv3179916410" style="word-wrap:break-word;"><div style="color:#000000;font-family:Helvetica;font-size:12px;font-style:normal;letter-spacing:normal;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;"><b class="yiv3179916410">Pamphile ROY</b><br clear="none" class="yiv3179916410">PhD candidate in Uncertainty Quantification<br clear="none" class="yiv3179916410">CERFACS - Toulouse (31) - France<br clear="none" class="yiv3179916410">+33 (0) 5 61 19 31 57<br clear="none" class="yiv3179916410">+33 (0) 7 86 43 24 22</div></div>
</div>
<br clear="none" class="yiv3179916410"></div></div>_______________________________________________<br clear="none">OpenTURNS users mailing list<br clear="none"><a rel="nofollow" shape="rect" ymailto="mailto:users@openturns.org" target="_blank" href="mailto:users@openturns.org">users@openturns.org</a><br clear="none"><a rel="nofollow" shape="rect" target="_blank" href="http://openturns.org/mailman/listinfo/users">http://openturns.org/mailman/listinfo/users</a><br clear="none"></div>
                </div>
            </div></div><br clear="none"></div></div></div></div></div></div>
                </div>
            </div></div></body></html>