"1仮想球の表面の上に多数の点をランダムに配置するステップ、2上記点に基づいて、仮想球の表面の上に多数の領域を想定するステップ、3上記領域の輪郭に基づいて、仮想球の表面にディンプルとランドとを割り当てるステップを含む、ゴルフボール表面の凹凸パターン設計方法。上記ステップ2において、上記領域の1つに上記点が1つ含まれるように、仮想球の表面の上に多数の領域が想定される請求項1に記載の設計方法。上記ステップ2において、ボロノイ分割によって領域が想定される請求項2に記載の設計方法。上記ステップ2が、2.1上記仮想球の表面の上に多数の微小なセルを想定するステップ、2.2それぞれのセルに最も近い点を母点として選定するステップ、2.3それぞれの母点に関し、この母点を最も近い母点とするセルの集合を想定するステップ、及び2.4それぞれの集合を、ボロノイ分割によって得られる領域とみなすステップを含む請求項3に記載の設計方法。上記ステップ2において、1つの点と他の点とを線で結ぶことにより、仮想球の表面の上に多数の領域が想定される請求項1に記載の設計方法。上記ステップ2において、ドロネー三角分割によって領域が想定される請求項5に記載の設計方法。上記ステップ3において、それぞれの輪郭の近傍にランドが割り当てられ、仮想球の表面のうちランド以外の部分にディンプルが割り当てられる請求項1から6のいずれかに記載の設計方法。上記ステップ1において、セルオートマトン法に基づいて多数の点がランダムに配置される請求項1から7のいずれかに記載の設計方法。上記ステップ1が、1.1複数の状態が想定されるステップ、1.2仮想球の表面の上に多数のセルが想定されるステップ、1.3それぞれのセルに、いずれかの状態が付与されるステップ、1.4上記セルの状態及びこのセルの近傍に位置する複数のセルの状態に基づいて、当該セルの属性として、インサイド、アウトサイド及び境界のいずれかが付与されるステップ、1.5境界のセルが結ばれることによって、ループが想定されるステップ、並びに1.6上記ループの中心点が決定されるステップを含む請求項8に記載の設計方法。その表面に多数のディンプルを備えており、これらのディンプルがランダムに配置されており、それぞれのディンプルが、他のいずれのディンプルとも異なる輪郭形状を有しているゴルフボール。上記ディンプルのパターンが請求項1から9のいずれかに記載の方法で設計された請求項10に記載のゴルフボール。下記ステップ1から16によって得られる変動幅Rh及び変動幅Roが、2.7mm以下である請求項10又は11に記載のゴルフボール。1ゴルフボールの両極を結ぶ線が、第一回転軸に想定されるステップ2ゴルフボールの仮想球の表面に存在し、かつ上記第一回転軸と直交する大円が想定されるステップ3ゴルフボールの仮想球の表面に存在し、上記第一回転軸と直交し、かつ上記大円との中心角の絶対値が30°である2つの小円が想定されるステップ4これらの小円によりゴルフボールが区画され、このゴルフボールの表面のうちこれら小円に挟まれた領域が特定されるステップ5上記領域に、軸方向において中心角度で3°刻みであり回転方向において中心角で0.25°刻みに、30240の点が決定されるステップ6それぞれの点から上記第一回転軸に下ろした垂線の長さL1が算出されるステップ7軸方向に並ぶ21個の垂線に基づいて算出された21個の長さL1が合計され、総長さL2が算出されるステップ8回転方向に沿って算出される1440個の総長さL2から、最大値と最小値とが決定され、最大値から最小値が減じられた値である変動幅Rhが算出されるステップ9上記ステップ1で想定された第一回転軸に直交する第二回転軸が想定されるステップ10ゴルフボールの仮想球の表面に存在し、かつ上記第二回転軸と直交する大円が想定されるステップ11ゴルフボールの仮想球の表面に存在し、上記第二回転軸と直交し、かつ上記大円との中心角の絶対値が30°である2つの小円が想定されるステップ12これらの小円によりゴルフボールが区画され、ゴルフボールの表面のうちこれら小円に挟まれた領域が特定されるステップ13上記領域に、軸方向において中心角度で3°刻みであり回転方向において中心角で0.25°刻みに、30240の点が決定されるステップ14それぞれの点から上記第二回転軸に下ろした垂線の長さL1が算出されるステップ15軸方向に並ぶ21個の垂線に基づいて算出された21個の長さL1が合計され、総長さL2が算出されるステップ16回転方向に沿って算出される1440個の総長さL2から、最大値と最小値とが決定され、最大値から最小値が減じられた値である変動幅Roが算出されるステップ上記変動幅Rhと上記変動幅Roとの差dRの絶対値が0.4mm以下である請求項12に記載のゴルフボール。特開平9-164223号公報に開示された方法では、ディンプルパターンを得るのに、試行錯誤が繰り返される。特開2000-189542公報に開示された方法でも、ディンプルパターンを得るのに、試行錯誤が繰り返される。本発明の目的は、空力的対称性に優れたゴルフボールの提供にある。本発明に係るゴルフボール表面の凹凸パターン設計方法は、1仮想球の表面の上に多数の点をランダムに配置するステップ、2上記点に基づいて、仮想球の表面の上に多数の領域を想定するステップ、3上記領域の輪郭に基づいて、仮想球の表面にディンプルとランドとを割り当てるステップを含む。好ましくは、ステップ2において、上記領域の1つに上記点が1つ含まれるように、仮想球の表面の上に多数の領域が想定される。好ましくは、ステップ2において、ボロノイ分割によって領域が想定される。好ましくは、ステップ2は、2.1上記仮想球の表面の上に多数の微小なセルを想定するステップ、2.2それぞれのセルに最も近い点を母点として選定するステップ、2.3それぞれの母点に関し、この母点を最も近い母点とするセルの集合を想定するステップ、及び2.4それぞれの集合を、ボロノイ分割によって得られる領域とみなすステップを含む。ステップ2において、1つの点と他の点とを線で結ぶことにより、仮想球の表面の上に多数の領域が想定されてもよい。好ましくは、ステップ2において、ドロネー三角分割によって領域が想定される。好ましくは、ステップ3において、それぞれの輪郭の近傍にランドが割り当てられ、仮想球の表面のうちランド以外の部分にディンプルが割り当てられる。好ましくは、上記ステップ1において、セルオートマトン法に基づいて多数の点がランダムに配置される。好ましくは、ステップ1は、1.1複数の状態が想定されるステップ、1.2仮想球の表面の上に多数のセルが想定されるステップ、1.3それぞれのセルに、いずれかの状態が付与されるステップ、1.4上記セルの状態及びこのセルの近傍に位置する複数のセルの状態に基づいて、当該セルの属性として、インサイド、アウトサイド及び境界のいずれかが付与されるステップ、1.5境界のセルが結ばれることによって、ループが想定されるステップ、並びに1.6上記ループの中心点が決定されるステップを含む。本発明に係るゴルフボールは、その表面に多数のディンプルを備える。これらのディンプルは、ランダムに配置されている。それぞれのディンプルは、他のいずれのディンプルとも異なる輪郭形状を有している。好ましくは、下記ステップ1から16によって得られる変動幅Rh及び変動幅Roは、2.7mm以下である。1ゴルフボールの両極を結ぶ線が、第一回転軸に想定されるステップ2ゴルフボールの仮想球の表面に存在し、かつ上記第一回転軸と直交する大円が想定されるステップ3ゴルフボールの仮想球の表面に存在し、上記第一回転軸と直交し、かつ上記大円との中心角の絶対値が30°である2つの小円が想定されるステップ4これらの小円によりゴルフボールが区画され、このゴルフボールの表面のうちこれら小円に挟まれた領域が特定されるステップ5上記領域に、軸方向において中心角度で3°刻みであり回転方向において中心角で0.25°刻みに、30240の点が決定されるステップ6それぞれの点から上記第一回転軸に下ろした垂線の長さL1が算出されるステップ7軸方向に並ぶ21個の垂線に基づいて算出された21個の長さL1が合計され、総長さL2が算出されるステップ8回転方向に沿って算出される1440個の総長さL2から、最大値と最小値とが決定され、最大値から最小値が減じられた値である変動幅Rhが算出されるステップ9上記ステップ1で想定された第一回転軸に直交する第二回転軸が想定されるステップ10ゴルフボールの仮想球の表面に存在し、かつ上記第二回転軸と直交する大円が想定されるステップ11ゴルフボールの仮想球の表面に存在し、上記第二回転軸と直交し、かつ上記大円との中心角の絶対値が30°である2つの小円が想定されるステップ12これらの小円によりゴルフボールが区画され、ゴルフボールの表面のうちこれら小円に挟まれた領域が特定されるステップ13上記領域に、軸方向において中心角度で3°刻みであり回転方向において中心角で0.25°刻みに、30240の点が決定されるステップ14それぞれの点から上記第二回転軸に下ろした垂線の長さL1が算出されるステップ15軸方向に並ぶ21個の垂線に基づいて算出された21個の長さL1が合計され、総長さL2が算出されるステップ16回転方向に沿って算出される1440個の総長さL2から、最大値と最小値とが決定され、最大値から最小値が減じられた値である変動幅Roが算出されるステップ好ましくは、変動幅Rhと変動幅Roとの差dRの絶対値は、0.4mm以下である。本発明に係る設計方法により、空力的対称性に優れたゴルフボールが容易に得られうる。以下、適宜図面が参照されつつ、好ましい実施形態に基づいて本発明が詳細に説明される。図1に示されたゴルフボール2は、球状のコア4と、カバー6とを備えている。カバー6の表面には、多数のディンプル8が形成されている。ゴルフボール2の表面のうちディンプル8以外の部分は、ランド10である。このゴルフボール2は、カバー6の外側にペイント層及びマーク層を備えているが、これらの層の図示は省略されている。コア4とカバー6との間に、中間層が設けられてもよい。このゴルフボール2の直径は、40mm以上45mm以下が好ましい。米国ゴルフ協会USGAの規格が満たされるとの観点から、直径は42.67mm以上が特に好ましい。空気抵抗抑制の観点から、直径は44mm以下がより好ましく、42.80mm以下が特に好ましい。このゴルフボール2の質量は、40g以上50g以下が好ましい。大きな慣性が得られるとの観点から、質量は44g以上がより好ましく、45.00g以上が特に好ましい。USGAの規格が満たされるとの観点から、質量は45.93g以下が特に好ましい。コア4は、ゴム組成物が架橋されることによって形成されている。ゴム組成物の基材ゴムとしては、ポリブタジエン、ポリイソプレン、スチレン-ブタジエン共重合体、エチレン-プロピレン-ジエン共重合体及び天然ゴムが例示される。2種以上のゴムが併用されてもよい。反発性能の観点からポリブタジエンが好ましく、特にハイシスポリブタジエンが好ましい。コア4の架橋には、共架橋剤が用いられうる。反発性能の観点から好ましい共架橋剤は、アクリル酸亜鉛、アクリル酸マグネシウム、メタクリル酸亜鉛及びメタクリル酸マグネシウムである。ゴム組成物には、共架橋剤と共に有機過酸化物が配合されるのが好ましい。好適な有機過酸化物としては、ジクミルパーオキサイド、1,1-ビスt-ブチルパーオキシ-3,3,5-トリメチルシクロヘキサン、2,5-ジメチル-2,5-ジt-ブチルパーオキシヘキサン及びジ-t-ブチルパーオキサイドが挙げられる。コア4のゴム組成物には、硫黄、硫黄化合物、充填剤、老化防止剤、着色剤、可塑剤、分散剤等の各種添加剤が、必要に応じて適量配合される。ゴム組成物に、架橋ゴム粉末又は合成樹脂粉末が配合されてもよい。コア4の直径は30.0mm以上、特には38.0mm以上である。コア4の直径は42.0mm以下、特には41.5mm以下である。コア4が2以上の層から構成されてもよい。コア4がその表面にリブを備えてもよい。カバー6に好適なポリマーは、アイオノマー樹脂である。好ましいアイオノマー樹脂としては、α-オレフィンと炭素数が3以上8以下のα,β-不飽和カルボン酸との二元共重合体が挙げられる。好ましい他のアイオノマー樹脂としては、α-オレフィンと炭素数が3以上8以下のα,β-不飽和カルボン酸と炭素数が2以上22以下のα,β-不飽和カルボン酸エステルとの三元共重合体が挙げられる。二元共重合体及び三元共重合体において、好ましいα-オレフィンはエチレン及びプロピレンであり、好ましいα,β-不飽和カルボン酸はアクリル酸及びメタクリル酸である。二元共重合体及び三元共重合体において、カルボキシル基の一部は金属イオンで中和されている。中和のための金属イオンとしては、ナトリウムイオン、カリウムイオン、リチウムイオン、亜鉛イオン、カルシウムイオン、マグネシウムイオン、アルミニウムイオン及びネオジムイオンが例示される。アイオノマー樹脂に代えて、又はアイオノマー樹脂と共に、他のポリマーが用いられてもよい。他のポリマーとしては、熱可塑性ポリウレタンエラストマー、熱可塑性スチレンエラストマー、熱可塑性ポリアミドエラストマー、熱可塑性ポリエステルエラストマー及び熱可塑性ポリオレフィンエラストマーが例示される。スピン性能の観点から、熱可塑性ポリウレタンエラストマーが好ましい。カバー6には、必要に応じ、二酸化チタンのような着色剤、硫酸バリウムのような充填剤、分散剤、酸化防止剤、紫外線吸収剤、光安定剤、蛍光剤、蛍光増白剤等が適量配合される。比重調整の目的で、カバー6にタングステン、モリブデン等の高比重金属の粉末が配合されてもよい。カバー6の厚みは0.1mm以上、特には0.3mm以上である。カバー6の厚みは2.5mm以下、特には2.2mm以下である。カバー6の比重は0.90以上、特には0.95以上である。カバー6の比重は1.10以下、特には1.05以下である。カバー6が2以上の層から構成されてもよい。図2は、このゴルフボール2が示された拡大正面図である。図3は、図2のゴルフボール2が示された平面図である。図2及び3から明らかなように、多数のディンプル8がランダムに配置されている。これらディンプル8とランド10とにより、ゴルフボール2の表面に凹凸パターンが形成されている。この凹凸パターンの設計方法には、ボロノイ分割Voronoitessellationが用いられる。この設計方法では、多数の母点が仮想球14図1参照の表面の上にランダムに配置される。この母点に基づいて、ボロノイ分割により、仮想球14の表面の上に多数の領域が想定される。本願明細書では、この領域は、ボロノイ領域と称される。このボロノイ領域に基づいて、ディンプル8が割り当てられる。この設計方法は、効率の観点から、コンピュータとソフトウエアとが用いられて実施されることが好ましい。もちろん、手計算でも本発明は実施されうる。本発明の本質がコンピュータソフトウエアにあるわけではない。好ましくは、母点の配置には、セルオートマトンCellularAutomaton法が用いられる。セルオートマトン法により、仮想球14の表面の上に多数のループがランダムに配置されたパターンが得られる。これらループの中心点が、母点とされる。ループの配置がランダムなので、母点の配置もランダムである。セルオートマトン法は、計算可能性理論、数学、理論生物学等の分野で、広く利用されている。セルオートマトン法のモデルは、多数のセルと単純な規則とからなる。このモデルにより、生命現象、結晶の成長、乱流等の自然現象が模擬されうる。このモデルでは、それぞれのセルが状態を有する。この状態は、ステージの進行に応じ、他の状態に変化しうる。あるセルのステージt1における状態は、ステージtにおける当該セルの状態及びこのセルの近傍にある複数のセルの状態によって決定される。決定は、ある規則に従ってなされる。全てのセルに、この規則が等しく適用される。この凹凸パターンの設計には、セルオートマトン法の反応拡散モデルが適している。このモデルは、獣、鳥、魚、昆虫等の体表面のループの模擬に用いられている。このモデルでは、複数の状態が想定される。状態の数は、通常は2以上8以下である。それぞれのセルにおいて、初期の状態が決定される。ステージの進行により、規則に基づいて、状態が更新される。更新により、その状態が変化するセルが存在する。更新により、状態が変化しないセルも存在する。セルオートマトン法は、セルオートマトン法複雑系の自己組織化と超並列処理加藤泰義ら著、森北出版株式会社発行の第25-28頁に開示されている。本発明に係る設計方法の特徴は、セルの状態が、当該セルの近傍に位置する他のセルの影響下にて更新される点にある。この更新により、多数のループがランダムに配置されたパターンが得られる。この特徴が維持される限り、如何なるモデルも用いられうる。以下、セルオートマトン法の反応拡散モデルを用いた設計方法が詳説される。図4は、ループのパターンの設計方法が示されたフローチャートである。図5は、図4の設計方法に用いられるメッシュ12が示された正面図である。このメッシュ12の形成には、球14が仮想されるSTEP1。この仮想球14の直径は、ゴルフボール2の直径と同一である。この仮想球14の表面が、多数の三角形に分割されるSTEP2。分割は、前進先端法advancingfrontmethodに基づいてなされている。前進先端法が、大学院情報理工学3計算力学伊藤耿一編、講談社発行の第195-197頁に開示されている。このメッシュ12において、三角形の数は176528個であり、頂点の数は88266個である。それぞれの頂点は、セル又はセルの中心と定義される。このメッシュ12では、セルの数は88266個である。他の手法によって仮想球14が分割されてもよい。この設計方法では、分化及び未分化の、2つの状態が想定される。それぞれのセルにおいて、いずれかの状態初期の状態が決定されるSTEP3。決定は、好ましくは、無作為になされる。無作為な決定には、乱数と剰余系とが用いられる。状態の数が2なので、基数が2である剰余系が用いられる。具体的には、コンピュータによって0以上1未満であり、下5桁の乱数が発生させられる。この乱数に100000が乗され、この積がさらに2で除される。この商の余りは、1又は0である。この余りに基づいて、当該セルの状態が決定される。例えば、余りが1の場合に分化が選定され、余りが0の場合に未分化が選定される。全てのセルについて、決定がなされる。決定の後のメッシュ12は、ステージ1にある。それぞれのセルにおいて、状態の変更の要否が判定されるSTEP4。判定は、規則に従ってなされる。図6は、この規則が説明されるためのグラフである。このグラフにおいて、縦軸は濃度であり、横軸は指数半径である。指数半径は、当該セルからの距離が基準値で除された値である。この基準値は、当該セルに最も近いセルと、当該セルとの距離である。濃度W1は正であり、濃度W2は負である。濃度W1の絶対値は、濃度W2の絶対値よりも大きい。指数半径R2は、指数半径R1よりも大きい。指数半径が0を超えてR1以下のエリアでは、濃度はW1である。指数半径がR1を超えてR2以下のエリアでは、濃度はW2である。図7は、図5のメッシュ12の一部が示された拡大図である。便宜上、図7では、メッシュ12が二次元で画かれている。図7の中心には、判定の対象であるセル16aが示されている。図7にはさらに、第一円18及び第二円20が示されている。第一円18は、セル16aを中心とする、指数半径がR1である円である。第二円20は、セル16aを中心とする、指数半径がR2である円である。塗りつぶされた円で示されているのは、第一円18に含まれる、セル16a以外のセル16である。塗りつぶされた四角形で示されているのは、第二円20に含まれ、第一円18に含まれないセル16である。塗りつぶされた三角形で示されているのは、第二円20に含まれないセル16である。この設計方法では、第一円18に含まれこの第一円18の中心に位置しない、特定の状態にあるセル16の数NR1がカウントされる。好ましい態様によれば、その状態が分化であるセル16の数がカウントされ、合計NR1が算出される。この設計方法ではさらに、第二円20に含まれ第一円18に含まれない、特定の状態にあるセル16の数NR1-R2がカウントされる。好ましい態様によれば、その状態が分化であるセル16の数がカウントされ、合計NR1-R2が算出される。この数NR1及びNR1-R2が下記数式1に代入され、値Eが算出される。この値Eに基づいて、セル16aの状態の変更の要否が決定される。E=W1NR1W2NR1-R21この判定に基づき、セル16aの状態の更新がなされるSTEP5。更新において、セル16a状態が変化することも、変化しないことも生じうる。好ましい態様によれば、値Eが正であるとき、セル16aの状態が分化であればこの状態が維持され、セル16aの状態が未分化であればこの状態が分化に変更される。この値Eがゼロであるとき、セル16aの状態が維持される。この値Eが負であるとき、セル16aの状態が分化であればこの状態が未分化に変更され、セル16aの状態が未分化であればこの状態が維持される。全てのセル16についての第一回目の更新が完了したメッシュ12は、ステージ2にある。以下に、判定及び更新の計算例が示される。条件第一濃度W1:1.00第二濃度W2:-0.60第一円18に含まれその状態が分化であるセル16の数当該セル16aを除く:8第二円20に含まれ第一円18に含まれない、その状態が分化であるセル16の数:13計算例E=1.008-0.6013=0.2この場合、値Eが正なので、セル16aの状態が分化であればこの状態が維持され、セル16aの状態が未分化であればこの状態が分化に変更される。以下に、判定及び更新の他の計算例が示される。条件第一濃度W1:1.00第二濃度W2:-0.60第一円18に含まれその状態が分化であるセル16の数:5第二円20に含まれ第一円18に含まれない、その状態が分化であるセル16の数:9計算例E=1.005-0.609=-0.4この場合、値Eが負なので、セル16aの状態が分化であればこの状態が未分化に変更され、セル16aの状態が未分化であればこの状態が維持される。この判定と更新とが、繰り返される。図4のフローチャートでは、繰り返し数はMである。M回の繰り返しが完了した後のメッシュ12は、ステージM1にある。ステージの進行に伴い、更新によって状態が変化するセル16の数が減少する。繰り返し数が小さいステージでは、更新によるパーターンの変化が激しい。多数回の更新がなされることで、パターンが収束する。繰り返し数は3以上が好ましく、5以上がより好ましい。繰り返し数が過大であると、コンピュータへの負荷が大きい。この観点から、繰り返し数は30以下が好ましく、10以下がより好ましい。判定と更新とがM回繰り返されることにより、それぞれのセル16の状態が確定する。この確定は、セル16への状態の付与である。図8は、状態の付与が完了した後のメッシュ12の一部が示された拡大図である。図8において、円で示されているのは分化のセル16であり、四角形で示されているのは未分化のセル16である。この状態に基づき、セル16にiflagが付与される。まず暫定的に、全てのセル16にiflagとして0が付与される。次に、状態が分化であるセル16に関し、iflagが変更される。図8において符号16bで示されたセル16は、6個のセル16c-16hと隣接している。本発明では、一方のセル16を頂点とする三角形の他の頂点に他方のセル16が存在するとき、一方のセル16が他方のセル16と隣接すると称される。これらのセル16c-16hの状態は、分化である。隣接する全てのセル16c-16hの状態が分化であるとき、当該セル16bのiflagが0から1に変更される。図8において符号16nで示されたセル16は、6個のセル16h-16mと隣接している。セル16h、16i、16l、16mの状態は、分化である。セル16j、16kの状態は、未分化である。状態が未分化である1又は2以上のセル16と隣接するとき、当該セル16nのiflagが0から2に変更される。状態が分化である全てのセル16に関し、そのiflagが変更される。状態が未分化であるセル16のiflagは、変更されない。このiflagに基づき、全てのセル16に属性が付与されるSTEP6。属性の付与は、下記のルールに基づいてなされる。iflag:0属性:アウトサイドiflag:1属性:インサイドiflag:2属性:境界属性の付与が完了したメッシュ12は、第一フェーズにある。属性が境界である複数のセル16を結ぶことにより、第一ループ21が完成する。図8において、第一ループ21が太線で示されている。多数の第一ループ21を有するパターンが、図9に示されている。このパターンは、以下のパラメータが用いられて得られたものである。W1:1.0W2:-0.6R1:4.5R2:8.0このパターンの占有率が算出されるSTEP7。この算出では、第一ループ21で囲まれた面積が算出される。全ての第一ループ21の面積が合計される。この合計の、仮想球14の表面積に対する比率が、占有率である。図5に示された多数の三角形が用いられ、占有率が近似的に算出されてもよい。近似的な算出では、第一ループ21に含まる三角形の合計面積が、全ての三角形の合計面積で除される。得られた占有率に基づき、判定がなされるSTEP8。このステップでは、占有率が所定値以上であるか否かが判定される。図4に示された実施形態では、占有率Yが65以上であるか否かが判定される。占有率Yが65未満である場合、属性の更新がなされるSTEP9。以下、この更新の方法が詳説される。図10は、属性の付与が完了した後のメッシュ12の一部が示された拡大図である。符号16nで示されたセル16は、第一ループ21の上に存在する。このセル16nには、6個のセル16hから16mが隣接している。セル16hのiflagは1であり、その属性はインサイドである。属性がインサイドであるセル16では、iflagの変更はなされない。セル16i、16l、16mのiflagは2であり、その属性は境界である。その属性が境界であり、かつ属性が境界である他のセル16と隣接するセル16では、iflagの変更はなされない。セル16j、16kのiflagは0であり、その属性はアウトサイドである。その属性がアウトサイドであり、かつ属性が境界である他のセル16と隣接するセル16では、iflagが0から3に変更される。第一ループ21の上に存在する全てのセル16に関し、このセル16と隣接するセル16のiflagが決定される。このiflagに基づき、属性の更新STEP9がなされる。属性の更新は、下記のルールに基づいてなされる。iflag:0属性:アウトサイドiflag:1-2属性:インサイドiflag:3属性:境界属性の更新が1回なされた後のメッシュ12は、第二フェーズにある。属性が境界である複数のセル16を結ぶことにより、第二ループ28が得られる。第二ループ28は、第一ループ21の面積以上の面積を有する。換言すれば、属性の更新STEP9により、占有率が大きくなる。多数の第二ループ28を有するパターンが、図11に示されている。図9及び11の対比から明らかなように、図11のパターンの占有率は、図9のそれよりも大きい。このパターンの占有率が算出されるSTEP7。得られた占有率に基づき、判定がなされるSTEP8。このステップでは、占有率が所定値以上であるか否かが判定される。図4に示された実施形態では、占有率Yが65以上であるか否かが判定される。以下同様に、占有率Yが65以上となるまで、属性の更新STEP9、占有率の算出STEP8及び判定STEP9が繰り返される。第N回目の属性更新に先立ち、その属性がアウトサイドであり、かつ属性が境界である他のセル16と隣接するセル16において、iflagが0からN2に変更される。第N回目の属性更新は、下記のルールに基づいてなされる。iflag:0属性:アウトサイドiflag:1からN1まで属性:インサイドiflag:N2属性:境界属性の更新がN回なされた後のメッシュ12は、第N1フェーズにある。2回の属性の更新がなされて得られたパターンが、図12に示されている。このパターンのメッシュ12は、第三フェーズにある。このパターンは、多数の第三ループ29を備えている。第三ループ29は、第二ループ28の面積以上の面積を有する。図9、11及び12の対比から明らかなように、図12に示されたパターンの占有率は大きい。このパターンの占有率は、79である。図13には、1つの第三ループ29が示されている。この第三ループ29は、属性が境界である25個のセル16を結ぶことにより得られている。この第三ループ29は、多数の頂点を有している。図14では、25個のセル16がスプライン曲線で結ばれている。スプライン曲線は、複数の点を通過するスムースな曲線である。スプライン曲線では、隣り合う2つのセル16の間の線が、多項式で定義される。一般的には、3次多項式が用いられる。図13と14との対比から明らかなように、スプライン曲線が用いられることにより、スムースなループが得られる。好ましくは、ループ上のセル16の座標にスムージングがなされ、このセル16に対応する基準点が得られるSTEP10。多数の基準点がスプライン曲線で結ばれることにより、新たなループが想定されるSTEP11。典型的なスムージングは、移動平均である。図15には、3点移動平均によって得られた基準点がスプライン曲線で結ばれて得られたループが示されている。図16には、5点移動平均によって得られた基準点がスプライン曲線で結ばれて得られたループが示されている。図17には、7点移動平均によって得られた基準点がスプライン曲線で結ばれて得られたループが示されている。図14から17の対比より明らかなように、移動平均によって輪郭のスムースが達成されうる。3点移動平均では、下記の3つのセル16の座標が平均される。1当該セル162ループの時計回りにおいて当該セル16に最も近いセル163ループの反時計回りにおいて当該セル16に最も近いセル165点移動平均では、下記の5つのセル16の座標が平均される。1当該セル162ループの時計回りにおいて当該セル16に最も近いセル163ループの反時計回りにおいて当該セル16に最も近いセル164ループの時計回りにおいて当該セル16に2番目に近いセル165ループの反時計回りにおいて当該セル16に2番目に近いセル167点移動平均では、下記の7つのセル16の座標が平均される。1当該セル162ループの時計回りにおいて当該セル16に最も近いセル163ループの反時計回りにおいて当該セル16に最も近いセル164ループの時計回りにおいて当該セル16に2番目に近いセル165ループの反時計回りにおいて当該セル16に2番目に近いセル166ループの時計回りにおいて当該セル16に3番目に近いセル167ループの反時計回りにおいて当該セル16に3番目に近いセル16ループの形成のとき、基準点の一部が間引かれてスプライン曲線が画かれてもよい。図18には、5点移動平均で得られた基準点が半分1点とばしに間引かれて得られたループが示されている。図19には、5点移動平均で得られた基準点が1/32点とばしに間引かれて得られたループが示されている。図19のループを備えたパターンが、図20及び21に示されている。このパターンは、多数のループ30を備えている。ループ30は、仮想球14の表面の上に、ランダムに配置されている。図20及び21に示されたループ30のパターンに基づき、ボロノイ分割により、ディンプルパターンが得られる。ボロノイ分割では、それぞれのループ30の中心点が求められる。中心点の座標は、ループ30の輪郭上のセル及び輪郭の内側に存在するセルの座標の平均値を求めることで得られる。中心点の座標が、ループ30の輪郭の内側に存在するセルのみの座標の平均値を求めることで得られてもよい。中心点の座標が、ループ30の輪郭上に存在するセルのみの座標の平均値を求めることで得られてもよい。以下、この中心点は、母点と称される。図22に、多数の母点32が示されている。ループ30がランダムに配置されているので、母点32も、仮想球14の表面の上にランダムに配置されている。図9に示された第一ループ21に基づいて、母点32が決定されてもよい。この場合も、ランダムに配置された多数の母点32が得られる。図11に示された第二ループ28に基づいて、母点32が決定されてもよい。この場合も、ランダムに配置された多数の母点32が得られる。図12に示された第三ループ29に基づいて、母点32が決定されてもよい。この場合も、ランダムに配置された多数の母点32が得られる。セル16がスプライン曲線で結ばれて得られるループ図14参照に基づいて、母点32が決定されてもよい。この場合も、ランダムに配置された多数の母点32が得られる。スムージングがなされて得られるループ図15-17参照に基づいて、母点32が決定されてもよい。この場合も、ランダムに配置された多数の母点32が得られる。これらの母点32に基づいて、多数のボロノイ領域が想定される。図23には、ボロノイ領域33が示されている。図23において、母点32aは6個の母点32bと隣接している。符号34で示されているのは、母点32aと母点32bとを結ぶ線分である。図23には、6本の線分34が示されている。符号36で示されているのは、それぞれの線分34の垂直二等分線である。母点32aは、6本の垂直二等分線36で囲まれている。図23において白抜き円で示されているのは、垂直二等分線36と他の垂直二等分線36との交点である。この交点は、球面多角形例えば球面六角形の頂点である。この球面多角形が、ボロノイ領域33である。それぞれの垂直二等分線36は、あるボロノイ領域33と、このボロノイ領域33に隣接するボロノイ領域33との境界である。この垂直二等分線36は、ボロノイ領域33の輪郭を画定する。垂直二等分線36の近傍のセル16にはランド10が割り当てられる。他のセル16には、ディンプル8が割り当てられる。これらの割り当てにより、ランド10とディンプル8とからなる凹凸パターンが得られる。垂直二等分線36の近傍のセル16にディンプルが割り当てられ、他のセル16にランドが割り当てられてもよい。垂直二等分線36に基づいてボロノイ領域33の輪郭を画定する計算は、複雑である。以下、簡便にボロノイ領域33が得られる方法が説明される。この方法では、それぞれのセル16における、このセル16と全ての母点32との距離が算出される。セル16毎に、母点32の数と同じ数の距離が算出される。これらの距離の中から、最も短い距離が選定される。この最も短い距離の対象となった母点32に、このセル16が関連づけされる。換言すれば、このセル16に最も近い母点32が選定される。なお、当該セル16からの距離が大きいことが明らかである母点32との距離の計算が、省略されてもよい。それぞれの母点32に関し、この母点32と関連づけされたセル16の集合が想定される。換言すれば、この母点32を最も近い母点32とするセル16の集合が想定される。この集合が、ボロノイ領域33とみなされる。こうして得られた多数のボロノイ領域33が、図24に示されている。図24では、当該セル16と隣接する他のセル16が、当該セル16が属するボロノイ領域33とは異なるボロノイ領域33に属する場合に、当該セルが黒く塗りつぶされている。図24から明らかなように、各ボロノイ領域33の輪郭は、ジグザグである。この輪郭に、セルオートマトン法に関して前述されたスムージング及び間引きがなされ、スプライン曲線が画かれる。このスプライン曲線により、ループが得られる。ループの拡大又は縮小がなされて、新たなループが得られてもよい。このループの外部に、ランド10が割り当てられる。換言すれば、ボロノイ領域33の輪郭の近傍にランド10が割り当てられる。一方、ループの内部にディンプル8が割り当てられる。このようにして、図2及び3に示された凹凸パターンが得られる。ディンプル8の割り当て方法の一例が説明される。この方法では、最深点が決定される。好ましくは、ループの中心と仮想球14の中心とを結ぶ線上に、最深点が想定される。ループの中心の座標は、このループを画定する全ての基準点の座標の平均である。この最深点が、仮想球14の表面に投影される。この投影された点を通過し、その両端がループ上にある、仮想球14の表面の上にある円弧が想定される。この円弧の両端と最深点とを通過し、ゴルフボール2の半径方向内向きに凸な、滑らかな曲線が想定される。好ましくは、滑らかな曲線は、円弧である。この滑らかな曲線とループとが、滑らかな曲面で結ばれる。これにより、ディンプル8が得られる。最深点とループとが滑らかな曲面で結ばれることで、ディンプル8が得られてもよい。図2及び3に示された凹凸パターンでは、多数のディンプル8がランダムに配置されている。それぞれのディンプル8は、他のいずれのディンプル8とも異なる輪郭形状を有している。この凹凸パターンを備えたゴルフボール2は、空力的対称性に優れる。ゴルフボール2のホップが抑制されるとの観点から、ディンプル8の深さは0.05mm以上が好ましく、0.08mm以上がより好ましく、0.10mm以上が特に好ましい。ゴルフボール2のドロップが抑制されるとの観点から、この深さは0.60mm以下が好ましく、0.45mm以下がより好ましく、0.40mm以下が特に好ましい。深さは、ディンプル8の最深点と仮想球14の表面との距離である。本発明においてディンプルの容積とは、仮想球14の表面とディンプル8の表面とに囲まれた部分の容積を意味する。全てのディンプル8の容積の合計総容積は、ゴルフボール2のホップが抑制されるとの観点から400mm3以上が好ましく、450mm3以上がより好ましく、500mm3以上が特に好ましい。ゴルフボール2のドロップが抑制されるとの観点から、この合計は750mm3以下が好ましく、700mm3以下がより好ましく、650mm3以下が特に好ましい。飛行性能の観点から、ディンプル8の面積の合計の、仮想球14の表面積に対する比率占有率は、65以上が好ましく、75以上がより好ましく、80以上が特に好ましい。実質的に球であるというゴルフボール2の本質が損なわれないとの観点から、ディンプル8の総数は250個以上が好ましく、300個以上が特に好ましい。それぞれのディンプル8が十分なディンプル効果を発揮するとの観点から、この総数は450個以下が好ましく、400個以下が特に好ましい。好ましくは、このゴルフボール2の差dRの絶対値は、0.50mm以下である。この絶対値は、ゴルフボール2の空力的対称性と相関するパラメータである。この絶対値が小さいほど、PH回転時の弾道とPOP回転時の弾道との差が小さい。この観点から、絶対値は0.45mm以下がより好ましく、0.40mm以下が特に好ましい。以下、差dRに基づく評価方法が説明される。図25は、この評価方法が説明されるための模式図である。この評価方法では、第一回転軸Ax1が想定される。この第一回転軸Ax1は、ゴルフボール2の2つの極点Poを通過する。それぞれの極点Poは、ゴルフボール2の成形に用いられるモールドの最深点である。一方の極点Poは上型の最深点であり、他方の極点Poは下型の最深点である。ゴルフボール2は、第一回転軸Ax1を中心として回転する。この回転は、PH回転である。このゴルフボール2の仮想球14の表面に存在し、かつ第一回転軸Ax1と直交する大円GCが想定される。ゴルフボール2の回転のとき、この大円GCの周速が最も速い。さらに、ゴルフボール2の仮想球14の表面に存在し、第一回転軸Ax1と直交する2つの小円C1、C2が想定される。図26には、図25のゴルフボール2の一部の断面が模式的に示されている。図26の左右方向は、軸方向である。図26に示されるように、小円C1と大円GCとの中心角の絶対値は、30°である。図示されていないが、小円C2と大円GCとの中心角の絶対値も、30°である。これらの小円C1、C2により上記仮想球14が区画され、ゴルフボール2の表面のうちこれら小円C1、C2に挟まれた領域が特定される。図26における点Pαは、ゴルフボール2の表面に位置し、かつ大円GCとの中心角がα°degreeである点である。点Fαは、点Pαから第一回転軸Ax1に下ろした垂線Peαの足である。矢印L1αで示されているのは、垂線Peαの長さである。換言すれば、長さL1αは、点Pαと第一回転軸Ax1との距離である。1つの断面において、21個の点Pαに関し、長さL1αが算出される。具体的には、-30°、-27°、-24°、-21°、-18°、-15°、-12°、-9°、-6°、-3°、0°、3°、6°、9°、12°、15°、18°、21°、24°、27°及び30°の角度αに関し、長さL1αが算出される。21個の長さL1αが合計され、総長さL2mmが得られる。総長さL2は、図26に示された断面における、表面の形状に依存するパラメータである。図27には、ゴルフボール2の一部の断面が示されている。図27において紙面垂直方向が、軸方向である。図27において符号βで示されているのは、ゴルフボール2の回転角度である。0°以上360°未満の範囲において、0.25°刻みに、回転角度βが設定される。それぞれの回転角度ごとに、総長さL2が算出される。この結果、回転方向に沿って1440の総長さL2が得られる。換言すれば、ゴルフボール2の1回転によって所定箇所に刻々と出現する表面の形状に依存するパラメータに関するデータ群が、算出される。このデータ群は、30240個の長さL1に基づいて算出されたものである。1440の総長さL2の最大値と最小値とが、決定される。この最大値から最小値が減じられ、変動幅Rhが算出される。変動幅Rhは、PH回転における空力特性を表す数値である。さらに、第一回転軸Ax1と直交する第二回転軸Ax2が決定される。第二回転軸Ax2を中心としたゴルフボール2の回転は、POP回転である。PH回転と同様、POP回転についても、大円GCと2つの小円C1、C2が想定される。小円C1と大円GCとの中心角の絶対値は、30°である。小円C2と大円GCとの中心角の絶対値も、30°である。ゴルフボール2の表面のうちこれら小円に挟まれた領域において、1440の総長さL2が算出される。換言すれば、ゴルフボール2の1回転によって所定箇所に刻々と出現する表面の形状に依存するパラメータに関するデータ群が、算出される。1440の総長さL2の最大値と最小値とが決定される。この最大値から最小値が減じられ、変動幅Roが算出される。変動幅Roは、POP回転における空力特性を表す数値である。第一回転軸Ax1と直交する直線は、無数に存在する。従って、大円GCも無数に存在する。ディンプル8に含まれる部分が最長である大円GCが選択され、変動幅Ro及び差dRが算出される。これに代えて、無作為に抽出された20個の大円GCに基づき、20個の変動幅が算出されてもよい。この場合、20個のデータの中の最大値が、Roとされる。変動幅Rhが小さいほど、PH回転時に大きな飛距離が得られうる。その理由は、変動幅Rhが小さいほど、乱流遷移が円滑に継続されるためであると推測される。この観点から、変動幅Rhは3.0mm以下が好ましく、2.7mm以下が特に好ましい。変動幅Roが小さいほど、POP回転時に大きな飛距離が得られうる。その理由は、変動幅Roが小さいほど、乱流遷移が円滑に継続されるためであると推測される。この観点から、変動幅Roは3.0mm以下が好ましく、2.7mm以下が特に好ましい。PH回転時及びPOP回転時のいずれにおいても大きな飛距離が得られうるとの観点から、変動幅Rh及び変動幅Roの両方が3.0mm以下であることが好ましく、2.7mm以下であることが特に好ましい。ボロノイ分割では、ボロノイ領域の1つに母点が1つ含まれるように、仮想球の表面の上に多数のボロノイ領域が想定される。ボロノイ分割以外の方法によって、領域の1つに点が1つ含まれるように、仮想球の表面の上に多数の領域が想定されてもよい。図28は、本発明の他の実施形態に係るゴルフボール102が示された拡大正面図である。図29は、図28のゴルフボール102が示された平面図である。図28び29から明らかなように、多数のディンプル108がランダムに配置されている。これらディンプル108とランド110とにより、ゴルフボール102の表面に凹凸パターンが形成されている。この凹凸パターンの設計方法には、ドロネー三角分割Delaunaytriangulationが用いられる。この設計方法では、多数の母点が仮想球の表面の上にランダムに配置される。この母点に基づいて、ドロネー三角分割により、仮想球の表面の上に多数の領域が想定される。本願明細書では、この領域は、ドロネー領域と称される。このドロネー領域に基づいて、ディンプル108が割り当てられる。好ましくは、母点の配置には、セルオートマトンCellularAutomaton法が用いられる。図2及び3に示されたディンプルパターンの設計に関して前述されたものと同様のセルオートマトンが、用いられうる。セルオートマトン法により、仮想球の表面の上に多数のループがランダムに配置されたパターンが得られる。これらループの中心点が、母点とされる。ループの配置がランダムなので、母点の配置もランダムである。セルオートマトン法によって得られた多数の第一ループ121を有するパターンが、図30に示されている。このパターンは、以下のパラメータが用いられて得られたものである。W1:1.0W2:-0.6R1:6.5R2:11.3セルの数:88266図30のパターンが更新されて、占有率が高められる。第一回の更新がなされた後のパターンが、図31に示されている。このパターンは、多数の第二ループ128を備えている。第二回の更新がなされた後のパターンが、図32に示されている。このパターンは、多数の第三ループ129を備えている。この第三ループ129に、前述のスムージングが施される。スムージング後のループの基準点が、スプライン曲線で結ばれて、新たなループが想定される。このループを備えたパターンが、図33及び34に示されている。図33及び34に示されたループ130のパターンに基づき、ドロネー三角分割により、ディンプルパターンが得られる。ドロネー三角分割では、それぞれのループ130の中心点が求められる。中心点の座標は、ループ130の輪郭上のセル及び輪郭の内側に存在するセルの座標の平均値を求めることで得られる。中心点の座標が、ループ130の輪郭の内側に存在するセルのみの座標の平均値を求めることで得られてもよい。中心点の座標が、ループ130の輪郭上に存在するセルのみの座標の平均値を求めることで得られてもよい。以下、この中心点は、母点と称される。図35に、多数の母点132が示されている。ループ130がランダムに配置されているので、母点132も、仮想球の表面の上にランダムに配置されている。図30に示された第一ループ121に基づいて、母点が決定されてもよい。この場合も、ランダムに配置された多数の母点が得られる。図31に示された第二ループ28に基づいて、母点が決定されてもよい。この場合も、ランダムに配置された多数の母点が得られる。図32に示された第三ループ129に基づいて、母点が決定されてもよい。この場合も、ランダムに配置された多数の母点が得られる。これらの母点132に基づいて、ドロネー三角分割により、仮想球の表面が多数の三角形ドロネー領域に分割される。以下、この分割の方法が説明される。この方法では、多数の母点132から、任意の3つの母点が選択される。これら3つの母点132を頂点とする三角形が、想定される。この三角形の外接円が、想定される。これら3つの母点以外の母点が、この外接円の中に含まれないとき、この三角形はドロネー領域と判定される。3つの母点の組み合わせの全てについて、判定がなされる。これにより、仮想球の表面の全体が、多数のドロネー領域に分割される。仮想球の表面のそれぞれのセルが属するドロネー領域に属するかが、決定される。以下、この決定方法が説明される。図36aに示された三角形P1P2P3ドロネー領域において、ベクトルP1P2、ベクトルP2P3及びベクトルP3P1は、反時計回りに連続している。従って、この三角形P1P2P3の、右手座標系における法線ベクトルは、紙面に垂直でかつ上向きである。図36bに示された三角形P4P5P6ドロネー領域において、ベクトルP4P5、ベクトルP5P6及びベクトルP6P4は、時計回りに連続している。従って、この三角形P4P5P6の、右手座標系における法線ベクトルは、紙面に垂直でかつ下向きである。法線ベクトルが下向きである三角形P4P5P6では、第二頂点P5と第三頂点P6とが、入れ替えられる。換言すれば、三角形P4P5P6は、三角形P4P6P5と再定義される。法線ベクトルが下向きである全ての三角形について、再定義がなされる。この再定義がなされた後の仮想球の表面では、全ての三角形の法線ベクトルは、上向きである。図36cには、三角形P1P2P3とセルS1とが示されている。この例では、3つの外積ベクトルP1P2×P1S1、P2P3×P2S1及びP3P1×P3S1は、紙面に対して上向きである。3つの外積ベクトルの、紙面に対する方向が揃っているとき、当該セルS1は、三角形P1P2P3に含まれる。図36dには、三角形P1P2P3とセルS2とが示されている。この例では、2つの外積ベクトルP1P2×P1S2及びP2P3×P2S2は、紙面に対して上向きである。一方、外積ベクトルP3P1×P3S2は、紙面に対して下向きである。3つの外積ベクトルの、紙面に対する方向が、揃っていないとき、当該セルS2は、三角形P1P2P3に含まれない。全てのセルについて、属するドロネー領域が決定される。図37では、当該セル16と隣接する他のセル16が、当該セル16が属するドロネー領域133とは異なるドロネー領域133に属する場合に、当該セルが黒く塗りつぶされている。図37から明らかなように、各ドロネー領域133の輪郭は、ジグザグである。この輪郭に、セルオートマトン法に関して前述されたスムージング及び間引きがなされ、スプライン曲線が画かれる。このスプライン曲線により、ループが得られる。ループの拡大又は縮小がなされて、新たなループが得られてもよい。このループの外部に、ランド110が割り当てられる。換言すれば、ドロネー領域133の輪郭の近傍にランド110が割り当てられる。一方、ループの内部にディンプル108が割り当てられる。このようにして、図28及び29に示された凹凸パターンが得られる。図28及び29に示された凹凸パターンでは、多数のディンプル108がランダムに配置されている。それぞれのディンプル108は、他のいずれのディンプル108とも異なる輪郭形状を有している。この凹凸パターンを備えたゴルフボール102は、空力的対称性に優れる。ゴルフボール102のホップが抑制されるとの観点から、ディンプル108の深さは0.05mm以上が好ましく、0.08mm以上がより好ましく、0.10mm以上が特に好ましい。ゴルフボール102のドロップが抑制されるとの観点から、この深さは0.60mm以下が好ましく、0.45mm以下がより好ましく、0.40mm以下が特に好ましい。深さは、ディンプル108の最深点と仮想球の表面との距離である。全てのディンプル108の容積の合計総容積は、ゴルフボール102のホップが抑制されるとの観点から400mm3以上が好ましく、450mm3以上がより好ましく、500mm3以上が特に好ましい。ゴルフボール102のドロップが抑制されるとの観点から、この合計は750mm3以下が好ましく、700mm3以下がより好ましく、650mm3以下が特に好ましい。飛行性能の観点から、ディンプル108の面積の合計の、仮想球の表面積に対する比率占有率は、65以上が好ましく、75以上がより好ましく、80以上が特に好ましい。実質的に球であるというゴルフボール102の本質が損なわれないとの観点から、ディンプル108の総数は250個以上が好ましく、300個以上が特に好ましい。それぞれのディンプル108が十分なディンプル効果を発揮するとの観点から、この総数は450個以下が好ましく、400個以下が特に好ましい。好ましくは、このゴルフボール102の差dRの絶対値は、0.50mm以下である。この絶対値が小さいほど、PH回転時の弾道とPOP回転時の弾道との差が小さい。この観点から、絶対値は0.45mm以下がより好ましく、0.40mm以下が特に好ましい。変動幅Rhが小さいほど、PH回転時に大きな飛距離が得られうる。その理由は、変動幅Rhが小さいほど、乱流遷移が円滑に継続されるためであると推測される。この観点から、変動幅Rhは3.0mm以下が好ましく、2.7mm以下が特に好ましい。変動幅Roが小さいほど、POP回転時に大きな飛距離が得られうる。その理由は、変動幅Roが小さいほど、乱流遷移が円滑に継続されるためであると推測される。この観点から、変動幅Roは3.0mm以下が好ましく、2.7mm以下が特に好ましい。PH回転時及びPOP回転時のいずれにおいても大きな飛距離が得られうるとの観点から、変動幅Rh及び変動幅Roの両方が3.0mm以下であることが好ましく、2.7mm以下であることが特に好ましい。ドロネー三角分割では、1つの母点と他の母点とを線で結ぶことにより、仮想球の表面の上に多数のドロネー領域が想定される。ドロネー三角分割以外の方法であって、1つの点と他の点とを線で結ぶ方法により、仮想球の表面の上に多数の領域が想定されてもよい。[評価1]図2及び3に示された、実施例1のパターンを設計した。このパターンは、391個のディンプルを有している。さらに、図20及び21に示された、比較例1のパターンを設計した。このパターンは、391個のディンプルを有している。前述の方法により、それぞれのパターンの変動幅Ro及びRhを算出した。この結果が、下記の表1に示されている。表1に示されるように、実施例1のパターンの変動幅Ro及びRhは小さい。この評価結果から、本発明の優位性は明らかである。[評価2]図28及び29に示された、実施例2のパターンを設計した。このパターンは、312個のディンプルを有している。さらに、図33及び34に示された、比較例2のパターンを設計した。このパターンは、158個のディンプルを有している。さらに、図46及び47に示された、比較例3のパターンを設計した。このパターンは、392個のディンプルを有している。前述の方法により、それぞれのパターンの変動幅Ro及びRhを算出した。この結果が、下記の表2に示されている。表2に示されるように、実施例のパターンの変動幅Ro及びRhは小さい。この評価結果から、本発明の優位性は明らかである。"